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ABSTRACT 


Sugarbeet  cylinders  were  subjected  to  a  number  of  different 
freezing  and  thawing  conditions.  Each  set  of  conditions 
resulted  in  a  marked  softening  of  the  beet  tissue.  The 
thawed  beets  decayed  rapidly.  Consequently,  if  beets  are 
stored  frozen,  they  should  not  be  allowed  to  thaw  prior  to 
processing . 

In  order  to  determine  the  effect  of  freezing  rate  on 
sugar  loss  in  storage,  samples  were  frozen  at  different 
rates  before  being  stored  at  263  K  or  248  K.  Experimental 
results  indicated  that  all  freezing  methods  produced  a 
negligible  change  in  sugar  content  after  six  months  storage 
at  either  of  the  storage  temperatures. 

A  numerical  model  was  developed  for  the  freezing  or 
thawing  of  an  infinite  cylinder  of  biological  material.  This 
model  predicted  the  variation  of  temperature  with  time  and 
with  radial  position.  The  model  predictions  were  in  good 
agreement  with  experimental  values. 
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1.  INTRODUCTION 


Sugar  is  an  important  commercial  product.  It  may  reach  the 
consumer  in  its  pure  form  or  as  a  component  of  various  food 
products  such  as  soft  drinks,  baked  foods  and  ice  cream.  It 
is  also  used  in  the  production  of  non-food  items  such  as 
pharmaceuticals  and  explosives.  Alberta's  only  sugar  factory 
is  located  at  Taber.  In  this  factory  sugar  is  extracted  from 
sugarbeets  which  contain  about  16%  sugar.  Approximately  90% 
of  this  sugar  is  recovered  as  bagged  sugar.  However,  with 
poor  beets  the  recovery  may  be  as  low  as  70%.  By-products  of 
the  factory  are  cattlefeed,  molasses,  and  invert  sugar.  The 
factory  operates  continuously  from  the  beginning  of  October 
to  the  end  of  January.  About  450,000  tonnes  of  beets  are 
processed  per  year  to  produce: 

1.  60,000  tonnes  of  sugar, 

2.  35,000  tonnes  of  beet  pulp,  which  is  a  cattle  feed,  and 

3.  20,000  tonnes  of  molasses,  which  is  used  in  the 
production  of  cattle  feed  or  yeast. 

At  harvest  time  the  leaves  and  tops  are  removed  from 
the  beets,  which  are  then  placed  in  large  storage  piles  to 
await  processing.  During  storage,  sugar  is  lost,  mainly  as  a 
result  of  respiration.  It  has  been  estimated  (229)  that  10% 
of  beet  sugar  in  the  U.S.  is  lost  during  storage.  This  loss 
may  be  reduced  by  lowering  the  storage  temperature.  Frozen 
storage  is  being  considered  in  the  U.S.,  and  Russia,  as  a 
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means  of  reducing  or  eliminating  sugar  loss.  This  method  of 
storage  becomes  economically  attractive  if  cold  outdoor  air 
can  be  used  to  freeze  the  beets.  However,  rapid 
deterioration  of  the  beets  may  occur  if  they  are  allowed  to 
thaw  in  the  storage  piles.  Deterioration  results  from  the 
injury  suffered  during  freezing  and  thawing.  It  may  be 
possible  to  minimize  this  damage  by  a  careful  evaluation  of 
freezing  and  thawing  conditions.  Texture  measurements  are 
useful  in  this  evaluation.  Freshly  harvested  beets  have  a 
firm  texture  whereas  freeze-damaged  beets  are  soft.  A  beet 
which  has  been  softened  by  freezing  and  thawing  is  probably 
damaged  and  is  likely  to  deteriorate  rapidly. 

Heat  transfer  models  can  also  be  used  to  estimate  the 
severity  of  the  freezing  injury  which  is  suffered  by  a 
biological  sample  when  it  is  subjected  to  specific  freezing 
and  thawing  conditions.  Freezing  injury  may  depend  on  the 
rate  at  which  a  sample  is  frozen  and  thawed.  Different 
regions  of  large  samples  may  experience  different  rates  of 
temperature  change  and,  consequently,  sustain  different 
degrees  of  injury.  To  obtain  an  overall  estimate  of  damage, 
it  is  useful  to  predict  how  quickly  the  temperature  changes 
throughout  the  sample  during  freezing  and  thawing. 
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1.0.1  Objectives  of  this  Thes i s 

In  this  thesis  it  was  desired  to: 

1.  examine  a  number  of  different  freezing  and  thawing 
conditions  in  order  to  produce  a  thawed  beet  with  a  firm 
texture , 

2.  determine  if  sucrose  losses  in  frozen  storage  depend  on 
the  freezing  conditions  to  which  a  beet  is  subjected 
prior  to  storage,  and 

3.  develop  a  numerical  model  which  would  predict  the 
transient  temperature  profile  in  an  infinite  cylinder  of 
biological  material  during  freezing  or  thawing. 


2.  PRODUCTION  OF  SUGAR  FROM  SUGARBEETS 
The  production  of  sugar  from  sugarbeets  is  briefly  described 
in  this  chapter.  This  description  is  intended  as  an  aid  in 
understanding  processing  problems,  which  arise  when 
sugarbeets  are  stored  incorrectly. 

Sugarbeets  are  removed  from  storage  piles  and  placed  in 
flumes  outside  the  factory.  Water  then  transports  the  beets 
into  the  factory.  The  beets  are  washed  with  water  in  an 
agitated  vessel.  After  washing  they  are  cut  into  thin  slices 
( "cossettes" ) .  Sugar  is  leached  from  the  cossettes  in  a 
diffuser  by  hot  water  which  flows  countercurrent ly .  The 
depleted  cossettes  ("pulp")  are  dried,  combined  with 
molasses,  and  pelletized  to  produce  cattle  feed.  The  extract 
solution  ("raw  juice")  proceeds  to  the  carbonation  stage. 

In  carbonation,  milk  of  lime  (a  mixture  of  calcium 
oxide  and  water)  and  carbon  dioxide  are  added  to  the  raw 
juice  to  precipitate  impurities,  such  as  organic  acids  and 
inorganic  salts.  The  precipitate  is  removed  by  pumping  the 
juice  through  filter  presses.  Sulfur  dioxide  is  then  added 
to  prevent  discolouration  of  the  finished  product.  The 
resulting  solution  is  "thin  juice".  This  juice  is 
concentrated  in  multiple  effect  evaporators  to  produce 
"thick  juice",  which  is  about  65%  solids. 

"Standard  liquor"  is  obtained  by  adding  recycled  sugar 
to  the  thick  juice.  The  liquor  is  pressure  filtered  (a 
filter  aid  is  used)  and  boiled  under  vacuum  in  the  "White 
Pan"  until  sugar  crystals  are  formed.  The  sugar  crystals 
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must  now  be  separated  from  the  syrup  in  which  they  are 
suspended.  This  is  accomplished  in  the  '’spinners1',  which  are 
large  circular  baskets  rotated  at  high  speed.  The  syrup  is 
forced  through  the  basket  by  centrifugal  action.  The  sugar 
remains  on  the  basket  and  is  washed  with  hot  water.  It  is 
then  dried  and  packaged.  The  syrup  is  processed  to  produce 
molasses  and  low  grade  sugars.  The  molasses  is  sold  as  a 
finished  product  or  is  added  to  pulp  to  form  cattle  feed. 

The  low  grade  sugars  are  dissolved  in  thin  juice  and  added 
to  the  thick  juice  to  form  standard  liquor. 
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3.  STORAGE  OF  UNFROZEN  SUGARBEETS 


When  beets  are  stored,  sucrose  is  lost  as  a  result  of 
aerobic  and  anaerobic  respiration,  biochemical 
transformations  and  attack  by  micro-organisms.  Wyse  (227) 
has  defined  ’’respiration  in  storage"  to  be  "the  process 
whereby  the  root  converts  sucrose  into  energy  to  maintain 
its  physiological  integrity".  In  aerobic  respiration, 
sucrose  is  degraded  according  to  the  exothermic  reaction: 

CxaHajOn  +  1 20  2  -->  12C0  2  +  11H20 
( sucrose ) 

Under  ideal  storage  conditions  aerobic  respiration  is 
responsible  for  60  to  75%  of  the  sucrose  loss  (213). 

Anaerobic  respiration  may  occur  in  regions  of  the 
storage  pile  where  air  ventilation  is  blocked  by  dirt, 
trash,  or  snow.  It  can  be  represented  by: 

C12H  2  201  2  +  H20  -->  4CO 2  +  4C2H5OH 

(ethanol ) 

This  form  of  respiration  can  rapidly  lead  to  cell  death  and 
invasion  by  micro-organisms. 

Micro-organisms  present  in  beet  tissue  use  sucrose  for 
their  respiration  and  growth  (88).  They  also  deposit  enzymes 
in  the  beet  (138).  These  enzymes  catalyze  reactions  which 
convert  sucrose  into  glucose,  fructose,  and  other 
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carbohydrates.  Damaged  beets  are  particularly  susceptible  to 
attack  by  micro-organisms. 

Enzymes  naturally  present  in  the  beet  tissue  are  also 
responsible  for  some  sucrose  loss.  One  enzyme,  invertase, 
catalyzes  the  hydrolysis  of  sucrose  to  glucose  and  fructose. 
This  mixture  of  glucose  and  fructose  is  known  as  invert 
sugar  because  it  rotates  plane  polarized  light  in  the 
opposite  direction  to  sucrose. 

Raffinose,  a  tr i sacchar ide ,  may  also  accumulate  during 
storage,  at  the  expense  of  sucrose.  The  mechanism  by  which 
raffinose  is  formed  is  still  unknown  (13)  but  it  is  thought 
that  it  results  from  the  combination  of  sucrose  and  a 
compound  (UDP  galactose)  formed  from  glucose. 

Nonsucrose  components  ("impurities")  of  sugarbeets 
cause  problems  in  the  sugar  factory.  Invert  sugars  are 
degraded  to  organic  acids,  such  as  formic,  acetic,  lactic, 
and  glucinic  acid.  These  acids  tend  to  increase  colour 
formation  (34)  and  reduce  the  pH  of  factory  sugar  solutions. 
This  decrease  in  pH  results  in  an  increase  in  the  amount  of 
lime  which  must  be  added  in  carbonation.  The  concentration 
of  lime  in  the  thin  juice  is  then  increased.  This  may  result 
in  an  increase  in  the  quantity  of  lime  which  is  deposited  in 
the  evaporators,  causing  scale  formation.  A  high  raffinose 
content  will  decrease  crystallization  rates  and  produce 
distorted  and  elongated  sugar  crystals.  Impurities,  in 
general,  tend  to  increase  the  amount  of  sugar  contained  in 
molasses  and  decrease  the  quantity  of  bagged  sugar  obtained. 
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3 . 1  Factors  Af f ect inq  Sugar  Loss  in  Storage 

3.1.1  Temperature 

The  quantity  of  sugar  lost  and  impurities  accumulated 
in  storage  will  depend  on  a  number  of  variables,  of  which 
temperature  is  probably  the  most  important.  The  sucrose  loss 
will  approximately  double  for  every  8  K  rise  in  temperature 
(48)  above  273  K.  The  invert  sugar  content  of  beets  also 
increases  as  the  temperature  is  increased.  On  the  other 
hand,  the  raffinose  content  increases  as  the  storage 
temperature  is  reduced  below  about  278  K. 

If  the  temperature,  anywhere  in  the  storage  pile, 
reaches  319  K  then  the  beet  structure  will  be  destroyed 
within  a  week  (156).  Micro-organisms  can  then  attack  the 
dead  cells.  This  results  in  a  further  rise  in  temperature 
and  a  rapid  loss  of  sugar. 


3.1.2  Time  in  Storage 

As  the  storage  time  is  increased  so  also  will  be  the 
sugar  loss.  However,  there  is  nonlinear  relationship  between 
sugar  loss  and  storage  time.  A  high  respiration  rate  is 
obtained  immediately  after  harvest.  This  may  result  from: 


1. 

2. 

3. 


injuries  sustained  during  the  harvest, 
metabolic  changes  resulting  from  the  leaves  being 
removed,  or 

the  large  exposed  surface  area,  which  is  available  for 
oxygen  diffusion,  before  the  surface  becomes  hardened  by 
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dehydration . 


Within  a  week  of  storage,  the  respiration  rate  drops  to 
approximately  half  its  initial  value  (227).  After  a  few 
months,  respiration  may  increase  again  as  a  result  of  mould 
growth  ( 229 )  . 

3,1.3  Condition  of  Beet  Piles 

Since  respiration  is  an  exothermic  reaction,  whose  rate 
increases  with  temperature,  heat  must  be  removed  or  high 
temperatures  will  be  attained  in  the  storage  pile.  The 
presence  of  dirt  and  trash  in  piles  impedes  air  circulation 
and,  consequently,  the  removal  of  heat.  Local  overheating  in 
the  pile  may  then  occur.  Beets  may  be  killed  by  the  high 
temperatures  obtained  or  by  a  deficiency  in  oxygen  supply. 
Microbial  fermentation  of  the  soil  and  trash  may  also  cause 
overheating  by  as  much  as  30  K. 


3.1.4  Condi t i on  of  Beets  Enter  in 
The  condition  of  the  beets 
also  important.  Beets  suffering 
Violet  Root  Rot,  Clamp  Rot,  and 
rapidly  during  storage  and  will 
diseases  are  caused  by  viruses, 
is  Phoema  betae.  Wilted  beets  lo 


£  Storage 

entering  the  storage  pile  is 
from  diseases,  such  as 
Crown  Rot,  will  decay 
infect  other  beets.  These 
the  most  prevalent  of  which 
se  sugar  quickly  and  are 


likely  to  undergo  fermentation.  Topped  beets  (crown  removed) 
show  greater  spoilage  and  sucrose  losses  than  untopped  beets 


ei  4  .10-  s  -:i  >•  :  1  :A  r  *  ' 


■ 


10 


(2).  Despite  this,  beets  are  generally  topped  because  the 
crown  contains  a  high  percentage  of  impurities.  Zielke  and 
Snyder  (237)  reported  that  the  concentration  of  impurities 
was  70%  higher  in  the  crown  than  in  the  remainder  of  the 
root.  If  leaf  buds  are  not  properly  removed  from  the  beet, 
then  shoots  will  grow  at  the  expense  of  sugar  during 
storage . 

3.1.5  Use  of  Pile  Covers 

Numerous  covers  have  been  used  for  storage  piles.  These 
include  storage  sheds,  covered  trenches,  straw,  beet  tops, 
earth,  plastic  canopies,  and  air  inflated  bubble  structures. 
Covers  have  the  advantage  of  minimizing  dehydration  and 
repeated  freezing  and  thawing  of  the  beets.  While  frozen 
storage  reduces  sugar  losses,  alternate  freezing  and  thawing 
increase  them  (58).  Covers  have  the  disadvantage  of  making 
it  more  difficult  to  cool  beet  piles  during  warm  weather. 

The  relative  merits  of  the  different  covers  will  depend  on 
the  weather  conditions  at  the  given  location. 

3.1.6  Controlled  Atmosphere  Storage 

In  controlled  atmosphere  storage,  beets  are  stored  in 
an  atmosphere  which  contains  about  5%  of  both  carbon  dioxide 
and  oxygen.  This  type  of  storage  reduces  both  the 
respiration  rate  and  the  rate  of  formation  of  impurities. 

The  oxygen  content  should  not  be  lower  than  5%  as  anaerobic 
respiration  will  then  occur  (45,227).  Despite  its 
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advantages,  controlled  atmosphere  storage  of  sugarbeets  has 
been  restricted  to  research  studies  (109,226,227). 

3.1.7  Use  of  Vent i lat i on 

Pile  ventilation  with  cold  air  has  been  used  to  lower 
the  temperature  in  the  pile  and  to  prevent  localized 
overheating.  If  ventilation  is  used,  humidifiers  may  be 
necessary  to  prevent  beet  dehydration. 

3.1.8  Grow i nq  Condi t ions 

Growing  conditions  influence  the  storage  behaviour  of 
sugarbeets.  For  example,  beets  generally  store  better  after 
a  wet  season  of  good  growth  than  after  a  dry  season  (104). 
Also,  beets  store  badly  if  large  quantities  of  nitrogen 
fertilizer  are  added  just  before  the  harvest  (104). 

3.1.9  Beet  Variety 

Nelson  and  Wood  (153)  observed  that  the  quantity  of 
sugar  lost,  during  storage,  was  dependent  upon  the  variety 
of  beet  which  was  stored. 

The  storage  and  processing  of  frozen  beets  is 
considered  in  the  next  chapter. 


4 


4.  STORAGE  AND  PROCESSING  OF  FROZEN  SUGARBEETS 
It  is  estimated  that  about  0.25  kg  of  sucrose/tonne  of 
beets/day  may  be  lost  in  storage  piles  in  the  Red  River 
Valley  of  North  Dakota  and  Minnesota  (15).  Based  on  a 
refined  sugar  price  of  $l/kg,  the  elimination  of  sugar  loss 
in  storage  piles  results  in  potentially  increased  revenues 
of  $0. 25/tonne  of  beets/day  stored.  If  this  figure  is 
applied  to  the  factory  at  Taber,  then  elimination  of  sugar 
loss  will  increase  operating  revenues  by  about  $7  million. 
This  is  approximately  10%  of  their  expected  gross  revenues 
for  the  1980-81  season. 

Frozen  storage  of  sugarbeets  has  been  observed  to 
drastically  reduce  the  amount  of  sugar  lost  during  storage 
(15,58,77,201,213,230).  It  has  also  been  concluded  (201) 
that  it  is  "easier”  to  process  a  frozen  beet,  in  March,  than 
one  which  has  been  stored  unfrozen.  The  decision  of  whether 
or  not  to  use  frozen  sugarbeet  storage  will,  of  course, 
depend  on  economic  considerations.  The  cost  of  frozen 
storage  must  be  balanced  against  the  extra  revenue  obtained 
from  a  reduction  in  sucrose  losses.  Frozen  storage  is 
economically  attractive  in  regions  where  cold  outdoor  air 
can  be  used  to  freeze  the  beets.  Use  of  expensive  freezing 
equipment  is  then  avoided. 
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4 . 1  Frozen  Storage  in  the  U . S . S . R, 

As  the  U.S.S.R.  is  the  world's  largest  sugarbeet 
producer  and  has  an  ideal  climate  for  freezing  beets 
outdoors,  it  is  not  surprising  that  frozen  storage  of 
sugarbeets  has  been  investigated  there.  In  1974,  Moroghan 
(150)  reported  that  about  200,000  tonnes  of  frozen  beets 
were  being  stored  in  the  U.S.S.R.  Forced  ventilation  of  cold 
outdoor  air  (at  a  temperature  less  than  261  K)  is  used  to 
freeze  the  beets  to  a  temperature  below  268  K,  which  is 
considered  to  be  a  "thawing  temperature".  The  air  is 
saturated  with  moisture  to  prevent  dehydration  of  the  beets. 
The  piles  are  covered  with  a  layer  of  ice  or  snow  which 
forms  a  heat  transfer  barrier  to  the  surroundings.  Ice 
covers  are  simply  formed  by  spraying  cold  water  on  top  of 
the  pile.  The  ice  layer  is  melted,  when  required,  by 
spraying  hot  water  on  the  pile. 


4 . 2  Frozen  Storage  i n  the  Red  River  Valley 

In  1974-75,  the  American  Crystal  Sugar  Company  (201) 
examined  two  methods  of  freezing  storage  piles.  One  pile  was 
frozen  by  blowing  cold  air  through  ducts  (metal  drums) 
underneath  the  pile.  A  second  pile,  the  "split  pile",  was 
formed  by  removing  beets  from  the  centre  of  a  regular  pile. 
It  was  then  allowed  to  freeze  naturally.  After  4  months  of 
storage,  each  pile  had  lost  an  average  of  0.1  kg  of 
sugar/tonne  of  beets/day  stored.  This  was  an  improvement  on 
unfrozen  piles  where  0.22  kg  of  sucrose  were  lost/tonne  of 
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beets/day  stored. 
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4 . 3  Problems  with  Frozen  Storage 

The  main  problem  with  frozen  storage  of  sugarbeets  is 
that  thawing  may  occur.  This  is  an  important  disadvantage  as 
thawed  beets  become  unprocessable  within  a  week  (157). 
Sugarbeet  cells  are  very  susceptible  to  freezing  injury,  a 
phenomenon  discussed  in  Appendix  B.  Dexter  et  al.  (59) 
stated  that  some  thawed  beets  "disintegrate  completely" 
during  fluming  and  washing.  Hull  (104)  observed  that  "when 
frozen  beets  begin  to  decay  they  heat  very  quickly,  have  a 
pungent,  vinegary  odour,  and  rapidly  collapse  into  a  putrid, 
slimy  mess". 

Some  specific  problems  are  encountered  with  thawed 
beets.  Pectic  substances  are  complex  carbohydrates,  which 
are  present  in  cell  walls.  They  occur  in  two  forms:  (1) 
water  soluble  pectins  and  (2)  water  insoluble  protopectin. 
The  soluble  pectins  form  calcium  pectate  gels  in  the  sugar 
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factory  (232)  and  these  gels  reduce  filtration  rates. 

Enzymes  degrade  the  pectic  substances  of  thawed  beets  and 
increase  the  percentage  of  soluble  pectins.  Consequently, 
thawed  beets  have  a  high  percentage  of  soluble  pectins.  They 
can,  therefore,  be  expected  to  produce  juices  which  are  more 
difficult  to  filter  than  those  produced  by  unfrozen  beets. 

Oldfield  et  al.  (156)  observed  that  thawed  beets, 
stored  at  293  K,  had  respiration  rates  (as  measured  by  the 
amount  of  carbon  dioxide  evolved)  more  than  ten  times  the 
rate  for  unfrozen  beets.  Most  of  the  sucrose  loss  was  due  to 
respiration  and  very  little  was  due  to  inversion. 

Micro-organisms  are  responsible  for  much  of  the 
deterioration  in  thawed  beets.  Damaged  cell  membranes  allow 
the  exodus  of  the  cell  solution  into  the  intercellular 
spaces  where  micro-organisms  are  present.  It  is  not  known 
whether  these  micro-organisms  enter  from  soil  after  topping 
or  are  naturally  present  in  the  beet.  Bacteria  (Leuconostoc 
mesenteroides )  act  on  the  cellular  solution  to  produce 
polysaccharide  gums  (dextran  and  levan),  which  cause 
filtration  problems  in  the  processing  plant  (13).  The 
polysaccharide  gums  lower  filtration  rates  because  they 
(particularly  dextran)  reduce  the  size  of  precipitated 
particles  and  increase  the  viscosity  of  the  juice,  which  is 
to  be  filtered.  Measurable  quantities  of  dextran  or  levan 
are  not  found  in  fresh,  mould  damaged,  heat  damaged  or 
waterlogged  beet  (157).  Consequently,  the  presence  of  these 
gums  is  taken  to  indicate  that  frost  injury  has  occurred. 
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Yeasts  also  act  on  the  cell  solution  to  produce  est 
the  characteristic  fruity  odour  of  thawed  beets  (13 
After  frozen  beets  have  been  thawed,  their  inv 
content  increases  rapidly  (157).  The  production  of 
sugars  starts  soon  after  thawing,  whereas  there  is 
delay  before  gums  are  formed.  The  lower  the  tempera 
reached  during  freezing,  the  faster  will  be  the  inv 
formation  on  thawing  (157).  The  increase  in  invert 
content  is  important  because: 


ers  and 
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1.  It  occurs  at  the  expense  of  sucrose. 

2.  The  invert  sugar  is  degraded  in  the  sugar  factory, 
during  carbonation  and  evaporation,  to  produce  organic 
acids.  Increased  amounts  of  lime  must  then  be  added, 
during  carbonation,  to  remove  these  acids. 

3.  The  higher  the  invert  sugar  content  of  a  beet  the  lower 
is  the  fraction  of  its  sucrose  which  is  recovered  as 
bagged  sugar  in  the  sugar  factory. 

4.  Degradation  of  invert  sugars  is  accompanied  by  increased 
colour  formation  in  factory  sugar  solutions. 


Other  disadvantages  of  thawed  beets  are: 

1.  They  may  lose  as  much  as  6%  of  their  sugar  content  into 
the  factory  transportation  and  cleaning  water  (196). 

2.  They  may  rot  in  the  piles,  causing  the  formation  of 
spots" . 
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3.  They  have  a  soft  texture  and  are,  therefore,  difficult 
to  slice. 

Because  of  the  difficulties  encountered  with  thawed 
beets,  frozen  beets  are  generally  processed  without  being 
allowed  to  thaw. 


4 . 4  Processing  Frozen  Sugarbeets 

Frozen  sugarbeets  are  more  difficult  to  process  than 
unfrozen  beets.  Additional  power  must  be  supplied  to  the 
slicers  (201).  Wider  knives,  producing  thicker  slices,  are 
required.  Poorly  formed  cossettes,  which  plug  the  diffuser, 
are  obtained.  Slicer  knives  must  be  replaced  often,  as  they 
can  be  damaged  by  dirt  and  stones,  which  are  frozen  to  the 
surface  of  the  beets. 

Extra  heating  is  needed  to  bring  the  frozen  cossettes 
to  the  diffuser  operating  temperature.  This  extra  heating  is 
not  as  much  as  would  at  first  be  expected  as  low  diffuser 
temperatures  can  be  used.  This  is  because  sucrose  enters  the 
diffuser  water  more  easily  from  freeze-damaged  than  from 
healthy  cells. 

When  frozen  beets  reach  the  hoppers,  above  the  slicers, 
their  surface  is  wet  as  a  result  of  washing.  If  there  is  any 
hold-up  in  the  sugar  factory,  this  water  will  freeze,  as  a 
result  of  losing  heat  to  the  frozen  beet.  Freezing  of  the 
surface  water  causes  the  beets  in  the  hopper  to  stick 
together.  Consequently,  when  frozen  beets  are  being 


_ 

. 


18 


processed,  fewer  beets  are  allowed  to  build  up  in  the 
hopper.  This  ensures  faster  movement  of  beets  and 
facilitates  the  separation  of  the  beets  if  they  do  stick 
together . 

Overall,  the  processing  of  frozen  beets  is  not 
considered  to  be  a  major  problem.  In  Taber,  the  throughput 
of  sugarbeets  only  falls  about  5%  when  frozen  beets  are 
used.  Swift  (201)  stated  that  the  greatest  problems  are 
encountered  when  mixtures  of  frozen  and  unfrozen  beets  must 
be  processed. 


5.  TEXTURE  EXPERIMENTS 


Following  the  frozen  storage  of  sugarbeets,  the  beets  may  be 
processed  in  either  the  frozen  or  the  thawed  state.  If  the 
beets  are  thawed,  prior  to  processing,  then  it  is  essential 
that  freeze-thaw  damage  be  kept  to  a  minimum.  Otherwise, 
rapid  deterioration  will  occur  after  thawing.  Beet  texture 
is  a  good  indicator  of  the  condition  of  a  thawed  beet.  A 
fresh  undamaged  sugarbeet  will  have  a  firm  texture  unless  it 
has  become  wilted  through  loss  of  moisture.  In  contrast,  a 
beet  which  has  become  damaged  by  freezing  and  thawing  will 
have  lost  its  rigidity  and  will  have  a  soft  texture. 

The  texture  of  a  thawed  beet  is  important  in  the  sugar 
manufacturing  process,  where  the  beet  must  be  cut  into 
cossettes,  before  the  sugar  is  extracted.  If  the  beet  is  too 
soft,  it  will  be  difficult  to  cut  cossettes  into  the 
required  shape. 

As  discussed  in  more  detail  in  Appendix  B,  the  severity 
of  freeze-damage  depends  on  the  following  factors: 

1.  the  rate  at  which  the  sample  is  frozen  and  thawed, 

2.  the  lowest  temperature  reached  by  the  sample  during 
freezing,  and 

3.  the  temperature  at  which  the  sample  is  stored. 

It  was  decided  to  investigate  the  effect,  which  a  number  of 
different  freezing  and  thawing  conditions,  would  have  on  the 
texture  of  sugarbeet  samples.  It  was  hoped  to  determine 
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freeze-thaw  conditions,  which  would  produce  a  thawed  beet 
with  a  firm  texture  and  little  indication  of  freeze-damage. 
The  freezing  conditions,  which  were  examined,  ranged  from 
rapid  freezing  in  liquid  nitrogen  to  very  slow  freezing 
(0.001  K/min).  Thawing  conditions  varied  from  rapid  thawing 
in  a  microwave  oven  to  very  slow  thawing  (0.001  K/min). 


EXPERIMENTAL  PROCEDURES' 


5 . 1  Preparation  of  Samples  Prior  to  Freezing 

The  sugarbeets,  which  were  used,  were  grown  for  6 
months  in  a  growth  chamber  and  then  stored  at  275  K  until 
required.  All  the  samples  were  cylinders,  which  were  removed 
from  along  the  length  of  a  beet  by  means  of  a  cork  borer. 

The  cylinders  were  then  cut  to  a  length  of  5  cm.  The 
cylinders  were  1  cm  in  diameter  unless  they  were  frozen  in 
liquid  nitrogen.  With  one  exception,  the  liquid  nitrogen 
frozen  samples  had  a  diameter  of  0.66  cm.  This  was  because 
larger  diameter  cylinders  tended  to  crack  when  placed  in 
liquid  nitrogen.  Cylinders,  which  were  to  be  subjected  to 
microwave  thawing  after  freezing  in  liquid  nitrogen,  were 
0.33  cm  in  diameter.  For  diameters  greater  than  0.33  cm,  the 
middle  of  a  microwave  thawed  sample  was  burned  before  the 
ends  were  thawed. 

Excluding  the  microwave-thawed  samples,  all  beet 
cylinders  were  covered  with  two  layers  of  aluminum  foil 
prior  to  freezing.  This  was  to  minimize  oxidation  and 
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dehydration  of  the  beet.  The  mic rowave-thawed  samples  were 
not  covered  during  freezing  and  thawing.  They  were,  however, 
covered  immediately  after  thawing. 

A  26  gauge  copper-constantan  thermocouple  was  used  to 
measure  the  centre  temperature,  during  freezing  and  thawing. 
The  thermocouple  was  pushed  through  one  end  and  along  the 
axis  of  the  cylinder  until  the  thermocouple  junction  was  at 
the  cylinder  centre.  To  facilitate  the  insertion  of  the 
thermocouple,  a  hole  had  been  previously  made  along  the  axis 
of  the  cylinder  with  an  18  gauge  syringe  needle. 


5 . 2  Freezing  Procedures 

5.2.1  Liquid  Nitrogen  Freezing 

Preliminary  experiments  were  performed  in  which  samples 
were  cut  into  pieces  and  visually  examined  immediately  after 
they  were  removed  from  liquid  nitrogen.  These  experiments 
indicated  that  0.33  cm  diameter  cylinders  were  completely 
frozen  after  15  s  in  liquid  nitrogen.  Consequently,  samples, 
which  were  to  be  subsequently  thawed  in  a  microwave  oven, 
were  frozen  for  15  s  in  liquid  nitrogen.  Samples,  to  be 
thawed  by  ambient  air  or  in  a  water  bath,  were  placed  in  a 
horizontal  position  on  a  sample  holder  (Figure  1)  and  then 
frozen  in  liquid  nitrogen  until  the  centre  temperature 
reached  268  K.  On  removal  from  the  liquid  nitrogen,  the 
centre  temperature  continued  to  fall  to  263  K  or  below. 
Samples  to  be  thawed  at  0.001  K/min  or  at  0.008  K/min  were 
frozen  in  a  vertical  position.  They  could  then  be 


' 


1 


22 


SUGARBEET  CYLINDERS 


TOP  VIEW 


STYROFOAM 


FRONT  VIEW  SIDE  VIEW 


FIGURE  h  SAMPLE  HOLDER  FOR  LIQUID 

NITROGEN  FREEZING 


transferred  quickly  to  the  centrifuge  tubes  (see  Figure  2) 
in  which  they  were  thawed.  They  too  were  removed  from  the 
liquid  nitrogen  when  the  centre  temperature  was  268  K. 
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It  was  observed  that  liquid  nitrogen  freezing  lowered 
the  centre  temperature,  of  0.66  cm  diameter  cylinders,  at  a 
rate  of  approximately  50  K/min  between  269  K  and  267  K. 

5.2.2  Freezing  in  243  K  Alcohol  Bath 

The  cylinders  were  placed  in  a  horizontal  position  on  a 
sample  holder,  similar  to  that  shown  in  Figure  3.  The 
cylinders  were  then  frozen  in  a  NESLAB  LT-50  bath,  which 
contained  95%  alcohol  at  243  K.  They  were  removed  from  the 
bath,  for  thawing,  when  their  centres  reached  263  K.  During 
freezing  in  the  243  K  bath,  the  centre  temperature  fell  at 
the  rate  of  2.5  K/min  between  269  K  and  267  K. 

5.2.3  Freezing  at  0.008  K/min  or  0.001  K/min 

The  samples  were  placed,  in  a  vertical  position  (Figure 
2),  in  closed  50  ml  centrifuge  tubes  (to  prevent  attack  by 
alcohol)  which  were  kept  in  the  NESLAB  LT-50  bath  at  273  K. 
When  the  samples  attained  a  temperature  of  273  K  the  bath 
temperature  was  programmed  to  fall  at  the  desired  rate 
(0.008  K/min  or  0.001  K/min).  When  the  bath  temperature 
reached  263  K,  the  samples  were  removed  and  thawed. 
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FIGURE  3  :  SAMPLE  HOLDER  FOR  FREEZING 


IN  ALCOHOL  BATH 


FRONT  VIEW  SIDE  VIEW 
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5 . 3  Thaw i nq  Procedures 

5,3.1  Thawing  i n  Microwave  Oven 

Microwave  thawing  was  performed  in  a  TOSHIBA  DELUXE 
SERIES  microwave  oven.  This  oven  has  nine  levels  of  power 
available  for  heating.  The  lowest  power  level  was  used  in 
these  experiments. 

The  0.33  cm  diameter  samples  were  placed  in  liquid 
nitrogen  for  15  s  and  then  transferred  immediately  to  the 
microwave  oven  for  75  s.  Samples  placed  in  the  microwave 
oven  for  less  than  75  s  were  still  partially  frozen  on  being 
removed.  In  contrast,  samples  left  in  the  oven  for  periods 
longer  than  75  s  were  too  hot  to  touch. 

5.3.2  Thawing  in  303  K  0£  310  K  Water  Bath 

The  samples  were  thawed  in  water,  maintained  at  either 
303  K  or  310  K,  in  a  BRAUN  THERMOMIX  1420  constant 
temperature  bath. 

Between  268  K  and  270  K,  the  temperature  at  the  centre 
of  the  sample  rose  at  a  rate  of  2.5  K/min. 

5.3.3  Slow  Thawing 

The  samples  were  placed  on  a  bench  in  the  laboratory 
until  they  were  thawed. 

The  centre  temperature  was  observed  to  rise,  at  a  rate 
of  0.5  K/min  between  268  K  and  270  K. 
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5.3.4  Thawing  at  0.008  K/min  or  0.001  K/min 

Closed  centrifuge  tubes  (capacity  50  ml)  were  submerged 
in  263  K  alcohol,  in  the  LT-50  bath.  The  tubes  contained 
sufficient  weights  to  allow  them  to  float  vertically  in  the 
bath.  Samples,  prefrozen  to  263  K,  in  a  prescribed  manner, 
were  placed  in  the  tubes.  The  bath  temperature  was  then 
programmed  to  rise  at  the  desired  rate  (0.008  K/min  or  0.001 
K/min).  When  the  bath  temperature  reached  273  K,  the  samples 
were  removed  and  brought  to  room  temperature  in  a  310  K 
water  bath. 
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5 . 4  Texture  Measurements 

After  thawing,  the  cylinders  were  cut  into  five  pieces, 
each  1  cm  long.  The  three  interior  pieces  were  then  tested 
on  an  INSTRON  UNIVERSAL  TESTING  MACHINE.  This  machine  has 
been  widely  used  (25,26,80,108,167)  in  Food  Science  to 
examine  the  effect  of  various  treatments  on  food  texture. 

The  Instron  was  arranged  as  in  Figure  4.  A  sugarbeet 
cylinder  (1  cm  diameter,  1  cm  long)  was  placed  on  the 
compression  table.  The  moving  crosshead  was  then  brought 
down  slowly,  using  manual  controls,  until  the  anvil  just 
touched  the  sample.  The  crosshead  was  then  lowered  at  a 
speed  of  0.5  cm/min  until  the  cylinder  was  compressed  to 
0.24  cm. 

When  a  sample  is  compressed  it  exerts  a  force,  which  is 
detected  by  the  load  cell.  The  load  cell  essentially 
consists  of  four  strain  gauges  attached  to  a  cantilever 
beam.  The  strain  gauges  are  arranged  as  a  Wheatstone  Bridge. 
Any  deflection  in  the  cantilever  beam,  due  to  a  force 
exerted  on  the  load  cell  by  a  sample  being  compressed,  will 
change  the  resistance  of  each  strain  gauge  in  proportion  to 
the  exerted  force.  The  change  in  the  strain  gauges  causes 
the  bridge  to  become  unbalanced  and  this  results  in  a  change 
in  the  output  voltage  to  the  recorder.  The  Instron  is 
calibrated  so  the  recorder  output  indicates  force  versus 
time.  As  the  crosshead  speed  is  known,  this  output  can  be 
easily  converted  to  a  stress-strain  relationship. 


' 

1 


COMPRESSIVE  STRESS 
(  N/cm2 ) 


29 


MOVING  CROSSHEAD 

ANVIL 

SAMPLE 

COMPRESSION  TABLE 
LOAD  CELL 


FIGURE  4*  SAMPLE  TESTING  ON  THE  INSTRON 


STRAIN 


FIGURE  5:  TEXTURE  PROFILES 
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5 . 5  Experimental  Results 

Prior  to  performing  the  freeze-thaw  experiments, 
preliminary  stress  measurements  were  taken  on  unfrozen 
samples.  This  was  to  test  the  variability  of  samples  within 
a  beet  and  between  different  beets.  The  results  are  listed 
in  Table  1.  This  table  shows  the  compressive  stress  of  1  cm 
long  cylinders,  after  they  have  been  compressed  to  a  length 
of  0.9  cm.  From  this  table,  it  can  be  seen  that  an  unfrozen 
beet  cylinder  has  a  compressive  stress  of  about  65  N/cm2  at 
a  strain  of  -10%.  However,  compressive  stress  measurements, 
for  different  beets,  can  vary  between  54.4  and  72.8  N/cm2. 
This  is  a  variation  of  about  35%.  Within  a  particular  beet, 
variations  of  the  order  of  20%  were  obtained. 

Typical  results  for  fresh,  and  freeze-damaged  beets  are 
shown  in  Figure  5.  Point  A  is  called  the  "brittleness 
point".  At  this  point  the  tissue  of  the  beet  collapses  and 
the  compressive  stress  drops  sharply.  The  brittleness  point 
of  a  freeze-damaged  beet  occurs  at  a  strain  of  about  -70%, 
whereas  for  a  fresh  beet  it  is  observed  at  a  strain  of  about 
-25%.  Also  the  compressive  stress  at  the  brittleness  point 
is  smaller  for  a  freeze-damaged  beet.  For  low  values  of 
strain,  the  curve  for  an  unfrozen  beet  is  much  steeper  than 
that  of  a  freeze-damaged  beet.  Consequently,  it  was  decided 
to  use  the  value  of  the  compressive  stress,  at  a  strain  of 
-10%,  as  an  indicator  of  freezing  injury.  Low  values  of 
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TABLE  1:  COMPRESSIVE  STRESS  OF  UNFROZEN  SAMPLES 


Beet  No.  of  Samples1  Compressive  Stress  (N/cm2) 


at  a  Strain 

of  - 

1 

14 

68.8  ± 

11.7 

2 

14 

72.8  ± 

11.5 

3 

14 

65.4  ± 

5.7 

4 

10 

67.1  ± 

3.1 

5 

15 

54.4  ± 

4.8 

i 


No.  of  Samples 


number  of  1  cm  long,  1  cm  diameter 
cylinders  which  were  examined. 
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this  stress  can  be  considered  synonymous  with  softness. 

Compressive  stress  measurements  (at  a  -10%  strain),  for 
unfrozen  and  thawed  beet  cylinders,  are  tabulated  in  Table 
2.  The  unfrozen  and  thawed  cylinders,  in  each  row  of  this 
table,  were  taken  from  the  same  beet.  A  wide  range  of 
freeze-thaw  conditions  was  considered.  They  all  produced 
soft  cylinders,  with  a  compressive  stress  which  was  an  order 
of  magnitude  less  than  that  of  an  unfrozen  cylinder.  Two 
trends  are  noticeable  in  Table  2.  The  first  one  is  for 
cylinders  frozen  in  liquid  nitrogen.  For  these  beets,  it 
appears  that,  with  the  exception  of  microwave  thawing,  the 
faster  the  thawing  rate  the  firmer  will  be  the  thawed  beet. 
With  microwave  thawing,  0.33  cm  diameter  cylinders  were  used 
and  this  may  have  affected  the  stress  measurements.  The 
second  trend  was  for  very  slowly  frozen  samples  (0.008  K/min 
or  0.001  K/min).  Firmness  appeared  to  increase  as  the 
thawing  rate  was  decreased.  Thus,  these  texture  experiments 
supported  the  general  rule,  explained  in  Appendix  B,  that 
fast  freezing  should  be  followed  by  fast  thawing  and  slow 
freezing  by  slow  thawing.  However,  the  experiments  were 
repeated  an  insufficient  number  of  times  to  conclusively 
prove  this  rule  for  sugarbeets.  It  was  not  considered 
worthwhile  to  carry  out  more  repeats  to  determine  which 
freeze-thaw  conditions  produced  the  "least  soft"  beets. 
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TABLE  2:  COMPRESSIVE  STRESS  OF  FROZEN  -  THAWED  SAMPLES 


FREEZE  -  THAW  CONDITIONS 


Liquid  Nitrogen  Freezing 
with  Thawing: 

-in  microwave  oven2 
-in  310  K  water  bath3 
-at  0.5  K/min3 
-at  0.008  K/min3 
-at  0.001  K/min3 
Freezing  in  243  K  Bath 
and  thawing: 

-in  303  K  water  bath 
-at  0.5  K/min 
Freezing  to  263  K  at 
0.008  K/min  and  thawing: 

-in  310  K  water  bath 
-at  0.5  K/min 
-at  0.008  K/min 


COMPRESSIVE  STRESS  (N/cm2 ) 


AT  A  STRAIN  OF  -10% 


UNFROZEN  THAWED 


70.213.5 

(2)  1 

2.311.2 

(16) 

- 

6.713.5 

(12) 

74.513.8 

(4) 

5.312.6 

(8) 

74.513.8 

(4) 

3.510.6 

(9) 

84.510 

(2) 

3.810.9 

(7) 

67, 

.114.7 

(9) 

4.110.4 

(9) 

76. 

.018.4 

(8) 

7.411.0 

(9) 

80.7+2.8 

(4) 

5.110.1 

(7) 

80.7+2.8 

(4) 

8.512.7 

(8) 

80.712.8 

(4) 

10.614.2 

(9) 
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Freezing  to  263  K  at 

0.001  K/min  and  thawing  UNFROZEN  THAWED 


-in  310  K  water  bath 
-at  0.001  K/min 


82 . 1±3 . 9  (2)  3 . 7±0 . 8  (9) 

82 . 1±3 . 9  (2)  5 . 6±1 . 3  (9) 


1  The  figure  inside  the  brackets  is  the  number  of 
1  cm  long  cylinders  tested. 


2  indicates  0.33  cm  diameter  cylinders  were  used. 


3  indicates  0.66  cm  diameter  cylinders  were  used. 
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The  above  experiments  indicated  that  thawed  beets  are 
soft,  irrespective  of  the  freezing  and  thawing  method  used. 
This  softening  of  beet  tissue  is  directly  related  to  water 
loss  from  the  beet  cells.  During  freezing,  water  leaves  the 
interior  of  the  cell  and  freezes  in  the  intercellular  region 
or  on  the  surface  of  the  cylinder.  On  thawing,  all  of  this 
water  may  not  re-enter  the  cell.  This  may  be  due  to  cellular 
damage  (discussed  in  Appendix  B),  or  to  nonequilibrium 
conditions,  where  the  water  has  insufficient  time  to 
re-enter  the  cell  during  thawing.  In  either  case,  the  water 
is  lost  as  "drip" .  Loss  of  water  reduces  the  firmness  of  the 
tissue.  Some  experiments  were  performed  to  test  whether  the 
softness  of  thawed  beets  was  due  to  nonequilibrium  thawing 
conditions  or  to  cellular  damage.  Thawed  cylinders  were  left 
uncovered  in  water,  for  periods  of  up  to  16  hours,  to 
determine  if  they  would  absorb  water  and  increase  their 
firmness.  It  was  observed  that  soaking  in  water  did  not 
appreciably  change  the  compressive  stress  of  the  samples. 
Consequently,  it  was  concluded  that  loss  of  firmness  in 
thawed  cylinders  is  caused  by  cellular  damage.  It  is, 
therefore,  expected  that  all  the  freeze-thaw  conditions, 
examined  in  this  thesis,  result  in  beets  which  deteriorate 
rapidly  and  become  unprocessable.  Therefore,  frozen  beets 
should  not  be  allowed  to  thaw  before  processing. 
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6.  STORAGE  EXPERIMENTS 

It  is  generally  accepted  (15,58,77,201,213,230)  that  if 
beets  are  stored  frozen  then  the  rate  of  sucrose  degradation 
is  decreased.  As  part  of  this  thesis,  experiments  were 
performed  to  determine  the  effect  of  a  number  of  variables 
on  sucrose  loss  in  frozen  storage. 

The  amount  of  sucrose  lost  during  frozen  storage  should 
depend  primarily  on  the  storage  temperature  and  the  storage 
time.  To  consider  the  effect  of  temperature  some  samples 
were  stored  at  263  K  and  others  at  248  K.  The  effect  of 
storage  time  was  removed  from  the  study  by  storing  all  the 
samples  for  six  months.  This  time  was  selected  because  the 
use  of  frozen  storage  will  probably  increase  the  factory 
operating  period  from  four  to  six  months.  Frozen  storage 
will  be  limited  to  six  months  because  outdoor  air 
temperatures  in  April  may  be  high  enough  to  cause  thawing  of 
the  beets.  If  shorter  storage  times  were  used  in  the 
experiments  then  smaller  changes  in  sugar  content  would 
result.  These  small  changes  would  be  more  difficult  to 
measure  experimentally. 

In  addition  to  storage  time  and  temperature,  the  rate 
at  which  beets  are  frozen  may  affect  the  amount  of  sugar 
lost  during  storage.  As  explained  in  Appendix  B,  the  rate  at 
which  a  plant  is  frozen  can  determine  the  amount  of  injury 
it  experiences  during  freezing.  It  would  appear  likely  that 
beets  experiencing  the  most  damage  during  freezing  would  be 
the  least  stable  during  storage  and  would  show  the  greatest 
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decline  in  sugar  content. 

Sucrose  loss,  during  frozen  storage,  may  also  depend  on 
the  lowest  temperature  attained  by  the  beets  during 
freezing.  For  example,  consider  the  effect  of  low 
temperature  on  invertase.  Fennema  (77)  stated  that  the 
activity  of  invertase  can  be  considerably  reduced  by 
exposing  it  to  a  temperature  of  233  K.  This  decrease  in 
activity  may  be  important  as  invertase  catalyzes  the 
hydrolysis  of  sucrose  to  glucose  and  fructose.  Consequently, 
if  beets  are  frozen  to  233  K  prior  to  storage,  sucrose 
losses  may  be  decreased.  However,  exposing  the  beets  to  this 
lower  temperature  may  increase  the  damage  inflicted  on  the 
beet  cells. 


6 . 1  First  Set  of  Storage  Experiments 

In  these  experiments,  beets  were  subjected  to  a  number 
of  treatments.  A  "treatment"  is  considered  to  be  defined  by 
the  freezing  and  storage  conditions  which  a  sample 
experiences.  The  eight  treatments  considered  are  now 
defined : 

Treatment  1 : Freeze  and  store  samples  in  liquid  nitrogen  at 
77  K. 

Treatment  2 :*  1  Freeze  samples  to  233  K  and  store  at  263  K  (± 

1.5  K)  . 


1Samples  subjected  to  treatments  2-5  were  frozen  in  a  233  K 
alcohol  bath. 
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Treatment  3  rFreeze  samples  to  233  K  and  store  at  248  K  (± 

IK)  . 

Treatment  4 ; Freeze  samples  until  the  centre  temperature 
reaches  263  K  and  store  at  263  K  (±  1.5  K). 

Treatment  5  rFreeze  samples  until  the  centre  temperature 
reaches  248  K  and  store  at  248  K  (±  1  K). 

Treatment  6  rFreeze  samples  at  a  rate  of  0.004  R/min  to  262  K 
and  store  at  263  K  (±  1.5  K). 

Treatment  7  rFreeze  samples  at  a  rate  of  0.004  K/min  to  247  K 
and  store  at  248  K  (±  1  K). 

Treatment  8_r_Freeze  samples  in  liquid  nitrogen  for  30  s  and 
store  at  248  K  (±  1  K). 

Any  sucrose  loss  resulting  from  treatments  2-8  was  to 
be  evaluated  by  comparing  them  to  treatment  1.  It  was 
assumed  that  beet  cylinders,  subjected  to  treatment  1,  would 
not  experience  a  decrease  in  sucrose  content  since  they  were 
stored  in  liquid  nitrogen.  Treatments  2  and  3  were 
undertaken  to  test  if  exposing  the  samples  to  233  K  would 
reduce  sucrose  losses  during  storage.  Treatments  4  and  5,  6 
and  7,  and  8  were  investigated  to  determine  the  effect  of 
fast  (10  K/min),  very  slow  (0.004  K/min),  and  very  rapid  (50 
K/min)  freezing  on  sucrose  losses  in  storage. 

6.1.1  Experimental  Procedures 

A  Randomized  Block  Experiment  was  set  up  for  the  eight 
treatments.  A  replicate  consisted  of  a  set  of  three  beet 
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cylinders,  each  5  cm  long  and  1  cm  in  diameter.  Every  set 
contained  one  cylinder  from  each  of  three  beets.  The  eight 
sets  were  then  assigned  at  random  to  the  eight  different 
treatments.  This  procedure  was  followed  twelve  times  to 
produce  twelve  replicates. 

The  beets,  used  in  replicates  4-6,  were  taken  from  a 
storage  pile,  outside  the  sugar  factory  in  Taber.  They  were 
stored  at  293  K  for  three  weeks,  prior  to  freezing.  Unless 
stated  otherwise,  all  the  other  replicates  were  taken  from 
beets  grown  for  six  months  in  a  growth  chamber  and  stored  at 
275  K  until  they  were  removed  for  freezing. 

Deta i led  Description  of  Treatments 
In  all  the  following  treatments,  the  beet  cylinders  were 
covered  by  two  layers  of  aluminum  foil  in  order  to  minimize 
dehydration  and  surface  oxidation  in  storage.  It  also 
afforded  some  protection  to  the  cylinders  from  the  alcohol 
in  which  they  were  frozen. 

For  the  four  intermediate  rate  freezing  treatments 
(2-5),  the  samples  were  put  in  the  sample  holder  (Figure  3) 
which  was  then  placed  in  a  233  K  alcohol  bath  for  various 
periods  of  time,  depending  on  the  treatment.  The  samples 
were  quickly  transferred,  after  being  shaken  to  remove 
alcohol  from  their  surface,  to  one  of  the  storage  freezers. 
Preliminary  experiments  indicated  the  centre  temperature 
would  be  at  263  K,  248  K,  and  233.5  K  after  2  min  ±  10  s,  2 
min  30  s  ±  15  s,  and  4  min  40  sec  ±  30  s,  respectively,  in 
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the  233  K  bath.  Based  on  these  preliminary  experiments,  the 
intermediate  rate  freezing  treatments  in  Table  3  were 
performed . 

The  ultraslow  freezing  treatments  (6  and  7)  were 
performed  by  placing  the  cylinders  in  closed  50  ml 
centrifuge  tubes  (Figure  2).  The  tubes  were  then  put  in  a 
273  K  alcohol  bath.  The  bath  temperature  was  programmed  to 
fall  at  a  rate  of  0.004  K/min.  When  the  bath  temperature 
reached  262  K  the  treatment  6  samples  were  transferred  to 
the  263  K  storage  freezer.  At  a  bath  temperature  of  247  K 
the  treatment  7  samples  were  moved  to  the  248  K  storage 

TABLE  3:  INTERMEDIATE  RATE  FREEZING  TREATMENTS 


TIME  IN 

STORAGE 

233  K 

TEMPERATURE 

BATH(min) 

(K) 

Treatment 

2 

10 

263 

Treatment 

3 

10 

248 

Treatment 

4 

2 

263 

Treatment 

5 

2.67 

248 

freezer.  The  samples  were  frozen  to  1  K  below  their  storage 
temperature.  This  was  because  preliminary  experiments  had 


« 
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shown  that  the  temperature  of  the  samples  rose  about  1  K, 
while  being  transferred  from  the  alcohol  bath  to  the  storage 
f  reezer . 

The  rapidly  frozen  samples  (treatment  8)  were  put  on  a 
holder  (Figure  1)  and  then  into  liquid  nitrogen  for  30  s. 

The  samples  were  next  transferred  to  the  248  K  storage 
freezer.  Preliminary  experiments  had  been  performed  in  which 
samples  were  placed  in  liquid  nitrogen  for  30  s  and  then 
warmed  in  ambient  air.  It  was  observed  that  the  centre 
temperature  was  between  279  K  and  271  K,  on  removal  from  the 
nitrogen.  It,  subsequently,  continued  to  fall,  to  a  value 
between  249  K  and  213  K. 

6.1.2  Experimental  Results 

The  sucrose  contents,  reproducible  to  ±  0.2%  sucrose, 
obtained  after  26  weeks  (±  2  days)  of  storage  are  shown  in 
Table  4.  A  detailed  description  of  the  procedure,  used  to 
estimate  sugar  contents,  is  listed  in  Appendix  C.  From  the 
results  in  Table  4,  an  Analysis  of  Variance  (AOV)  table  was 
constructed  and  it  is  shown  as  Table  5.  From  the  F  values  in 
this  table  it  can  be  seen  that  there  is  a  significant 
difference  between  treatments  at  the  1%  level  i.e.  there  is 
a  99%  probability  that  a  difference  between  treatments 
exists . 

In  Table  6,  the  8  treatments  are  listed  in  order  of  the 
mean  sugar  contents  they  produced. 
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TABLE  4:  %  SUCROSE  AFTER  STORAGE  (First  Set  of  Experiments) 
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TABLE  5:  ANALYSIS  OF  VARIANCE  TABLE  FOR  FIRST  SET  OF  STORAGE  EXPERIMENTS 
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TABLE  6:  MEAN  SUGAR 

CONTENTS  OF  TREATMENTS  1-8 

MEAN 

TREATMENT 

SUGAR  CONTENT 

( %sucrose ) 

Trt  .1 

17.6 

Trt  .7 

17.6 

Trt  .8 

17.5 

Trt .  6 

17.5 

Trt .  5 

16.9 

Trt .  3 

16.8 

Trt  .4 

16.6 

Trt .  2 

16.4 
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Duncan's  Multiple  Range  Test  (69,70)  was  used  to  determine 
which  treatments  were  significantly  different  from  treatment 
1  (the  "control”).  The  detailed  analysis  is  given  in 
Appendix  D.  The  results  may  be  summarized  as  follows: 

5%  level  of  s i gn i f icance 

Trt.2  Trt.4  Trt.3  Trt.5  Trt.6  Trt.8  Trt.7  Trt.l 


1%  level  of  sign i f icance 

Trt.2  Trt.4  Trt.3  Trt.5  Trt.6  Trt.8  Trt.7  Trt.l 


Any  two  treatments  which  are  underlined  with  the  same  line 
are  not  significantly  different.  If  two  treatments  are  not 
underlined  with  the  same  line  then  they  are  significantly 
different.  With  one  exception,  the  use  of  either  a  1  or  5% 
level  of  significance  lead  to  the  same  conclusions. 

Treatment  2  was  significantly  different  from  treatments  3 
and  5  at  the  5%  but  not  at  the  1%  level. 

It  can  be  reasonably  concluded  that  ultraslow  freezing 
(0.004  K/min)  produced  a  negligible  loss  of  sucrose  for  both 
storage  temperatures  (treatments  6  and  7).  Liquid  nitrogen 
freezing,  followed  by  248  K  storage  (treatment  8)  also  lead 
to  an  insignificant  change  in  sucrose  content.  However,  an 
intermediate  freezing  rate  (treatments  2-5)  caused  a 
significant  sucrose  loss.  This  was  true  for  both  storage 
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temperatures.  Freezing  the  samples  to  233  K  (treatments  2 
and  4),  instead  of  to  the  storage  temperature,  did  not 
appear  to  reduce  the  amount  of  sucrose  lost  during  storage. 
However,  this  sucrose  loss  could  have  been  due  either  to  the 
freezing  rate  or  to  attack  by  alcohol.  In  treatments  2-5  the 
beet  cylinders,  covered  in  aluminum  foil,  were  placed  in  a 
233  K,  95%  ethanol,  bath  for  2-10  minutes  depending  on  the 
treatment.  The  samples  were  then  removed  from  the  bath, 
shaken  to  remove  alcohol,  and  placed  in  50  ml  containers  in 
a  263  K  or  248  K  freezer.  Despite  this  shaking,  a  few  drops 
of  alcohol  did  enter  the  container.  As  these  containers  were 
closed  during  storage,  a  few  drops  of  alcohol  would  have 
been  sufficient  to  saturate  the  air  inside  the  containers. 
Also,  a  small  amount  of  liquid  alcohol  did,  despite  the 
aluminum  foil,  reach  the  surface  of  the  cylinders.  Prior  to 
an  examination  of  the  experimental  results,  it  was  assumed 
that  the  minute  amounts  of  alcohol  involved  would  not  affect 
the  beets,  especially  as  they  were  being  stored  at  low 
temperatures.  However,  as  the  four  treatments  which  showed  a 
loss  in  sucrose  were  the  only  treatments  in  which  the 
cylinders  were  physically  contacted  by  alcohol,  this  may  not 
have  been  a  valid  assumption. 

The  low  sugar  contents  may  also  have  resulted  from  the 
experimental  procedure  which  was  used  to  measure  them.  As 
part  of  this  procedure,  the  cylinders  were  cut  into  minute 
pieces  in  an  Oster  blender.  Before  blending,  the  aluminum 
foil  cover  was  removed  from  the  cylinders.  The  surfaces  of 
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the  intermediate  rate  frozen  samples  were  moist. 
Consequently,  a  small  amount  of  liquid  ("drip"),  containing 
sugar,  may  have  been  lost  when  the  foil  was  removed.  The  ice 
layer  on  the  surfaces  of  ultraslowly  frozen  samples  may  have 
prevented  any  significant  decrease  in  measured  sugar  content 
due  to  drip  loss. 

As  well  as  being  part  of  the  Randomized  Block  Design, 
treatments  2-7  also  formed  a  Complete  Factorial  Design.  This 
latter  design  consisted  of  two  levels  of  storage  temperature 
(263  K  and  248  K)  and  three  levels  of  freezing  method 
(ultraslow  freezing,  intermediate  rate  freezing  to  the 
storage  temperature,  and  intermediate  rate  freezing  to  233 
K) .  For  these  six  treatments  the  experimental  results  are 
rewritten  in  Table  7  in  the  form  of  a  Factorial  Design. 

Table  8  is  a  summary  of  Table  7.  It  simply  shows  the  column 
totals  of  Table  7  in  a  more  succinct  form.  The  AOV  table  for 
the  Factorial  Design  is  shown  as  Table  9.  From  the  F  values 
in  this  table  it  can  be  seen  that  there  is  a  significant 
difference  between  different  freezing  methods  and  between 
different  storage  temperatures. 
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Duncan’s  Multiple  Range  Test  (Appendix  D)  produced  the 
following  results  for  5  and  1%  levels  of  significance: 

5%  Level  of  S i qn i f icance 

Trt.2  Trt.4  Trt.3  Trt.5  Trt.6  Trt.7 


1%  Level  of  Siqni f icance 

Trt.2  Trt.4  Trt.3  Trt.5  Trt.6  Trt.7 


From  the  results  of  Duncan's  Test,  it  can  be  concluded  that 
when  using  an  intermediate  rate  of  freezing  it  is  immaterial 
whether  the  sample  is  frozen  to  233  K  (treatments  2  and  4) 
or  merely  to  the  storage  temperature  (treatments  3  and  5). 
Also,  ultraslow  freezing  leads  to  a  significantly  higher 
sucrose  content  than  intermediate  rate  freezing. 

The  low  F  value  for  the  Freezing  X  Storage  interaction 
indicates  that  the  interaction  between  the  freezing  method 
and  the  storage  temperature  is  not  significant.  Thus,  the 
freezing  conditions  which  produce  the  highest  sugar  content 
after  storage  at  263  K  will  also  be  the  best  conditions  for 
storage  at  248  K.  Similarly,  248  K  storage  following  any 
given  freezing  conditions  will  produce  higher  sugar  content 
beets  than  those  subjected  to  263  K  storage,  following  the 
same  freezing  conditions. 

The  conclusions  from  the  Randomized  Block  and  Factorial 


' 

' 
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Designs  are  now  summarized: 

1.  It  is  significantly  better  to  store  beets  at  248  K  than 
at  263  K. 

2.  There  is  an  insignificant  interaction  between  freezing 
method  and  storage  temperature. 

3.  Ultraslow  freezing  produces  higher  sucrose  content  beets 
than  intermediate  rate  freezing. 

4.  Freezing  to  233  K,  instead  of  to  the  storage 
temperature,  has  an  insignificant  effect  on  the  sucrose 
content . 

However,  as  noted  in  the  Randomized  Block  analysis,  the 
lower  sucrose  contents  resulting  from  intermediate  rate 
freezing  may  have  been  due  to  alcoholic  contamination  or  to 
drip  loss  during  the  sugar  content  determination. 

6.1.3  Physical  Condi t i on  of  Beets  After  Storage 

A  number  of  beet  cylinders  were  visually  examined 
immediately  after  storage.  It  was  apparent  that  the 
condition  of  these  cylinders  depended  on  the  freezing  method 
and  storage  temperature  used.  Liquid  nitrogen  frozen  samples 
stored  at  77  K  or  at  248  K  were  as  white  (immediately  after 
thawing)  as  when  they  were  placed  in  storage. 

Ultraslowly  frozen  cylinders,  stored  at  248  K,  had 
slightly  off-white  coloured  surfaces  but  their  interior  was 
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white.  Ultraslowly  frozen  cyl inders , stored  at  263  K,  had 

some  surface  browning  but  their  interiors  were  white  or  very 

light  brown.  These  samples  had  a  much  thicker  surface  layer 

* 

of  ice  than  cylinders  frozen  more  quickly. 

Intermediate  rate  frozen  samples  stored  at  248  K 
(treatments  3  and  5)  were  generally  as  white  as  the  unfrozen 
samples.  Occasionally  some  browning  occurred.  Intermediate 
rate  freezing,  to  263  K  or  to  233  K,  usually  resulted  in  a 
marked  deterioration  of  at  least  part  of  the  surface 
following  storage  at  263  K.  The  surfaces  of  these  cylinders 
had  either  black,  brown,  or  grey  areas.  This  surface 
colouring  was  also  present  in  the  interior  beneath  it, 
although  generally  in  a  lighter  form.  Usually,  deterioration 
started  at  one  or  both  ends  of  the  cylinder.  This  could  have 
been  due  to  the  faster  freezing  rates  and  lower  temperatures 
which  were  obtained  at  the  ends  during  freezing.  It  could 
also  have  resulted  from  the  fact  that  alcohol  may  have  been 
trapped  under  the  aluminum  foil  at  the  ends  of  the  cylinder. 

6.1.4  Addi t i onal  Experiments 

Some  experimental  difficulties  were  experienced  with 
replicates  1-3  of  treatment  6  and  with  replicates  1-6  of 
treatment  7.  Consequently,  it  was  decided  to  carry  out  an 
extra  three  replicates  of  treatment  6  and  an  extra  six 
replicates  of  treatment  7.  The  results  obtained  are  given  in 
Table  10.  These  results  support  the  earlier  conclusion  that 
ultraslow  freezing,  followed  by  six  months  storage,  at 
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TABLE  10:  SUCROSE  CONTENTS  AFTER  6  MONTHS  STORAGE  -  ADDITIONAL  RESULTS 
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either  263  K  or  248  K,  produces  insignificant  changes  in 
sucrose  content.  The  mean  sucrose  content  of  replicates  4-15 
for  treatments  1  and  6  were  17.2  and  17.1%,  respectively. 

For  replicates  7-18,  both  treatments  1  and  7  had  a  mean 
sucrose  content  of  17.8%. 

The  behaviour  of  refrozen  samples  during  storage  was 
tested  for  three  replicates.  It  was  reasoned  that  refreezing 
the  samples  would  cause  extensive  tissue  damage  and  would  be 
a  good  test  to  determine  if  badly  damaged  beets  would  lose 
sugar  when  stored  at  248  K.  The  samples  were  frozen  slowly 
(0.004  K/min)  to  249  K  (±  2  K)  and  then  immediately  thawed, 
at  a  rate  of  0.5  K/min,  to  271.5  K.  They  were  refrozen, 
without  delay,  at  a  rate  of  0.25  K/min,  to  247  K  and  stored 
at  248  K  for  six  months.  These  freezing  and  storage 
conditions  were  considered  to  be  treatment  9.  The  results 
are  summarized  in  Table  11.  From  this  table,  it  would  appear 
that  even  refrozen  samples,  stored  at  248  K,  will  have  a 
negligible  loss  in  sucrose.  However,  additional  replicates 
would  be  required  to  validate  this  conclusion. 


% 


TABLE  11:  SUCROSE  CONTENTS  OF  REFROZEN  SAMPLES1 


REPLICATE 

CONTROL 
( Trt . 1 ) 

REFROZEN 
(Trt. 9) 

( %sucrose ) 

( %sucrose ) 

4 

15.7 

15.1 

5 

17.6 

17.7 

6 

16.9 

17.2 

MEAN 

16.7 

16.7 

Taber  beets  were  used  in  these  experiments 
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6 . 2  Second  Set  of  Storage  Experiments 

These  experiments  were  designed  to  test  if  the  sucrose 
lost  in  the  intermediately  fast  frozen  samples  of  the  first 
set  of  experiments  was  actually  due  to  the  freezing  rate.  At 
the  same  time  it  was  decided  to  determine  if  liquid  nitrogen 
frozen  samples  would  lose  sucrose  when  stored  at  263  K.  The 
following  four  treatments  were  considered  in  this  set  of 
experiments : 

Treatment  1 :  Freeze  and  store  samples  in  liquid  nitrogen  at 
77  K. 

Treatment  10 :  Freeze  cylinders  in  liquid  nitrogen  and  store 
them  at  264  K  (±  1.5  K). 

Treatment  11 :  Freeze  cylinders  in  a  233  K  alcohol  bath  until 
the  centre  reaches  248  K  and  then  store  at  250  K  (±1  K) . 
Treatment  12 :  Freeze  cylinders  in  a  233  K  alcohol  bath  until 
the  centre  reaches  263  K  and  then  store  them  at  264  K  (±1.5 
K)  . 


Prior  to  freezing,  the  treatment  11  and  12  samples  were 
wrapped  in  aluminum  foil  and  placed  in  sealed  polystyrene 
tubes  to  prevent  contamination  by  alcohol.  The  polystyrene 
tubes  reduced  the  freezing  rate  by  a  factor  of  three.  The 
time  for  the  geometric  centre  to  go  from  269  K  to  267  K 
increased  from  15  s  to  45  s. 


TX 
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For  each  replicate,  four  groups  of  three  cylinders  were 
taken  from  two  beets.  Each  group  contained  two  cylinders 
from  the  first  beet  and  one  from  the  second.  The  four  groups 
were  then  assigned  at  random  to  the  four  treatments. 

6.2.1  Experimental  Procedures 

Treatment  10  consisted  of  freezing  the  samples  (covered 
with  aluminum  foil)  in  liquid  nitrogen  for  30  s  and  then 
transferring  them  rapidly  to  the  264  K  storage  freezer. 

In  treatments  11  and  12,  the  samples  were  wrapped  in 
aluminum  foil  and  placed  in  polystyrene  tubes  (7.5  cm  long, 

1  cm  I.D.)  which  were  then  sealed.  An  epoxy  glue,  which 
dried  in  5  min,  was  used  to  prevent  alcohol  leaking  into  the 
tube  (Figure  6).  The  tubes  were  then  put  in  the  holder 
(Figure  7)  which  kept  them  in  a  horizontal  position  in  the 
233  K  alcohol  bath.  In  treatment  11,  the  tubes  were  in  the 
bath  for  4.5  min  before  being  moved  to  the  264  K  storage 
freezer.  In  treatment  12,  the  tubes  were  in  the  bath  for  5.5 
min  before  being  transferred  to  248  K  storage.  It  had 
earlier  been  observed  that  the  centre  of  the  sample  would 
reach  258  ±  3  K  and  244  ±  2  K  after  4.5  and  5.5  min, 
respectively  in  the  bath. 

The  procedure,  used  to  determine  sugar  content,  was 
modified  for  these  experiments.  The  sugarbeet  cylinders  were 
placed  in  the  blender  jars  without  first  having  their 
aluminum  foil  cover  removed.  This  modification  prevented  any 
drip  loss  of  sugar. 
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FIGURE  6:  ARRANGEMENT  TO  PREVENT 


ALCOHOL  CONTAMINATION 
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BOTTOM  VIEW 
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FRONT  VIEW 
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SIDE  VIEW 


FIGURE  7:  SECOND  SAMPLE  HOLDER  FOR 

FREEZING  IN  ALCOHOL  BATH 
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6.2.2  Experimental  Results 


Physical  Condition  of  Beets  After  Storage 

Cylinders  frozen  in  liquid  nitrogen  and  stored  at  263  K 
(treatment  10)  usually  contained  golden  and  white  regions, 
both  at  the  surface  and  in  the  interior. 

Cylinders  frozen  at  an  intermediate  rate  and  stored  at 
250  K  (treatment  11)  were  all  white. 

The  treatment  12  cylinders  (frozen  at  an  intermediate 
rate  and  stored  at  264  K)  had  a  similar  appearance  to  those 
subjected  to  treatment  10.  However,  dark  brown  regions  were 
sometimes  obtained.  Generally,  the  treatment  12  cylinders 
showed  more  sign  of  deterioration  than  treatment  10 
cylinders . 

Sucrose  Content 


The  sucrose  contents,  which  were  obtained  after  180  ±  4  days 
of  storage,  are  summarized  in  Table  12.  This  data  was  used 
to  construct  Table  13,  which  is  an  Analysis  of  Variance 
table.  From  the  F  values  in  this  table,  it  can  be  concluded 
that  a  significant  difference  between  treatments  existed  at 


the  1%  level. 
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TABLE  12:  SUCROSE  CONTENTS1  AFTER  180  DAYS  OF  STORAGE 


-  RANDOMIZED  BLOCK  DESIGN 


Replicate 

Trt .  1 

Trt  .10 

Trt .11 

Trt  .12 

1 

17.4 

17.6 

17.7 

17.2 

2 

18.5 

18.3 

18.2 

18.0 

3 

17.2 

17.3 

17.6 

17.2 

4 

17.4 

17.8 

17.9 

17.8 

5 

16.9 

16.6 

16.7 

16.0 

6 

17.6 

17.7 

17.9 

17.5 

7 

17.1 

16.9 

16.8 

16.2 

8 

17.5 

17.4 

17.7 

17.0 

i 


Results  expressed  as  %sucrose  by  weight. 
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TABLE  13:  ANALYSIS  OF  VARIANCE  TABLE  FOR  SECOND  SET  OF  STORAGE  EXPERIMENTS 
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The  four  treatments  are  now  listed  in  order  of  their 
mean  sugar  contents: 

Trt . 12  Trt.10  Trt.l  Trt.ll 
17.1  17.5  17.5  17.6 

The  results  of  a  Duncan's  Multiple  Range  Test,  performed  in 
Appendix  D,  are  summarized  below: 

5%  Level  of  S iqn i f icance 

Trt. 12  Trt.10  Trt.l  Trt.ll 


1%  Level  of  S i qn i f icance 

Trt. 12  Trt.10  Trt.l  Trt.ll 


From  the  above  results,  it  is  concluded  that  treatments  10 
and  11  caused  a  negligible  loss  in  sugar.  Treatment  12  is 
significantly  different  from  treatment  1  at  the  5%  level  but 
not  at  the  1%  level  of  significance.  However,  it  must  be 
mentioned  that  replicates  5-8  of  treatment  12  were  stored  at 
264  K,  whereas  replicates  1-4  were  stored  at  263  K.  This  was 
because  of  the  temperature  variation  within  the  storage 
freezer.  For  replicates  1-4,  treatments  1  and  12  both  had  an 
average  sucrose  content  of  17.6%.  For  replicates  5-8, 
treatment  1  had  an  average  sucrose  content  of  17.3%,  whereas 


, 


65 


it  was  only  16.7%  for  treatment  12.  It  would,  therefore, 
appear  that  a  negligible  loss  in  sucrose  is  obtained,  for 
treatment  10,  if  the  storage  temperature  is  maintained  at 
263  K.  However,  a  small  rise  in  storage  temperature  may 
cause  a  significant  loss  in  sucrose.  This  theory  is 
supported  by  the  experiments  of  Wyse  (230).  He  observed  that 
the  respiration  rate  of  sugarbeets  decreased  gradually  with 
decreasing  temperature  until  a  temperature  of  265  K  was 
reached.  At  this  temperature  the  respiration  rate  ’’declined 
sharply  to  near  zero". 

From  the  above  experiments,  it  is  concluded  that  an 
insignificant  loss  of  sucrose  results  from: 

1.  six  months  storage  at  250  K  or  at  263  K,  following 
freezing  in  a  233  K  alcohol  bath,  or 

2.  six  months  storage  at  263  K,  following  freezing  in 
liquid  nitrogen. 

6.2.3  Supplementary  Experiments 

The  above  conclusions  are  supported  by  the  results  of 
supplementary  experiments.  The  results  of  additional 
replicates  of  treatments  10  and  11  are  listed  in  Tables  14 
and  15.  The  average  of  the  sugar  contents  of  fourteen 
replicates  (eight  from  Table  12  and  six  from  Table  14)  is 
17.8%  for  both  treatments  1  and  11.  Similarily,  eleven 
replicates  (eight  from  Table  12  and  three  from  Table  15)  of 
treatment  1  and  of  treatment  10  have  the  same  average 
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TABLE  14:  SUCROSE  CONTENTS1 


Replicate  Trt.l 

1  18.5 

2  18.5 

3  18.8 


1  Results  expressed  as  %sucr 


TABLE  15:  SUCROSE  CONTENTS1 


Replicate  Trt.l 

1  17.4 

2  18.7 

3  18.4 

4  18.2 

5  18.5 

6  18.9 


EXTRA  TREATMENT  10  SAMPLES 


Trt.10 

19.0 

18.3 

18.6 

se  by  weight. 


EXTRA  TREATMENT  11  SAMPLES 


Trt .11 

16.8 

18.7 
17.9 
18.1 
18.5 

18.8 


1 


Results  expressed  as  ^sucrose  by  weight. 


67 


sucrose  content  (17.8%).  One  extra  replicate  of  treatment  12 
(stored  at  263  K)  produced  a  sugar  content  of  18.3%,  whereas 
the  treatment  1  sample  had  a  sugar  content  of  18.5%. 


6 . 3  Conclusions  from  Storage  Experiments 

Based  on  all  the  storage  experiments,  it  may  be 
concluded  that  sucrose  losses  will  be  negligible  in  beets, 
stored  for  six  months  at  248  K  or  263  K,  irrespective  of  the 
freezing  method  used.  However,  ultraslow  freezing  has  the 
advantage  of  providing  a  thick  surface  layer  of  ice.  This  is 
important  as  it  may  (at  least  until  the  ice  melts)  reduce 
the  amount  of  sugar  which  is  lost  in  the  beet  factory  to  the 
water,  which  is  used  for  the  transportation  and  washing  of 
the  frozen  beets. 

A  storage  temperature  of  248  K  is  preferred  to  one  of 
263  K  as  it  minimizes: 

1.  the  probability  of  the  beets  becoming  thawed  prior  to 
processing , 

2.  the  quantity  of  sugar  which  is  lost  to  the 
transportation  and  cleaning  water,  and 

3.  the  beet  deterioration  which  occurs  during  storage. 

Of  course  if  mechanical  refrigeration  must  be  used,  storage 
at  248  K  will  be  the  more  expensive. 
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7.  SUGGESTED  EXPERIMENTS  ON  STORAGE  PILES 
The  laboratory  experiments  indicated  that  frozen  storage  can 
reduce  the  sugar  loss  in  storage  piles  from  a  value  of  0.25 
kg/tonne/day  (obtained  in  conventional  outdoor  piles)  to  a 
negligible  amount.  Using  a  refined  sugar  price  of  $l/kg, 
frozen  storage  could  increase  operating  revenues  by 
$0. 25/tonne  of  beets/day  stored.  Frozen  storage  can  be 
expected  to  increase  the  gross  operating  revenues  of  the 
Taber  factory  by  about  $7  million.  It,  therefore,  appears 
worthwhile  to  extend  the  experiments  to  the  frozen  storage 
of  sugarbeet  piles  at  Taber.  A  pilot  scale  study  on  one 
storage  pile  should  first  be  performed.  If  the  results  from 
this  study  are  satisfactory  then  frozen  storage  of  beets 
could  be  carried  out  on  a  large  scale. 

There  are  two  alternative  methods  of  storing  frozen 
piles  of  beets.  The  first  method  entails  the  use  of  cold 
night  air  to  freeze  the  beets,  with  subsequent  outdoor 
storage.  The  main  problem  with  this  method  is  that  it  may  be 
difficult  to  prevent  the  beets  from  thawing  during  periods 
of  warm  weather.  This  is  a  particularly  valid  concern  at 
Taber  because  of  the  presence  of  warm  Chinook  winds.  The 
second  frozen  storage  method  involves  freezing  and  storing 
the  beets  inside  a  building,  cooled  by  mechanical 
refrigeration.  It  should  not  be  necessary  to  use  a 
Quickfreeze  unit  to  freeze  the  beets.  This  is  because 
sucrose  losses  are  negligible  in  frozen  storage  even  if  the 
beets  are  frozen  slowly. 
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7 . 1  Outdoor  Frozen  Storage  of  Beet  Piles 

Research  efforts  should,  at  first,  be  concentrated  on 
the  outdoor  storage  of  beets.  This  is  because  the  cost  of 
freezing  and  storage  is  minimal.  Pile  temperatures  should  be 
monitored  to  determine  if  the  beets  will  thaw  during 
storage.  The  use  of  a  pile  cover  (ice,  snow,  straw,  plastic 
etc.)  in  the  prevention  of  thawing  should  be  investigated. 


7 . 2  Use  of  Mechanical  Refrigeration 

If  it  appears  that  thawing  of  beets  cannot  be  avoided 
in  outdoor  piles  then  the  use  of  mechanical  refrigeration 
should  be  considered.  This  would,  of  course,  significantly 
increase  the  cost  of  storage.  However,  frozen  storage  may 
still  be  economically  attractive. 

7.2.1  Pilot  Study  Using  Mechanical  Refrigeration 

A  pilot  study  could  be  performed  by  storing  2,500 
tonnes  of  beets  at  253  K  in  an  insulated  (10  cm  thick 
styrofoam)  archrib  structure  of  dimensions: 

20  m  wide  X  40  m  long  X  8m  high. 

The  building  would  have  a  six  inch  concrete  floor.  The 
archrib  structure  was  chosen  because  it  is  sturdier  than  an 
inflatable  building.  This  is  especially  relevant  because  of 
the  high  wTind  velocities  obtained  in  the  Taber  area. 
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7.2.2  Requirements  of  Refrigeration  System 

The  refrigeration  system  will  have  its  greatest  load 
requirement  in  October  when  outdoor  conditions  are  warmest. 
The  following  assumptions  were  made  in  estimating  the 
requirements  of  the  refrigeration  system: 

1.  The  average  outdoor  temperature  in  October  is  288  K, 

2.  the  average  outdoor  temperature,  during  the  whole 
storage  period,  is  273  K, 

3.  the  overall  heat  transfer  coefficient  is 

1  J/(s  .  m2  .  K)  for  both  heat  transfer  through  the 
floor  and  through  the  walls, 

4.  the  effect  of  the  sun’s  radiation  is  accounted  for  by 
adding  3  K  to  the  AT  across  the  walls. 

The  capacity  and  average  requirements  of  the  refrigeration 
system  were  calculated  to  be  76.7  kW  (22  tons  of 
refrigeration)  and  46.7  (13  tons),  respectively. 

7.2.3  Refrigeration  Requi rement  to  Cool  Beets  to  253  K 

It  was  assumed  .that  the  2,500  tonnes  of  beets  were  at 
283  K  when  placed  in  storage  in  October.  The  heat,  which 
must  be  removed  to  cool  them  to  253  K,  was  then  calculated 
(using  the  enthalpy  data  in  Table  18)  to  be  7.65  X  105  kj. 
Assuming  the  beets  are  cooled  in  one  week,  the  rate  of  heat 
removal  required  to  freeze  the  beets  to  253  K  is  1.3  kW. 
This  rate  is  negligible  relative  to  the  capacity  of  the 
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refrigeration  system. 

7.2.4  Removal  of  Heat  of  Respiration 

The  heat  of  respiration  of  sugarbeets  at  283  K  is  0.019 
kW/tonne  of  beets  (48).  The  respiration  rate  drops  to  near 
zero  at  about  265  K  (230).  It  is,  therefore,  assumed  that 
the  heat  of  respiration  is  negligible  below  265  K.  The  heat 
of  respiration  to  be  removed  is  greatest  when  the  beets  are 
placed  in  storage  (their  temperature  is  then  283  K) .  The 
maximum  value  of  the  heat  of  respiration,  evolved  by  2,500 
tonnes  of  beets,  was  calculated  to  be  47.2  kW.  This 
refrigeration  load,  although  much  greater  than  that  required 
to  cool  the  beets  to  253  K,  is  less  than  the  capacity  of  the 
refrigeration  system.  The  heat  of  respiration  does  not  alter 
the  required  capacity  of  the  system  as  it  is  negligible  when 
the  refrigeration  unit  is  operating  at  full  capacity.  This, 
of  course,  occurs  when  the  storage  temperature  is  253  K. 

7.2.5  Economics  of  Using  Mechanical  Refrigeration 

This  economic  assessment  applies  to  the  pilot  scale 
equipment  already  described.  It  basically  estimates  the  net 
cost  of  performing  the  experiment.  It,  also,  gives  an 
indication  of  the  economic  feasability  of  storing  frozen 
beets  on  a  large  scale,  using  mechanical  refrigeration. 

The  potential  increase  in  annual  revenues  obtained  by 
storing  the  2,500  tonnes  of  beets  in  the  frozen  state  was 
calculated  after  making  the  following  assumptions: 


V 


1.  the  beets  are  stored  for  120  days,  and 

2.  the  price  of  refined  sugar  is  $l/kg. 
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Annual  operating  costs  were  calculated  after  making  the 
following  assumptions: 

1.  no  additional  labour  will  be  required  to  pile  the  beets 
in  the  archrib  building  over  that  required  for 
conventional  outdoor  piles, 

2.  the  electrical  power  required  to  run  the  refrigeration 
system  is  0.4  kW  per  kW  of  refrigeration,  and 

3.  the  cost  of  electricity  is  1.5  c^/kW-hour. 

The  results  of  the  economic  calculations  are  summarized  in 
Table  16.  From  this  table  it  can  be  quickly  determined  that 
the  increased  revenues  would  pay  for  the  frozen  storage 
system  in  about  3  to  4  years. 

A  capital  investment  of  $200,000  would  be  required  to 
perform  the  pilot  scale  experiment.  However,  if  it  was 
discovered  that  this  form  of  storage  was  not  feasible  then 
the  investment  would  not  be  wasted.  The  archrib  building 
could  always  be  used  as  a  general  storage  building.  The 
refrigeration  system  could  be  sold  for  $50,000. 

As  shown  above,  frozen  storage  of  beets  may  be 
economically  feasible  when  applied  to  an  existing  sugar 
factory.  The  possibilities  for  a  new  sugar  factory  are  even 
better.  Frozen  storage  could  then  be  used  to  extend  the 
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TABLE  16:  SUMMARY  OF  ECONOMIC  CALCULATIONS 


CAPITAL  COSTS 

Archrib  building1 
Refrigeration  system 

OPERATING  COSTS 

Additional  labour 
Electric ity 
Ma  i ntenance 


OPERATING  REVENUES 


Increase  in  refined 
sugar  production 


$100,000 

$100,000 


0 

$1,000 

$10,0002 


$75,000 


1  includes  cost  of  installation  and 
concrete  floor. 


2 


5%  of  capital  cost. 
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operating  period  of  the  factory  from  four  to  six,  or  more, 
months.  This  would  result  in  a  significant  reduction  in  the 
required  daily  capacity  of  the  factory.  The  savings  in 
capital  costs  would  then  be  substantial. 

In  conclusion,  the  frozen  storage  of  sugarbeets  in 
outdoor  piles  is  economically  attractive.  A  pilot  study 
should  be  undertaken  to  determine  if  the  beets  can  be 
maintained  frozen  during  storage  in  outdoor  piles  in  the 
Taber  area.  If  the  beets  are  observed  to  thaw  during  storage 
in  these  piles  then  the  frozen  storage  of  beets  in 
mechanically  cooled  archrib  structures  should  be 
investigated . 


8.  FREEZING  MODELS  -  LITERATURE  REVIEW 
Heat  transfer  models  are  useful  in  estimating  the  freezing 
injury  incurred  by  biological  samples.  The  severity  of  this 
injury  depends  on  the  rate  at  which  the  biological  material 
is  frozen  and  thawed.  Different  regions  of  large  samples  may 
be  injured  to  different  extents  since  they  may  experience 
different  rates  of  temperature  change.  Thus,  in  estimating 
the  overall  freeze-damage  in  an  entire  sample,  it  is 
advantageous  to  be  able  to  predict  how  the  temperature 
varies  with  time  throughout  the  sample  during  freezing  and 
thawing . 

Freezing  of  biological  materials  is  difficult  to  model 
because  all  the  heat  of  fusion  is  not  removed  at  one 
temperature  i.e.  they  do  not  have  a  unique  freezing  point. 
Assuming  there  is  no  supercooling,  the  water  in  a  biological 
material  starts  to  freeze  when  the  temperature  is  lowered  to 
the  initial  freezing  point.  More  water  is  frozen  as  the 
temperature  is  further  reduced.  Approximately  75%  of  a 
food's  water  content  is  frozen  by  268  K. 

Thermal  properties  change  rapidly  during  freezing. 

Water  is  the  largest  component  of  most  foods.  For  example, 
lean  beef  contains  74%  water,  potatoes  78%,  and  apples  84%. 
When  this  water  freezes,  the  thermal  conductivity  of  the 
food  increases  rapidly.  In  fact,  the  thermal  conductivity 
may  triple  in  the  first  5  K  below  the  initial  freezing 
point.  This  dramatic  increase  results  from  the  fact  that  the 
thermal  conductivity  of  ice  is  four  times  that  of  water.  In 
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contrast,  the  specific  heat  of  a  food  falls  during  freezing. 
This  also  stems  from  the  behaviour  of  water,  which 
experiences  a  50%  reduction  in  its  specific  heat  on 
freezing . 

Another  complication  in  food  freezing  is  the  fact  that 
the  heat  of  fusion  is  removed  over  a  range  of  temperatures, 
beginning  at,  or  just  below  (supercooling),  the  initial 
freezing  point.  This  is  sometimes  allowed  for  by  the  use  of 
an  "effective  specific  heat",  which  includes  a  heat  of 
fusion  term.  The  effective  specific  heat  is  a  complicated 
function  of  temperature  just  below  the  initial  freezing 
point . 

A  literature  review  was  undertaken  to  examine  existing 
two  phase  heat  transfer  models.  Some  of  these  models  will 
now  be  discussed. 


8 . 1  Analytical  Solut ions 

Analytical  solutions  are  based  on  a  number  of 
simplifying  assumptions.  For  example,  Plank's  formula  (27), 
which  is  perhaps  the  most  well-known  food  freezing  equation, 
is  based  on  the  following  assumptions: 

1.  the  food  freezes  at  a  unique  temperature,  rather  than 
over  a  range  of  temperatures, 

2.  physical  properties  are  independent  of  temperature, 
although  they  do  have  different  values  in  the  frozen  and 


unfrozen  regions,  and 
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3.  the  food  is  initially  at  its  freezing  point. 


Plank’s  formula  may  be  written  as 


'f 


T,-T 
f  a 
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=  time  for  centre  of  food  sample  to  reach  the  freezing 
point , 

thickness  of  infinite  slab  or  diameter  of  sphere  or 
cylinder , 

=  freezing  point  of  food, 

=  temperature  of  cooling  medium, 
heat  of  fusion  per  unit  mass, 
density , 

heat  transfer  coefficient  at  the  surface, 

1/2  for  infinite  slab  or  infinite  cylinder, 

1/3  for  sphere, 

1/4  for  infinite  slab, 

1/2  for  infinite  cylinder  or  sphere, 

=  thermal  conductivity  of  frozen  food. 


Many  researchers  (27,40,42,49,68,152)  have  modified 
Plank’s  equation.  Cleland  and  Earle  (40,43,44,175)  produced 
perhaps  the  best  adaptation  of  Plank's  formula.  Based  on  a 
least  squares  fit  of  empirical  data,  the  shape  factors  D  and 
G  (in  Plank's  equation)  were  modified  to  allow  for  the 
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situation  where: 

1.  the  food  is  initially  at  a  temperature  above  its 
freezing  point,  and 

2.  the  final  temperature  at  the  centre  of  the  food  is  below 
the  freezing  point. 

The  heat  of  fusion,  L,  was  replaced  by  the  enthalpy  change 
in  going  from  the  freezing  point  to  the  final  temperature. 
Freezing  times  (time  for  the  centre  of  the  food  to  reach 
263  K)  were  predicted  for  Karlsruhe  Test  Substance,  meat  and 
potatoes,  which  were  within  ±  10%  of  experimental  values. 
They  concluded  (42)  that  their  equations  predicted  food 
freezing  times  as  accurately  as  Bonacina  and  Comini's 
numerical  model  (22-24,46).  Their  method  was  also  shown  (44) 
to  be  a  significant  improvement  on  the  equations  of  Plank 
and  Nagaoka  (152). 

Nagaoka  et  al.  (152),  also,  modified  Plank's  equation 
to  allow  for  cooling  above  and  below  the  initial  freezing 
point.  However,  their  equation  was  based  on  experimental 
data  for  fish  freezing  (elliptical  shape).  Consequently, 
Nagaoka' s  equation  is  not  as  generally  applicable  as  that  of 
Cleland  and  Earle. 

Mott  (246)  proposed  an  equation,  without  stating  his 
assumptions,  for  calculating  freezing  times  of  irregularly 
shaped  foods.  However,  this  equation  predicted  freezing 
times  which  were  about  80%  too  long. 
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Starting  with  the  freezing  point  depression  equation, 
Schwartzberg  et  al.  (185)  developed  equations  to  predict  the 
effective  specific  heat  as  a  function  of  temperature  during 
freezing  or  thawing.  They  then  assumed  that  one  phase 
solutions  to  Fourier's  heat  equation  could  be  extended,  for 
low  Biot  numbers,  to  the  two  phase  problem  by  replacing  the 
constant  one  phase  specific  heat  with  the  effective  specific 
heat.  They  also  assumed  that,  during  freezing,  the  thermal 
conductivity  was  the  average  value  for  the  frozen  and 
unfrozen  food.  An  equation  was  then  obtained  which  allowed 
the  calculation  of  the  time  required  for  the  average  food 
temperature  to  reach  a  given  value.  This  equation  only 
applied  to  the  freezing  range.  The  time  to  go  from  the 
starting  temperature  to  the  initial  freezing  point,  or  vice 
versa  during  thawing,  was  calculated  from  an  analytical 
solution  for  single  phase  heat  transfer.  Schwartzberg  tested 
his  equations  against  experimental  thawing  times  for  slabs 
of  chopped  beef.  The  agreement  was  generally  within  about  5% 
but  a  deviation  as  large  as  22%  was  obtained.  He  concluded 
that  his  equations  produced  better  freezing  times  than  the 
methods  of  Plank,  Tao,  and  Mott,  all  of  which  gave  maximum 
deviations  from  the  experimental  values  of  the  order  of  40%. 
The  maximum  deviation  obtained  with  Nagaoka's  equation  was 
27%. 

Bakal  and  Hayakawa  (16,17,96)  developed  semi -empi r ical 
formulae  which  predicted  temperature- t ime  distributions  in 
an  infinite  slab.  Their  method  was  more  realistic  than  that 
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of  Plank  in  that  they  allowed  for  the  fact  that  foods  freeze 
over  a  range  of  temperatures.  However,  their  solutions  were 
based  on  many  other  assumptions.  For  example,  they  assumed 
that : 

1.  Thermal  properties  are  independent  of  temperature, 
although  they  do  have  different  values  in  frozen  and 
unfrozen  regions.  The  value  of  a  thermal  property  in  a 
partially  frozen  region  is  the  average  of  its  values  in 
the  frozen  and  unfrozen  states. 

2.  Before  the  temperature  at  the  centre  of  the  slab  has 
reached  the  initial  freezing  point,  the  fraction  of  an 
elemental  volume  which  is  frozen  is  inversely 
proportional  to  the  distance  between  the  element  and  the 
surface.  After  the  centre  starts  to  freeze,  this 
fraction  becomes  linearly  related  to  the  surface 
temperature . 

Their  model  results  were  generally  within  3  K  of  their 
experimental  values  but  deviations  as  great  as  11  K  were 
observed . 

8 . 2  Numerical  Models 

Semi-empirical  equations  are  useful  in  estimating 
freezing  times.  However,  in  predicting  the  transient 
temperature  distributions  in  a  food  during  a  phase-change, 
one  should  use  a  numerical  approach.  Increased  flexibility 
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is  also  gained  by  using  a  numerical  method.  For  example,  a 
numerical  model  can  allow  for  a  heat  transfer  coefficient 
which  varies  with  time.  This  is  especially  important  for 
free  convection  situations,  such  as  the  still  air  thawing 
carried  out  in  this  thesis. 

Charm  (37,38)  presented  a  numerical  method  which 
predicted  transient  temperature  distributions  in  an  infinite 
cylinder  or  an  infinite  slab.  This  method  applied  to  both 
freezing  and  thawing.  In  developing  his  equations,  he 
assumed  foods  had  a  unique  freezing  point.  It  was  also 
assumed  that  specific  heat  and  thermal  conductivity,  in  both 
the  frozen  and  unfrozen  phases,  were  independent  of 
temperature.  However,  the  value  of  these  thermal  properties 
changed  during  the  transition  from  the  frozen  to  the 
unfrozen  state.  Charm  tested  his  model  against  experimental 
centre  temperature  measurements,  during  freezing  of  cod  and 
haddock  slabs.  For  haddock,  his  model  predicted  a  freezing 
plateau  which  lasted  about  15  min.  In  contrast,  no  plateau 
was  obtained  in  his  experimental  results.  The  long  plateau 
in  the  model  predictions  resulted  from  the  use  of  a  unique 
freezing  point  approach.  Differences  as  great  as  10  K  were 
obtained  between  model  predictions  and  experimental 
measurements . 

Fleming  (83)  developed  a  finite  difference  method  for 
determining  temperature  distributions  during  freezing  of 
bodies  of  arbitrary  geometry.  He  allowed  for  thermal 
properties  which  change  drastically  in  the  freezing  range. 
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The  variation  of  these  properties  with  temperature  was 
approximated  by  a  discontinuous  function  in  order  to  reduce 
the  computing  time  required.  He  tested  his  model  against 
experimental  data  for  foodstuffs  and  concluded  that  his 
model  precision  was  limited  by  "an  adequate  knowledge  of  the 
thermal  properties". 

Cullwick  and  Earle  (54)  used  a  numerical  approach  to 
predict  temperature  as  a  function  of  time  in  an  infinite 
slab.  They  allowed  for  the  variation  of  specific  heat  and 
thermal  conductivity  with  temperature.  They  compared  their 
results  to  experimental  data  for  minced  beef.  For  slabs  7.5 
cm  and  10  cm  thick,  the  experimental  and  model  results  were 
essentially  the  same.  However,  it  should  be  noted  that  the 
heat  transfer  coefficient  at  the  surface,  which  was  used  in 
the  model  equations,  was  not  measured.  Instead  it  was 
selected  so  the  calculated  freezing  time  ("measured  from  the 
beginning  of  the  cooling  and  freezing  cycle  to  the  time  at 
which  the  centre  temperature  falls  from  the  latent  heat 
plateau")  was  equal  to  the  experimental  value.  For  a  12.5  cm 
thick  slab,  model  predictions  deviated  from  experimental 
results  by  as  much  as  10  K. 

Bonacina  and  Comini  (22-24,46)  numerically  solved  the 
one-dimensional  heat  equation: 


[k(T)  ^-T] 



8T  9 


(2) 


p  ( T )  c  ( T ) 


3 1  3y 


by  replacing  3T/3t  by  a  central  difference  value.  The 
derivative  3T/3y  was  averaged  over  three  time  periods,  j-1, 
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j  r  and  j+1,  where  j.t  is  the  time  since  the  start  of 
cooling.  Their  model  is  thus  referred  to  as  a  "three  tier 
time  level  implicit  scheme"  (22).  By  using  this  approach  the 
thermal  conductivity  can  be  evaluated  at  the  present  time 
(time  increment  =  j)  and  this  results  in  a  set  of  linear 
equations.  Linearity  can  also  be  obtained  by  using  an 
explicit  model.  However,  convergence  problems  can  then 
arise.  In  contrast  to  explicit  models,  Bonacina  and  Comini 
(22)  proved  that  their  system  was  unconditionally  stable  and 
convergent.  Clelana  and  Earle  (41)  showed  that  Bonacina  and 
Comini's  model  predicted  freezing  times  which  were  accurate 
to  ±  8%. 

All  the  numerical  models,  which  have  been  discussed, 
used  specific  heat  in  the  fundamental  energy  balance 
equations.  As  stated  by  Fennema  (77)  "it  is  often  more 
accurate  to  deal  with  enthalpy"  rather  than  specific  heat. 
Also,  the  effective  specific  heat  is  difficult  to  use  as  it 
increases  rapidly  when  the  temperature  is  lowered  below  the 
initial  freezing  point.  It  then  falls  quickly  again  as  the 
temperature  is  further  reduced.  A  two-phase  numerical  model 
is  developed  in  the  next  chapter  which  is  based  on  enthalpy 
rather  than  on  spec i f ic.  heat .  Shamsundar  and  Sparrow 
(186,187)  did  derive  an  implicit  model  which  was  based  on 
enthalpy.  However,  their  model  was  for  a  material  with  a 
unique  freezing  point. 
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8.  DERIVATION  OF  MODEL  EQUATIONS 

In  estimating  the  freeze-damage  which  will  occur  in  a  biological 
sample,  it  is  useful  to  know  the  rate  at  which  every  region  of  the  sample 
is  frozen  and  thawed.  A  numerical  heat  transfer  model  which  accomplishes 
this  is  now  developed.  The  model  predicts  the  transient  temperature 

profile  in  an  infinite  cylinder  during  freezing  and  thawing.  The  model 
derivation  is  based  on  the  following  assumptions: 

1.  the  cylinder  is  homogeneous, 

2.  the  density  remains  unchanged,  even  during  a  phase  change,  and 

3.  the  sample  does  not  have  a  unique  freezing  point  but  rather 
freezes  over  a  range  of  temperatures. 

The  cylinder  is  divided  into  N  +  1  concentric  elements  as  shown 
in  Figure  8.  The  surface  element  has  a  thickness  of  AR/2,  where 
AR  =  R/N.  All  the  other  elements  are  AR  units  thick. 

Energy  balance  equations  are  constructed  for  the  various  elements. 

3T 

In  the  following  equations  -  (k  ^-)  is  evaluated  at  the  boundary 

between  element  i  and  element  i  +  1. 


CENTRAL  ELEMENT  (i=l) 


HEAT  IN  -  HEAT  OUT  =  ENERGY  ACCUMULATED 

9H-|  (T) 


0  +  2tt 


AR 

2 


/ |  8T \  _  ( AR \ 2 

(k  a?)  ih  ^  p  9t 


(3) 


where: 

k  =  thermal  conductivity* 
T  =  temperature, 
r  =  radial  posi tion  , 
p  =  density. 
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Surface  Element 


Element  I 


t^-AR/2 


FIGURE  8=  DIVISION  OF  CYLINDER  INTO 

CONCENTRIC  ELEMENTS 


H-j  =  enthalpy  per  unit  mass  of  element  1,  and 
t  =  time 
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Equation  (3)  is  rearranged  to  yield: 

,|8L  _  AR  3H1(T) 

^K9r'i+J$  4  p  9t 


(4) 


SURFACE  ELEMENT  (i=N+l) 

HEAT  IN  -  HEAT  OUT  =  ENERGY  ACCUMULATED 

-2,(R  -  f)  (k  f)  ^  -  h  2ttR  (Tn+1  - 

=  ptt[R2  -  (R  - 


V 


f)2] 


3WT> 

9t 


where: 

h  =  heat  transfer  coefficient  at  the  surface, 
T  =  environment  temperature,  and 

a 

T^+-|  =  surface  temperature. 

The  next  equation  follows  from  Equation  (5) 


?>  '  hR  (TN+1  *  V 


=  |  [  R  AR  - 


(AR) 


3HN+1 <T> 

at 


INTERIOR  ELEMENTS  (i:  2,N) 

HEAT  IN  -  HEAT  OUT  =  ENERGY  ACCUMULATED 

-  2 7i  (i-Us)  AR  (k  |I)  +  2 7i  (i-%)  AR  (k  |£)  i+H 

3H.(T) 

=  P77  [  { ( 1  -Js)  AR}2  -  AR  )2]  TTj: 


(5) 


(6) 


(7) 
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Equation  (7)  is  rearranged  to  produce 
8T\  \  /,.  3T< 


3  H  .  ( T ) 


(kff)  i_H  +  (1-%)  (k^)i+j5=  pAR(i-l)  -L 


(3) 


i  =2  ,N 


Equations  (4),  (6)  and  (8)  are  rewritten  in  terms  of  dimension¬ 
less  temperatures,  times  and  distances  to  yield: 


If)  1*2 
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Tn-T 
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3Hi (u) 
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(9) 


-  (1  - 


Ax  w  ,  3  U  \ 

T>  (k  s7> 


3x'  N+Jj  *  hR  UN+1 


pa 

2(VTa> 


[ax  -  Mr1] 


8Hn+i ^ 

3t 


(10) 


and 

'  (t-1*)  ‘k  If*  1-%  +  <k  If)  Uk 


apAX  i  3Hi^J) 

1VTT  (1‘1) 


3t 


i  =  2,N 


(id 


where 


u  = 


T-T 

a 

T  -T 
o  a 


and 


(12) 


Tq  =  initial  uniform  temperature  of  the  cylinder, 
x  =  r/R, 
ax  =  aR/R,  and 
t  =  at/R2 

FIRST  TIME  INCREMENT 


(13) 

(14) 

(15) 


In  the  first  time  increment,  (f^O-j+jg  is  rePlaced  in  equations 
(9)  -(11)  by  a  two  tier  time  level  approximation,  which  is: 
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C— ) 

^3x'  i+k 


1 

2 


1 

i+1 


-  u 


AX 


+ 


0 

'i+1 


AX 


(16) 


and  where: 

u°.  =  initial  dimensionless  temperature  at  the  centre  of  element  i, 
and 

u]  =  dimensioness  temperature  at  the  centre  of  element  i  at  the  end  of 
the  first  time  increment. 


From  the  definition  of  u  (Equation  12)  it  is  apparent  that: 


0  _  0  _  -i 

ui+l  ui 


(17) 


Equation  (16)  can  then  be  simplified  to  yield: 

1  1 

„  u .  .  -i  -  u . 

/BUn  _  1+1  1 

i+h 


2ax 


(18) 


The  --tj-  u--i  terms  are  replaced  by  a  forward  difference 

3  T 


approximation,  which  is 


3H(u). 

3t 


H ( u ) 1  -  H(u)? 

At 


(19) 


and  where: 

Ax  is  the  length  of  a  dimensionless  time  increment. 

When  the  finite  difference  approximations  are  substituted  into  Equations 
(9)  -(11)  the  following  equations  are  obtained. 


(20) 


o  o  a 


K  (uL  -  urli)  +  hR  K+l  +  UN+1  ) 

(ax) 


2- Ax 
2  Ax  "o  V“N+1 
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+  at(t  -tT  (ax  "  *  (h(u)n+i 


H( u )^+l )  =  0,  and  (21) 
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(i-Us)  k0  (u)  -  ul_1 )  -  (i-%)  kQ  (u]+1  -  u!) 


+  2ap^xli_0ijl  (H(u)l  _  H(u)0)=  °  1=2,N 

0  a 


(22) 


and  where: 

kQ  =  initial  value  of  thermal  conductivity. 


SUBSEQUENT  TIME  INCREMENTS 


g  u 

For  all  times,  after  the  first  time  increment,  (— )  .  x  is 
replaced  in  Equations  (9)  -(11)  by  a  three  tier  time  level  approximation 
For  time  increment  j,  the  approximation  is: 

1 


(— )J 
'  dx.  i+*i  3 


uJ+1  -  ul+1  u+3  -  uf1 
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The 


3H( u) . 
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approximati  on: 

3H(u) 


terms  are  replaced  by  a  forward  difference 
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Equations  (9)  -  (11) then  become: 


(24) 
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-  ( i_1)s)  [(uf1  -  uft)  +  (uJ  -  u^)  +  (uf1  -  u+])] 


(i  -  h)  k^  [(<]  -  uf1)  +  (UJ+1  -  uJ.)  +  (u^;]  -  uf1)] 


•  — •e4r4r^TT  (H<u>i+1  -  h(uW)  =  o  i=2,n 

v  o  a' 


(27) 


where: 


kj 

ki-^ 


kJ 

Ki+^ 


=  k  evaluated  at 


=  k  evaluated  at 


u'?  +  u3.  , 
i  i-l 


u'?  +  u^ 
i  l+l 


(28) 

(29) 


To  solve  the  Eouations  (20)  -  (22)  and  (25)  -  (27)  it  is 


necessary  to  know  how  the  enthalpy  changes  with  temperature.  As 

in  the  next  chapter,  H(T)  can  be  adequately  represented  by: 

expl ai ned 

H(T)  = 

a0  +  aiT 

T  >  Tc 

C1 

(30) 

H(T)  = 

exp  [bg  +  b^  T  +  b2T2  +  b^T3] 

T  <  T  <  T 

C2  C1 

(31) 

H(T)  = 

exp  [cQ  +  C]T] 

T<Tc, 

(32) 

where: 

aQ,  a^,  bQ,  b-j ,  b2»  by  cQ,  and  c-j  are  constants,  and 

T  =  initial  freezing  point  of  the  sample,  and 
C1 

T  =  temperature  at  which  In  H  vs  T  becomes  linear  (during  freezing) 

C2 

The  dependence  of  thermal  conductivity  on  temperature  is  also 
required.  The  k(T)  profile  is  approximated  by: 

T  >  T 


k(T)  =  dQ  +  d^ 
k(T)  =  eQ  +  e-j  T 
k(T)  =  fQ  +  fp  +  f2T2  +  fsT3 


U1 

Tc  ±T<T 

3 

T  <  T. 


(33) 

(34) 
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(35) 
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cIq,  ,  e^,  e^ ,  fg,  f-j  ,  and  are  constants,  and 

T  =  temperature  at  which  k  vs  T  is  no  longer  linear  (during  freezing). 
c3 

Prior  to  being  used  in  the  model  equations  the  H(T)  and  k(T)  functions 
are  first  converted  to  H(u)  and  k(u)  functions.  This  results  in  the 
following  set  of  equations: 


H(u)  =  a^  +  a^u 

H(u)  =  exp  (b^  +  b^u  +  b^u2  +  b^u3) 

H(u)  =  exp  (c 2  +C3U) 
and 


u  >  u  (36) 

C1 

u  <  u  <  u  (37) 

c2  “  ~  C1 

u  <  u  (38) 

c2 


k(u)  =  d^  +  d3u 
k(u)  =  e2  +  e3u 
k(u)  =  f4  +  f  5  u  +  f  6  u2  +  fyu3 

where : 
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The  Newton-Raphson  procedure  is  used  to  solve  the  non¬ 
linear  Equations  (20)  -  (22)  and  (25  )  -  (27  ). 

The  temperature  at  the  midpoint  of  an  element  is  used  in 
determining  its  enthalpy.  There  are  two  exceptions  to  this  rule.  In 
the  centre  element, the  temperature  at  the  midpoint  is  greater  than  in 
the  remainder  of  the  element  during  freezing  and  is  lower  during  thawing. 
Thus,  the  midpoint  temperature  cannot  be  considered  to  be  a  representative 
temperature  for  this  element.  Instead,  it  is  assumed  that  the  temperature 
varies  linearly  with  radial  position  between  the  midpoints  of  the  centre 


and  adjacent  elements.  It  then  follows  that  the  mass  average  temperature 

uJ  +  2u" 


in  the  centre  element  is 


1  .  This  temperature  is  used  in  the 


determination  of  the  enthalpy  of  the  centre  element.  In  calculating 


' 

1 
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the  enthalpy  of  the  surface  element  it  is  assumed  that  the  temperature 
at  the  inner  boundary  of  the  surface  element  is  (u^+-|  +  u^)/2.  The  average 
of  the  two  boundary  temperatures  of  the  surface  element  is  then 
(3u^+i  +  u^)/4.  This  temperature  is  used  to  calculate  the  enthalpy  of 
the  surface  element. 

A  computer  program  was  written  for  the  above  heat  transfer 
model.  It  is  listed  in  Appendix  6.  The  program  was  run  on  a  HEWLETT 
PACKARD  computer  (21  MXE-series).  Twenty  one  space  elements  (N=20)  were 
used  in  the  model  computer  runs  for  both  freezing  and  thawing  of  a  1  cm 
diameter  sugarbeet  cylinder.  For  freezing  a  time  increment  of  0.005  min 
was  used.  Under  these  conditions  9  hours  of  CPU  time  were  required  to 
execute  the  program.  However,  for  N=5,  the  CPU  time  was  only  1.5  hours. 
Reducing  N  from  20  to  5  changed  temperatures  at  a  given  time  by  less  than 
0.5  K.  For  thawing,  a  time  increment  of  0.1  min  was  used. 


10.  ESTIMATION  OF  THERMAL  PROPERTIES 


10.1  Literature  Review 

In  order  to  use  the  numerical  model  developed  in  the 
last  chapter,  it  is  necessary  to  know  the  dependence  of  both 
thermal  conductivity  and  enthalpy  on  temperature. 

Thermal  properties  of  various  foods  have  been 
experimentally  determined  by  a  number  of  authors 
(10,11,63,77,115,222,225) . 

Literature  values  for  thermal  conductivity  should  be 
carefully  examined  before  being  used.  Fennema  (77)  pointed 
out  a  number  of  inconsistencies  in  the  literature  values 
which  he  examined.  The  literature  on  thermal  conductivity 
and  enthalpy  will  now  be  reviewed. 


10.2  Thermal  Conductivity 

Cuevas  and  Cheryan  (53)  summarized  a  number  of 
different  methods  of  experimentally  measuring  or 
theoretically  predicting  the  thermal  conductivity  of  foods. 

Experimental  methods  for  the  measurement  of  thermal 
conductivity  are  basically  divided  into  two  classes: 


1.  Steady  state  methods  (e.g.  parallel  plate  and  concentric 
cylinders  or  spheres)  in  which  the  temperature  within 
the  food  is  kept  constant  with  time  (12,194). 

2.  Transient  methods  in  which  a  heated  probe  is  allowed  to 
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change  the  temperature  of  the  food  for  a  short  period  of 
time  (114,61). 


The  thermal  conductivity  may  alternatively  be  determined 
indirectly  from  experimentally  measured  diffusivity  (a) 
values  using  : 

k  =  pea 

A  number  of  simple  formulae  can  also  be  used  to 
estimate  the  thermal  conductivity.  Comini  (47)  obtained 
kj.  and  from  : 

k„=wki+(l-w)k  (61) 

j  w  w  d 

k.=wk+(l-w)k  (62) 

x,  w  w  w  d 

where  : 

w  =  mass  fraction  of  water. 

V7 

kp  =  thermal  conductivity  of  the  unfrozen  food, 

j/(s  .  m  .  K) 

k„  =  thermal  conductivity  of  the  frozen  food, 

J  J/( s  .  m  .  K) 

ki  =  thermal  conductivity  of  ice,  J/(s  .  m  .  K) 

k  =  thermal  conductivity  of  water,  J/(s  .  m  .  K) 

w 

k.  =  thermal  conductivity  of  dry  matter, 

d 

=  0.2602  J/(s  .  m  .  K) . 

In  the  freezing  range  k  is  assumed  to  vary  linearly  with 
temperature  from  k^  to  k^. . 
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Alternatively,  Fikiin's  (78)  formula  can  be  used  for 
temperatures  below  the  initial  freezing  point.  It  is: 


r  Tci  -  273.2“| 

k  ( T )  =  1 . 74 5w  1  -  -  +  0.233  (63) 

WL  T  -  273.2  J 


where : 

Tc !  =  initial  freezing  point  (K). 

By  assuming  that  the  thermal  conductivity  was 
linearly  dependent  on  the  percent  of  water  frozen, 
Schwartzberg  (184)  obtained  the  following  equation 
for  the  freezing  range: 


k (T )  =  k  +  ( k 

f  * 


-  k  )  ( Tw  Tc  i  )  / ( Tw  -  T) 

f 


(64) 


where : 

T  =  freezing  point  of  water. 

W 

He  concluded  that  this  equation  predicted  thermal 

conductivities  which  were  within  about  ±  10%  of 

most  experimental  values  he  examined  in  the  literature. 

Riedel's  empirically  obtained  formula  (176)  for 
thermal  diffusity  may  also  be  applied, 
for  an  unfrozen  food,  to  give: 


k  ( T )  =  6.978p(T)c(T)  (0.053+U  (T)-0.053)(w  )) 

w  w 


(65) 
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where : 

a  (T)  =  thermal  diffusity  of  water  at  temperature  T 
in  cm1  2 /min. 

Lentz  and  Van  den  Berg  (115,116,211)  developed  a  model 
for  predicting  the  thermal  conductivity  of  foodstuffs.  Their 
model  assumes  that  an  unfrozen  food  consists  of  nonaqueous 
components  dispersed,  as  spheres,  in  water.  They  assumed 
that  the  thermal  conductivity  of  frozen  foods  was  the 
average  of  two  k  values.  The  first  k  value  was  obtained  by 
assuming  water  was  the  continuous  phase  in  which  nonaqueous 
components  and  ice  were  dispersed.  The  second  k  value  was 
obtained  by  assuming  ice  was  the  continuous  phase.  They 
observed  that  their  calculated  thermal  conductivities  at 
253  K  were  within  ±15%  of  their  experimental  data. 


10.3  Enthalpy 

The  enthalpy-tempera ture  ( H ( T ) )  relationship  can  be 
experimentally  evaluated  (131,174)  or  theoretically 
predicted.  Two  theoretical  methods  will  now  be  discussed. 
Starting  with  the  freezing  point  depression  equation, 
Schwartzberg  (184)  derived  equations  for  the  calculation  of 
H(T)  by  making  the  following  assumptions  : 

1.  Enthalpies  of  food  components  are  additive  i.e.  heat  of 
solution  =  0. 

2.  The  solution  contained  in  the  food  is  ideal  i.e.  water 
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activity  (a  )  =  mole  fraction  of  water. 

w 

3 .  In  (a  )  =  a  -  1 . 

w  w 

For  lean  beef,  Schwar tzberg ' s  predicted  enthalpy  values 
which  deviated  by  about  5%  from  the  experimental  results  of 
Riedel . 

10.3.1  Riedel 1 s  Method 

Riedel  (174)  proposed  a  method  for  calculating  H(T)  for 
fruits  and  vegetables  during  freezing.  Riedel's  method  was 
used  in  this  thesis  and  thus  is  considered  here  in  some 
deta i 1 . 

Based  on  experimental  observation  Riedel  (174) 
concluded  that  the  specific  heat  of  vegetable  juices  (c  ) 
could  be  calculated  from  : 

c  =  4.186  -  2.386w0  (66) 

J 


where  : 

w0  =  mass  fraction  of  soluble  solids  in  unfrozen  juice. 
From  further  experimental  data  it  was  concluded  that  the 
mass  fraction  of  soluble  solids  in  the  residual  solution 
(w  )  during  freezing  is  related  to  the  temperature  by  : 

T  =  273.2  -  lOw  -  50w 3  (67) 

s  s 

Riedel  used  Randall's  formula  to  calculate  the  heat  of 
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fusion  of  ice  (L)  at  the  temperature  T. 


L  =  334.12  +  2.118 (T-273 . 2 )  -  0 . 00783 (T-273 . 2 ) 2 


(68) 


where : 

L  in  J/g,  and 
T  in  K. 

The  mass  fraction  of  the  juice  which  has  become  ice  (wj )  can 
be  calculated  from  : 


wi  =  1  -  w0/ws 


(69) 


Riedel  assumed  that,  in  warming  a  sample  from  T  to  293.2  K, 
all  the  ice  is  melted  at  T  and  the  unfrozen  solution  is 
warmed  to  293.2  K.  He  justified  this  assumption  on  the  basis 
that,  since  enthalpy  is  a  state  function,  the  enthalpy 
change  in  going  from  one  temperature  to  another  is 
independent  of  the  path.  Thus,  the  enthalpy  change  (AHj(T)) 
in  going  from  T  to  293.2  K  is  given  by  : 


AH  (T)  =  H  (293.2)  -  H  (T) 

=  wiL  +  c  (293.2  -  T) 


(70) 


where : 


H  (T)  =  enthalpy  (J/g)  of  juice  at  temperature  T. 

J 

Equations  (66)  and  (69)  are  substituted  into  Equation 


(70)  to  yield: 
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A  H  ( T )  =  (L+4. 186(293. 2-T) )-(L/w  +2 . 386 ( 293 . 2-T) ) w 0  (71) 

J  s 

To  obtain  the  H  (T)  profile  select  a  w  value.  Equation 

J  S 

(67)  determines  a  T  corresponding  to  this  w  .  Equation  (68) 

S 

then  calculates  an  L  value.  Using  these  w  ,  T  and  L  values 

S 

Equation  (71)  predicts  a  AHJ(T)  value.  By  selecting  a  number 
of  w  values  the  AH  (T)  profile  is  obtained.  Assuming  the 

S  J 

enthalpy  at  233.2  K  is  zero  then  the  H  (T)  profile  is 

J 

calculated  from  : 

H  (T)  =  AH  (233.2)  -AH  (T)  (72) 

J  J  J 

The  enthalpy  versus  temperature  profile  for  the  vegetable  is 
then  determined  from: 

H ( T )  =  (1-w  )H(T)  +  c  w  (T-233.2)  (73) 

u  J  u  u 


where  : 

H ( T )  =  enthalpy  (J/g)  of  the  vegetable  at  T. 

w  =  mass  fraction  of  insoluble  solids, 
u 

c  =  specific  heat  of  insoluble  solids 

u 

=  1.21  J/(g  .  K)  (value  suggested  by  Riedel) 
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11.  DETERMINATION  OF  THERMAL  PROPERTIES  OF  SUGARBEETS 


11.1  Estimation  of  k (T)  of  Suqarbeets 

Karlsruhe  Test  Substance  (175)  is  a  gel  containing  77% 
water,  methylcellulose ,  and  salt.  Sugarbeets  have  a  water 
content  of  78%  (133)  which  is  very  close  to  that  of 
Karlsruhe  Test  Substance.  In  fact,  a  sugarbeet  examined  in 
this  study  had  a  76%  water  content  (determined  using  a 
vacuum  oven  at  343  K  and  a  pressure  of  6.8  kPa).  The 
densities  of  Karlsruhe  Test  Substance  and  sugarbeets  are 
very  similar;  1.06  and  1.05  g/cm  ,  respectively.  Thus,  these 
two  materials  should  have  almost  equal  water  contents  on  a 
volume  as  well  as  on  a  weight  basis.  Van  den  Berg  and  Lentz 
(211)  observed  that  the  measured  thermal  conductivities  of  a 
number  of  foods,  with  the  same  volumetric  water  content,  at 
253  K  were  within  a  range  of  ±  15%.  As  the  thermal 
conductivity  of  a  food  is  primarily  dependent  on  water 
content  (211)  it  was  considered  reasonable  to  use  the 
published  k  versus  T  graph  (22)  for  Karlsruhe  Test  Substance 
to  represent  the  k(T)  profile  for  sugarbeets.  This  graph  was 
approximated  by  the  following  set  of  equations: 

k ( T )  =  0.5143  +  0.002199(T-273.2)  T  >  271.0  K  (74) 

k(T)  =  -0.7620  -  0 . 57657 (T-273 . 2 ) 

269.7  K  <  T  <  271.0  K  (75) 

k ( T )  =  0.9504  -  0 . 10697 (T-273 . 2 )  -  0 . 005726 (T-273 . 2 ) 2 

-0.0000984(T-273.2) 3  T  <  269.7  K  (76) 
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The  k(T)  equations  used  for  sugarbeets  were  not  quite 
identical  to  those  used  for  Karlsruhe  Test  Substance  (listed 
in  Appendix  F).  They  were  altered  slightly  to  allow  for  the 
fact  that  the  initial  freezing  point  of  sugarbeets  is  271  K, 
whereas  for  Karlsruhe  Test  Substance  it  is  272  K.  Thus,  the 
unfrozen  region,  with  its  lower  k(T)  values,  extends  for 
sugarbeets  to  271  K  instead  of  to  272  K.  Below  270  K  this 
modification  only  changes  the  k(T)  values  by  about  2%. 
However,  at  271  K  it  has  a  noticeable  effect  in  that  the 
k(T)  value  used  for  sugarbeets  is  0.510  J/(s  .  m  .  K) 
whereas  for  Karlsruhe  Test  Substance  it  is  equal  to 
0.935  J/(s  .  m  .  K) .  The  published  k(T)  values  for  Karlsruhe 
Test  Substance  and  the  predicted  k(T)  values  for  sugarbeets 
are  listed  in  Table  17.  Assuming  Bonacina  and  Comini's 
measured  k(T)  values  were  accurate  to  ±  15%,  the  k(T) 
values,  used  in  this  thesis  for  sugarbeets,  were  accurate  to 
±  30%. 
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TABLE  17:  k(T)  VALUES  FOR  KARLSRUHE  TEST  SUBSTANCE 

AND  SUGARBEETS 


TEMPERATURE 

THERMAL 

CONDUCTIVITY 

(K) 

(J/(s 

.  m  .  K )  ) 

KARLSRUHE 

1  SUGARBEETS 

243.2 

1.67 

1.67 

246.2 

1.65 

1.60 

248.2 

1.64 

1.59 

253.2 

1.61 

1.59 

257.2 

1.59 

1.60 

263.2 

1.53 

1.55 

265.2 

1.49 

1.50 

267.2 

1.44 

1.40 

268.2 

1.40 

1.35 

269.2 

1.29 

1.29 

270.2 

1.20 

0.97 

271.2 

0.85 

0.51 

272.2 

0.50 

0.51 

283.2 

0.53 

0.53 

293.2 

0.57 

0.57 

1  Read  from  a  published  k  versus  T  graph  (22). 

Approximate  values  based  on  the  published  graph 
for  Karlsruhe  Test  Substance. 
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11.2  Enthalpy  of  Sugarbeet s 

Riedel’s  method  was  used  to  predict  the 
enthalpy-temperature  relationship  for  sugarbeets.  The  H(T) 
values  obtained  were  checked  against  experimental 
measurements  made  by  personnel  of  the  Food  Science 
Department  of  the  University  of  Alberta.  The  experimental 
enthalpy  values  were  obtained  by  means  of  a  bomb  calorimeter 
(see  Figure  9).  A  small  sugarbeet  piece  ,  weighing  about  50 
mg,  was  placed  in  the  sample  holder.  The  calorimeter  was  put 
in  an  alcohol  bath  until  it  reached  the  bath  temperature  of 
213  K.  The  bath  was  then  warmed  at  a  rate  of  1  K/min. 

The  temperature  of  the  block  ( T b )  and  the  temperature 
difference  between  the  block  and  the  sample  (AT)  were 
continually  recorded  by  means  of  chr omel-alumel 
thermocouples  attached  to  a  Cole  Parmer  recorder. 

During  warming  the  rate  of  heat  flow  (Q(t))  from  the 
block  to  the  sample  and  sample  holder  is  given  by: 

Q(t)  =  B  A  T ( t )  (77) 

where : 

B  =  a  constant  which  is  obtained  from  a  calibration  run  with 
a  material  of  known  specific  heat. 

A  ’’blank”  run  is  performed  to  remove  the  effect  of  the 
sample  holder.  This  run  is  undertaken  without  a  sample  in 
the  holder  and  ATt(t)  is  measured.  The  rate  at  which  the 
sample  obtains  heat  during  a  normal  run  is  then  given  by: 
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FIGURE  9 :  CALORIMETER  FOR  MEASURING  H(T) 
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Q(t)  =  B  (aT  ( t )  -  ATi(t)) 


(78) 


However  , 
Q(t) 


me  (T) 
s  s 


dT 

dt 


(79) 


where : 

mg  =  mass  of  sample, 

c  (T)  =  effective  specific  heat  of  sample  (includes 
heat  of  fusion)  at  T,  and 
T  =  temperature  of  sample  =  Tb  -  AT. 

Equations  (78)  and  (79)  are  combined  to  yield: 

m  c  (T)dT  =  B (A  T ( t )  -  ATb(t))dt  “  (80) 

s  s 

Integration  of  this  equation  leads  to: 

T  t 

/  c  s  ( T )  dT  =  B_  /  (A  T  ( t )  -  ATb(t))dt  (81) 

m 

233.2  s  0 

or 

t 

H  ( T )  =  B_  /  (AT(t)  ~  A  T  b  ( t )  )  dt  (82) 


From  experiments,  using  acetic  acid  as  the  sample,  it  was 
determined  that: 

B  =  0.223  J/(min  .  K)  (83) 

A  blank  run  resulted  in  a  graph  of  ATt  versus  time.  Thus, 
the  experimentally  determined  values  of  AT(t)  and  ATb  (t) 
produced  (by  means  of  numerical  integration  of  Equation 


. 
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(82))  the  required  H(T)  data  for  the  sample. 

As  can  be  seen  from  Table  18  and  Figure  10,  results 
obtained  using  Riedel’s  method  agree  well  with  the 
experimental  data  below  the  initial  freezing  point  (271  K) . 
Above  this  temperature,  the  experimental  results  were 
apparently  in  error  as  they  produced  specific  heats  greater 
than  that  of  water  (4.2  J/(g  .  K) ) . 

To  facilitate  the  use  of  the  enthalpy  data  in  the 
model,  it  was  decided  to  approximate  the  H(T)  function  by 
mathematical  expressions  which  could  be  easily 
differentiated.  From  Figure  10  it  can  be  seen  that  the 
theoretical  H(T)  function  is  linear  above  the  initial 
freezing  point  but  has  a  non-linear  dependence  on  T  below 
it.  In  Figure  11,  In  H(T)  is  plotted  as  a  function  of  T. 
Based  on  this  figure,  In  H  was  represented  by  a  cubic 
function  of  T  between  271  K  and  266  K  and  by  a  linear 
function  of  T  below  266  K.  From  the  results  of  a  least 
squares  curve  fit,  in  the  region  266  K  to  271  K,  it  was 
determined  that  H(T)  could  be  approximated  by: 

H ( T )  =  321.98+3.695(T-273.2)  T  >  271  K  (84) 

H ( T )  =  exp[6. 7584+0. 6093(T-273. 2) 

+0.0775(T-273.2)2+0.00399(T-273.2)3] 

266  K  <  T  <  271  K  (85) 

H ( T )  =  exp[5.420+0.0722(T-273.2) ]  T  <  266  K  (86) 


. 


. 


108 


TABLE  18:  ENTHALPY  -  TEMPERATURE  DATA 


TEMPERATURE  ENTHALPY  (j/g) 


(K) 

RIEDEL1 

EXPERIMENTAL 

RUN  1 

RUN  2 

RUN  3 

RUN  4 

233.2 

0.0 

0.0 

0.0 

0.0 

0.0 

243.2 

23.4 

38.1 

30.1 

29.3 

25.5 

245.2 

31.0 

41.0 

36.4 

34.7 

31.0 

247.2 

33.9 

48.6 

43.1 

40.6 

37.3 

249.2 

40.2 

56.5 

50.7 

46.5 

43.5 

251.2 

46.9 

64.5 

58.6 

52.7 

50.2 

253.2 

53.6 

72.4 

67.0 

59.0 

56.9 

255.2 

61.5 

80.4 

76.6 

66.1 

63.6 

257.2 

70.3 

88.7 

87.1 

73.3 

69.5 

259.2 

80.0 

98.0 

98.0 

81.6 

76.2 

261.2 

91.7 

109.7 

109.3 

90.8 

85.8 

263.2 

106.3 

126.0 

122.6 

103.4 

104.2 

265.2 

125.6 

148.6 

140.2 

118.0 

121.4 

267.2 

154.0 

175.8 

162.8 

135.2 

136.8 

269.2 

220.2 

235.7 

212.6 

172.9 

181.7 

271.2 

316.0 

329.4 

303.1 

272.9 

285.1 

273.2 

323.2 

357.5 

360.8 

344.1 

316.0 

277.2 

337.8 

385.5 

367.1 

354.1 

351.6 

283.2 

359.6 

432.4 

416.1 

406.0 

391.4 

. 


• 

. 
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TEMPERATURE 

RIEDEL 

RUN  1 

RUN  2 

RUN  3 

RUN  4 

287.2 

374.2 

453.3 

435.8 

426.1 

412.3 

293.2 

396.0 

481.0 

460.9 

448.7 

439.9 

Calculated  using  Riedel's  method  (174). 
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FIGURE  10:  ENTHALPY  VS.  TEMPERATURE 
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FIGURE  11:  Ln  H  VS.  TEMPERRTURE 
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From  Figure  12  it  can  be  seen  that  the  above  equations 
adequately  represent  the  experimental  data  below  the  initial 
freezing  point. 


11.3  Determination  of  Initial  Freezing  Point 

At  the  initial  freezing  point  the  fraction  of  the  beet 
juice  which  has  become  ice  is  zero  i.e.: 

w  i  =  0  ( 87 ) 

Therefore  from  Equation  (69): 

w  =  w  0  ( 88 ) 

s 

and  Equation  (67)  becomes: 

Tci  =  273.2  -10wo  -  50w  0  3  (  89  ) 

where : 

Tcl  =  initial  freezing  point  of  beet  juice  (K). 

One  g  of  beet  is  composed  of  0.78  g  water,  0.18  g  soluble 
solids  and  0.04  g  insoluble  solids.  The  beet  juice  consists 
of  water  and  soluble  solids. 

Therefore , 


w  0  =  0.18/(0.78  +  0.18) 


0.1875 


(90) 


. 
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TEMPERATURE  (K) 


FIGURE  12:  PLOT  OF  H(T)  FROM  EQUATIONS 
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Thus , 


T  c  i  =  273.2  -  10  X  0.1875  -  50  X  0.18753 

=  271  K  (91) 

The  initial  freezing  point  was  also  measured 
experimentally  using  a  Differential  Thermal  Analyzer  and  was 
determined  to  be  271  K.  From  the  recorded  T(t)  curves  during 
air  thawing  of  beet  cylinders  it  was  estimated  to  be  between 
270.5  K  and  271  K. 


11.4  Heat  Transfer  Coef f ic ient  at  the  Surface 

The  determination  of  the  heat  transfer  coefficient 
involved  subjecting  a  copper  cylinder  to  the  same 
experimental  conditions  as  those  experienced  by  beet 
cylinders  during  freezing  and  thawing.  Analytical  solutions 
for  a  single  phase  material  were  then  used  to  calculate  the 
heat  transfer  coefficient  from  the  T(t)  data  which  was 
recorded  for  the  copper  cylinder. 

11.4.1  Experimental  Procedures 

The  ends  of  a  copper  cylinder  (15  cm  long,  1  cm 
diameter)  were  insulated  with  1  cm  long,  1  cm  diameter  cork 
cylinders.  A  sheathed  copper-constantan  thermocouple 
(THERMOELECTRIC  model  no.  T125G304-0-12-12-5H ,  wire  diameter 
=  0.13  cm,  sheath  diameter  =  0.102  cm)  was  placed  at  the 
geometric  centre  of  the  cylinder.  For  thermal  insulation, 
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the  part  of  the  thermocouple  sheath,  which  was  outside  the 
sample,  was  covered  with  0.64  cm  o.d.  POLYFLOW  tubing.  The 
thermocouple  lead  wire  was  insulated  with  asbestos  tape.  The 
voltage  produced  by  the  thermocouple  was  recorded  by  a 
Hewlett  Packard  recorder  (model  7100B),  maintained  at  a  full 
scale  range  of  5  mV.  Published  charts  showing  mV  versus 
temperature  for  a  copper-constantan  thermocouple  were  used 
to  prepare  a  table,  showing  temperature  as  a  function  of 
recorder  reading.  This  table  was  then  used  to  convert  the 
recorded  mV  readings  to  temperatures.  Measured  temperatures 
were  observed  to  be  accurate  to  ±  0.5  K. 

Imitating  the  situation  for  the  beet  cylinders 
(discussed  in  the  next  chapter)  the  middle  of  the  copper 
cylinder  was  covered  with  aluminum  foil.  The  ends  were 
covered  with  Scotch  plastic  tape.  A  schematic  diagram  of  the 
experimental  set-up  is  shown  in  Figure  13. 

Air  Warming 

The  copper  cylinder,  at  known  initial  temperatures  of  243  K, 
253  K,  or  263  K,  was  placed  on  the  sample  holder  (see  Figure 
14)  and  heated  by  ambient  air.  The  temperature  at  the 
cylinder  centre  was  measured  as  a  function  of  time. 

Cooling  in  Alcohol  Bath 

The  copper  cylinder,  at  a  temperature  of  approximately 
298  K,  was  placed  in  the  LT-50  alcohol  bath  which  was 
maintained  at  either  243  K,  253  K,  or  263  K.  The  variation 
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of  the  cylinder  centre  temperature  with  time  was  then 
recorded . 

11.4.2  Experimental  Results 

The  T(t)  data  obtained  for  air  warming  are  listed  in 
Tables  (E-l)  to  ( E  —  3 ) .  The  results  for  liquid  cooling  are 
summarized  in  Tables  (E-4)  to  (E-6). 

11.4.3  Calculation  of  h 


Air  Warming 

As  copper  has  a  very  high  thermal  conductivity,  the  internal 
resistance  to  heat  transfer  is  negligible  and  a  uniform 
temperature  is  maintained  in  the  cylinder.  For  this 
situation,  it  can  be  shown  that: 
m  c  T(t~At )-T 

h  =  -  In  - 2  (92) 

A  t  T  ( t )  -  T 

a 

where : 

m  =  mass  of  cylinder  =  101.4  g, 

c  =  specific  heat  of  copper  =  0.383  J/(g  .  K) , 

A  =  surface  area  =  nDL  =  0.004668  m2 , 

At  =  time  difference  between  two  temperature 
measurements  ( s ) , 

T(t)  =  temperature  of  copper  at  time  t  (K), 

Ta  =  temperature  of  cooling  or  heating  medium  (K) . 
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During  air  warming  of  the  copper  cylinder,  Biot  numbers 
(hR/k)  of  the  order  of  0.0001  were  obtained.  Consequently, 
Equation  (92)  is  applicable  for  air  thawing.  Using  the  data 
in  Tables  (E-l)  to  (E-3),  the  variation  of  h  with  time  is 
calculated  from  Equation  (92)  and  is  shown  in  Tables  (E-7) 
to  ( E—  9 ) .  The  value  of  the  dimensionless  temperature  at  the 
surface  (uN+1)  where: 

T(t )-Ta 

uN  +  !  = -  (  93  ) 

T0  -  Ta 

is  calculated  at  the  various  times  and  is  also  listed  in 
Tables  (E-7)  to  (E-9).  The  data  in  Tables  (E-7)  to  (E-9)  are 
plotted  in  Figures  15  -  17  in  the  form  of  h  versus  uN + x . 
Least-squares  curve  fitting  of  the  h(uN+1)  data  produced  the 
following  second  degree  approximations  for  the  heat  transfer 
coef  f ic ient : 


To  =  243  K 


h  =  0.8953  +  49.3939un+1  -  30 . 5409 ( uN+ x ) 2  (94) 


To  =  253  K 


h  =  0.0363  =  39.0551un+1  -  19. 3128 (uN+ x ) 2  (95) 


To  =  263  K 


-  9 . 1 4  7  9  (  uN  +  !  )  2 


h  =  1.9937  +  26.9908un+  x 


(96) 
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N  +1 


FIGURE  15  :  h  VERSUS  u 


(T0  =243  K) 
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FIGURE  16  :  h  VERSUS  uN+1  (T0  =253  K) 


, 


N+1  (T0  =263  K) 


FIGURE  17  :  h  VERSUS  u 
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Error  Analysis  of  h  Used  in  Air  Thawing 
The  measured  temperature  of  the  copper  cylinder  is 
accurate  to  ±  0.5  K.  For  0.25  <  uN+1  <  1.0,  this  leads  to  a 
maximum  error  of  approximately  10%  in  the  heat  transfer 
coefficient,  calculated  from  Equation  (92).  The  maximum 
percentage  error  in  h  increases  as  uN+1  is  further  reduced. 
This  is  because  T(t)  -  T  and  T(t  -  At)  -  T  in  Equation 

El  El 

(92)  are  decreasing  with  uN+1  and  are  affected  to  a  greater 
degree  by  an  error  in  the  temperature  measurement.  When  uN  +  i 
reaches  about  0.1  the  maximum  error  in  h  will  have  risen  to 
about  50%.  However,  for  uN * x  values  less  than  0.25  an 
accurate  estimate  of  h  is  not  very  important  for  air 
thawing.  This  is  because  the  sample  temperature  is  above  283 
K  and  has  already  been  thawed. 
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Cooling  in  Alcohol  Bath 

The  heat  transfer  coefficient  is  now  evaluated  for  cooling 
in  the  LT-50  bath.  The  alcohol  in  this  bath  is  continually 
circulated.  Consequently,  forced  convection  is  the  primary 
method  of  heat  transfer  between  the  cylinder  and  the 
alcohol.  It  is,  therefore,  assumed  that  the  heat  transfer 
coefficient  remains  constant  during  cooling.  The  h  values 
are  calculated  from  Carslaw  and  Jaeger's  (35)  analytical 
solution  for  one  phase  heat  transfer  in  an  infinite 
cylinder.  This  solution  is: 

-(£  ) 2 at /R 2 

T  ( r ,  t )  -  T  n  J,(r«  /E) 

- §/  =  2P>  - 2 -  (97) 

T.  -  Ta  ^  ( (0n) 2+P2 )J„ (0n) 


where : 

±R  ,  n  =  1,2, - are  the  roots  of 

n 

(3  Jx  (f3  )  =  PJo  (0  )  (98) 

n  n  n 

J0(/8  )  and  J  x  ( >8  )  are  zero  and  first  order 
n  n 

Bessel  functions,  and 

P  =  Rh/k  (99) 


A  trial  and  error  procedure  is  used  to  calculate  h  from 
Equation  (97).  Different  values  of  h  are  substituted  into 
this  equation  until  the  calculated  and  measured  values  of 
T(0.5)  are  equal.  The  temperature  after  0.5  min  was  selected 
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because  at  least  50%  of  the  temperature  change  inflicted  on 
the  copper  cylinder  occurs  during  the  first  0.5  min.  It  can 
thus  be  considered  to  yield  an  "average"  h  value  for  the 
cooling  process.  For  times  greater  than  0.5  min,  the 
calculated  heat  transfer  coefficient  becomes  more  sensitive 
to  errors  in  the  measured  temperature. 

For  a  constant  heat  transfer  coefficient.  Equation  (97) 
yields  more  accurate  h  values  than  Equation  (92).  This  is 
because  it  allows  for  a  nonuniform  temperature  profile 
inside  the  copper  cylinder.  Equation  (97)  was  not  used  for 
air  thawing  because  the  heat  transfer  coefficient  changed 
dramatically  with  uN+1. 

The  second  computer  program  in  Appendix  H  was  used  to 
determine  h.  The  calculated  heat  transfer  coefficients  are 
listed  in  Table  19.  The  mean  values  in  this  table  are  used 
as  the  heat  transfer  coefficients  for  freezing  of  sugarbeet 
cylinders  in  the  LT-50  bath.  They  are  summarized  in  Table 
20. 


Error  Analys i s  of  h  used  for  Alcohol  Bath 
The  h  values  for  cooling  were  also  calculated,  as  a 
function  of  time,  using  Equation  (92).  These  values  are 
listed  in  Tables  (E-10)  to  (E-12).  From  these  values,  it  is 
apparent  that  errors  in  h  are  of  the  order  of  50%  and  h 
generally  decreases  with  uN+1.  Equation  (92)  was  used  to 
estimate  the  error  involved  in  measuring  h  because  the 
variation  of  h  with  time  could  be  easily  obtained  from  this 
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TABLE  19: 

CALCULATED  h 

VALUES  1 

FOR 

THE  FIRST 

0.5  min  OF 

COOLING 

T 

RUN  1 

RUN  ; 

l  RUN  3 

MEAN 

S.D. 

<* 

2  S.D. /MEAN 

(K) 

243 

238 

283 

328 

283 

45 

0.159 

253 

285 

317 

370 

324 

43 

0.133 

263 

432 

377 

423 

411 

30 

0.072 

1  The  h  values  (J/s  .  m2  .  K)  were  calculated 
using  Carslaw  and  Jaeger’s  analytical  solution. 


2  S.D.  =  standard  deviation. 


TABLE  20:  HEAT  TRANSFER  COEFFICIENTS  FOR  FREEZING 


IN  THE  LT-50  BATH 


BATH  TEMPERATURE  h  (J/s  .  m2  .  K) 


243 

K 

283 

253 

K 

324 

263  K 


411 


i* » 
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equation.  This  variation  gives  some  indication  of  the 
accuracy  involved  in  the  h  measurements.  For  cooling  of  the 
copper  cylinder  in  the  alcohol  bath,  Biot  numbers  of  the 
order  of  0.001  were  obtained.  Thus,  for  constant  h  and  a 
given  T0 ,  T  and  At,  Equation  (92)  predicted  a  h  value 

d. 

within  5%  of  that  calculated  from  Equation  (97). 
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12.  EVALUATION  OF  FREEZING  MODEL 


12.1  ANALYTICAL  SOLUTION 

The  numerical  model  was  first  tested  for  a 
one-dimensional  heat  transfer  problem.  The  thermal 
conductivity  and  specific  heat  were  assumed  to  be  constant. 
An  analytical  solution  (Equation  (97))  is  given  in  Carslaw 
and  Jaeger  (35)  for  this  problem.  The  second  computer 
program  in  Appendix  H  was  used  to  calculate  the  T(t)  profile 
from  the  analytical  solution.  As  shown  in  Table  21,  the 
model  results  are  within  0.1  K  of  the  analytical  solution. 


12 . 2  KARLSRUHE  TEST  SUBSTANCE 

The  numerical  model  was  next  tested  against 
experimental  data,  obtained  by  Cleland  and  Earle  (44),  for 
cylinders  of  Karlsruhe  Test  Substance.  This  material  is  a 
methylcellulose  gel  which  has  been  used  (22,41-44)  for 
testing  food  freezing  models.  Cleland  and  Earle  considered 
the  "freezing  time"  to  be  the  time  required  for  the  centre 
of  a  food  to  reach  263  K.  They  (44)  published  experimental 
data  for  the  freezing  time  of  a  5.2  cm  diameter  cylinder  of 
Karlsruhe  Test  Substance  under  various  freezing  conditions. 
This  data  consisted  of  freezing  times  which  corresponded  to 

specified  h,  T0 ,  and  T  values.  The  h  values  were  reported 

a 

to  be  correct  to  within  ±  4%. 
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TABLE  21:  RESULTS  FOR  ONE-PHASE  HEAT  TRANSFER 


TEMPERATURE  (K) 


TIME 

SURFACE 

r/R 

=  0.5 

CENTRE 

(min ) 

ANAL. 

MODEL 

ANAL. 

MODEL 

ANAL. 

MODEL 

0.375 

274.1 

274.1 

287.7 

287.6 

291.5 

291.5 

0.75 

265.8 

265.8 

279.6 

279.5 

284.4 

284.4 

1.125 

261.3 

261.3 

272.7 

272.6 

276.9 

276.9 

LO 

• 

1 - 1 

257.8 

257.8 

267.0 

267.0 

270.5 

270.5 

1.875 

255.0 

255.0 

262.5 

262.4 

265.3 

265.3 

2.25 

252.8 

252.7 

258.8 

258.7 

261.1 

261.0 

2.625 

250.9 

250.9 

255.8 

255.8 

257.6 

257.6 

CONDITIONS 

h 

=  293  J/(s  .  m2  .  K) 

T0  =  293.2  K 

c 

s 

=  8.37  j/(g  .  K) 

R  =  0 . 005  m 

k 

=  1.0  J/(s  .  m  .  K) 

p  =  1.05  X  106  g/m3 

T 

=  243.2  K 

At  =  0.225  s 

a 


N  =  20 


. 
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In  this  thesis  published  H ( T )  and  k(T)  data,  for 
Karlsruhe  Test  Substance,  were  put  in  a  form  suitable  for 
use  in  the  numerical  model.  More  details  are  given  in 
Appendix  F.  Cleland  and  Earle’s  h,  T0 ,  and  T  values  were 

cl 

then  entered  and  a  freezing  time  calculated.  The  model  and 
experimental  results  are  both  listed  in  Table  22.  It  can  be 
seen  that  the  model  freezing  times  deviated  from  the  average 
experimental  value  by  between  -2  and  +11.6%.  This  is 
essentially  the  same  deviation  as  reported  by  Cleland  and 
Earle  (44)  for  the  best  finite  difference  scheme  which  they 
tested . 

The  model  predicts  freezing  times  which  are  generally 
around  the  upper  limit  of  experimental  values.  It  thus 
predicts  freezing  times  which  err  on  the  safe  side  i.e.  it 
ensures  nearly  all  the  samples  are  frozen  to  263  K,  rather 
than  just  the  average  sample.  However,  in  one  case  (T0=298.5 

K,  T  =243  K,  h=23.9  J/(s  .  m2  .  K))  the  model  actually 

a 

predicted  a  freezing  time  2.5  min  shorter  than  the  average 
experimental  value.  This  indicates  that  there  is  not  a 
fundamental  error  in  the  model  which  would  always  make  it 
predict  freezing  times  longer  than  those  which  actually 


occur . 
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TABLE  22:  FREEZING  TIMES  FOR  KARLSRUHE  TEST  SUBSTANCE 


RUN  CONDITIONS  FREEZING  TIME  (min)1 


EXPERIMENTAL2 

MODEL 

%  DIFF 

T  e  =  3  01 . 2  K  ; 
h=35 . 5  J/(s 

T  =253.3 
.  m2  .  K) 

K 

129 . 6±10 . 2 

140.9 

8.7 

T0  =292 . 9 
h=  3  5 . 5 

K  ; 

Ta  =253.2 

d 

K 

118 . 8±9 . 6 

132.6 

11.6 

T  o  =  290 . 4 
h=35 . 5 

K  ? 

T  =240 . 0 

d 

K 

70 . 8±4 . 8 

77.0 

8.8 

T0  =  298 . 5 
h=23 . 9 

K  ? 

T  =243.0 

a 

K 

126 . 0±11 . 4 

123.5 

-2.0 

T0  =  288 . 6 
h=23 . 9 

k  ; 

T  =233.0 

a 

K 

80 . 4±7 . 8 

85.2 

6.0 

T0  =  277 . 5 
h=35 . 5 

K  ; 

T  =239.7 

a 

K 

63 . 6±4 . 2 

65.6 

3.1 

T0  =  277 . 2 
h=35 . 5 

K  ; 

T  =239.5 

a 

K 

62 . 4±4 . 2 

65.0 

4.2 

1  FREEZING  TIME  =  time  for  centre  of  sample  to  reach  263  K 

2  Cleland  and  Earle's  experimental  results 


MODEL  RESULT-MEAN  EXPERIMENTAL  VALUE 


%  DIFF 


MEAN  EXPERIMENTAL  VALUE 


X  100% 
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12 . 3  FREEZING  AND  THAWING  OF  SUGARBEETS 

The  model  was  also  tested  against  experimental  results 
obtained  by  the  author  for  1  cm  diameter,  15  cm  long 
sugarbeet  cylinders. 

12.3.1  Experimental  Procedures  Used  to  Test  Model 

The  temperature  at  the  geometric  centre  of  a  sugarbeet 
cylinder  (1  cm  diameter,  15  cm  long),  insulated  at  both  ends 
by  cork  (1  cm  diameter,  1  cm  long),  was  measured  as  a 
function  of  time  during  freezing  and  thawing.  The  15  cm  long 
cylinder  was  obtained  by  glueing  together  three  5  cm  long 
cylinders.  This  was  necessary  because  of  the  small  size  of 
the  beets.  The  middle  cylinder  was  covered  with  two  layers 
of  aluminum  foil  to  prevent  dehydration.  For  ease  of 
handling,  the  two  outside  cylinders  and  cork  insulation  were 
covered  with  plastic  tape  instead  of  aluminum  foil. 

A  sheated  copper-constantan  thermocouple  (wire  diameter 
0.013  cm,  sheath  diameter  0.102  cm)  was  used  to  measure  the 
temperature  at  the  geometric  centre  (see  Figure  18).  Heat 
transfer  from  the  thermocouple  wire  to  the  surrounding  fluid 
was  minimized  by  covering  the  sheath  with  0.635  cm  (1/4 
inch)  polyflo  tubing.  The  thermocouple  lead  wire  was 
insulated  with  asbestos  tape. 

In  order  to  test  the  model  results,  the  sugarbeet 
cylinder  was  frozen  in  the  LT-50  bath,  which  was  maintained 
at  a  constant  temperature  of  either  243  K,  253  K,  or  263  K, 


I  I  ; 
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FIGURE  I8:  EXPERIMENTAL  ARRANGEMENT 

FOR  TESTING  MODEL 


■  \  H.  P.PA  JATM3MW31X3  81  IflUOR 
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and  the  centre  temperature  recorded  as  a  function  of  time. 
The  sugarbeet  cylinder  was  later  thawed  in  still  air,  by 
natural  convection.  For  still  air  thawing,  the  cylinder  was 
placed  on  a  wooden  holder  (Figure  14).  This  holder  was 
located  beside  the  LT-50  bath  so  the  frozen  sample  could  be 
quickly  placed  on  the  holder.  The  temperature  of  the  ambient 
air  was  measured  with  a  thermometer,  accurate  to  ±  0.5  K. 
During  both  freezing  and  thawing  the  cylinder  was  in  a 
horizontal  position. 

The  experimental  results  and  model  predictions  are 
listed,  in  Appendix  G,  for  various  freezing  and  thawing 
conditions . 

The  freezing  times  (time  for  centre  of  sample  to  reach 
263  K)  in  the  243  K  and  253  K  alcohol  baths  are  summarized 
in  Table  23.  For  the  243  K  bath  the  average  of  the  three 
experimental  freezing  times  was  2.2  min,  whereas  the  model 
predicted  a  freezing  time  of  2.5  min.  In  contrast,  the 
predicted  freezing  time  (3.5  min)  in  the  253  K  bath  was  less 
than  the  average  experimental  value  (4.1  min). 
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TABLE  23:  FREEZING  TIMES  IN  243  R  AND  253  K  BATHS 


T  =  243  K 

a. 


FREEZING  TIME  (min) 1 


RUN  1 

2.1 

RUN  2 

2.1 

RUN  3 

2.3 

MODEL  ( h=283 ) 2 

2.5 

T  =  253  K 

a 


FREEZING  TIME  (min) 


RUN  1 

4.1 

RUN  2 

4.0 

RUN  3 

4.0 

RUN  4 

4.2 

MODEL  ( h=317 ) 

3.5 

1  FREEZING  TIME  =  time  for  centre  of  sample  to  reach  263  K. 


2 


h  is  in  units  of  J/(s  .  m2  .  K) 


« 
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Another  indication  of  how  close  the  predicted 
temperatures  are  to  experimental  values  is  given  by  the 
’’standard  deviation  from  the  predicted  curve".  It  is  defined 
by : 


standard  deviation 


(t) 

(n 


-  Ve  <*)]» 

-  1) 


(98) 


where : 

n  =  number  of  time  increments  at  which  the  temperatures 
are  compared, 

Texp  (t)  =  average  experimental  value  at  time  t,  and 
Tpre  (t)  =  predicted  value  at  time  t. 

The  denominator  is  n  -  1  because  one  degree  of  freedom  is 

lost  as  a  result  of  setting  T  (0)  equal  to  T  (0).  The 

pre  ^  exp 

standard  deviations  were  calculated  for  the  different 
freezing  and  thawing  conditions  from  the  data  in  Tables  G-l 
to  G-6.  They  are  summarized  in  Table  24  along  with  the 
maximum  differences  between  experimental  and  predicted 
temperatures . 
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TABLE  24:  COMPARISON  OF  EXPERIMENTAL  AND  PREDICTED  RESULTS 


FREEZING 

MAXI MUM 1  S  T ANDARD 2 

CONDITIONS 

DIFFERENCE  DEVIATION 

(K)  (K) 

T  =  243  K 

9.2  4.1 

T  =  253'  K 

4.8  2.2 

T  =  263  K 

3.2  1.9 

THAWING 

MAXIMUM  STANDARD 

CONDITIONS 

DIFFERENCE  DEVIATION 

(K)  (K) 

T0  =  243  K 

6.3  2.4 

T0  =  253  K 

4.2  1.3 

To  =  263  K 

4.9  1.6 

1  MAXIMUM  DIFFERENCE  is  the  maximum  value 

Of  |T  (t)  -  T  (t)  I  . 

1  exp  pre 


2 

OT  (t)  -  T  (t)]2 

/  exp  pre 

STANDARD  DEVIATION 


( n  -  1 ) 


' 
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The  tables  in  Appendix  G  are  summarized  in  Figures  19  - 
26.  For  freezing  in  the  243  K  bath,  Figure  19  shows  that  the 
model  predicts  the  correct  general  trend  but,  for 
h  =  283  J/(s  .  m2  .  K) ,  lags  the  experimental  data  by  about 
0.5  min,  after  the  initial  freezing  point  has  been  reached. 
It  was  thought  that  the  difference  between  predicted  and 
experimental  curves  might  result  from  an  error  in  the 
determination  of  h.  Consequently,  an  h  value  of  417  was  also 
used.  From  the  dashed  line  in  Figure  19  it  is  apparent  that 
increasing  the  value  of  h  actually  caused  the  predicted 
results  to  precede  the  experimental  values  by  about  0.5  min. 
Thus,  an  incorrect  h  value  could  be  responsible  for  the 
differences  between  predicted  and  experimental  results. 

The  results  for  the  253  K  bath  are  somewhat  improved. 
After  the  initial  freezing  point  is  reached,  the  theoretical 
curve,  in  Figure  20,  initially  lags  the  experimental  curve 
but  the  situation  is  later  reversed.  This  may  be  due  to  the 
fact  that  the  actual  h  value  is  decreasing  whereas  a 
constant  h  is  used  in  the  model. 

The  experimental  data  in  Figure  21  appears  to  suggest 
the  existence  of  two  plateaus,  a  short  one  at  270.5  K  and  a 
longer  one  at  269  K.  Pereira  (167)  discussed  this  phenomenon 
in  some  detail.  Numerous  cells  are  damaged  when  a 
thermocouple  is  placed  inside  a  solid  food.  Liquid  then 
exudes  from  these  cells  and  surrounds  the  thermocouple 
junction.  According  to  Maximov’s  "wound  sap  theory", 
freezing  of  this  liquid  is  responsible  for  the  first 
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FIGURE  19:  FREEZING  IN  243  K  BRTH 
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FIGURE  20:  FREEZING  IN  253  K  BRTFI 
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FIGURE  21:  FREEZING  IN  263  K  BRTFi 
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FIGURE  22:  EFFECT  OF  CHANGING  k 
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FIGURE  23:  AIR  THAWING  -  EFFECT  OF  CONSTANT  h 
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FIGURE  24:  RIR  THRUING  OF  243  K  BEET 
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FIGURE  25:  BIR  THRWING  OF  253  K  BEET 
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FIGURE  26: 


RIR  THRUING  OF  263  K  BEET 
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plateau.  Pereira  supported  Maximov's  theory.  His 
experimental  results  also  indicated  that  the  second  plateau 
is  caused  by  freezing  within  individual  cells  close  to  the 
thermocouple.  The  fact  that  the  second  plateau  is  about  2  K 
below  the  initial  freezing  point  may  be  a  result  of 
supercooling . 

To  estimate  how  sensitive  the  model  is  to  errors  in 
k(T),  it  was  decided  to  increase  all  the  k(T)  values  by  20 
or  100%  and  then  examine  the  changes  in  the  model  results. 
From  Figure  22  it  can  be  seen  that  a  20%  increase  in  k(T) 
values  reduces  the  time  to  reach  a  given  temperature  by  less 
than  0.2  min.  Figure  22  shows  that  even  a  100%  increase  in 
k(T)  does  not  prevent  a  plateau  occurring  in  the  model  T(t) 
profile.  The  existence  of  this  plateau,  in  the  model 
predictions,  may  result  from  the  use  of  an  incorrect  heat 
transfer  coefficient.  Alternatively,  the  lack  of  a  plateau 
in  experimental  results  may  be  caused  by  supercooling. 

The  model  is  based  on  equilibrium  enthalpy  values  and 
thus  does  not  allow  for  supercooling.  Initially,  an  attempt 
was  made  to  consider  supercooling  but  then  it  was  realized 
that  different  parts  of  the  cylinder  would  have  different 
degrees  of  supercooling.  To  predict  the  amount  of 
supercooling  throughout  the  cylinder  it  would  be  necessary 
to  allow  for  water  diffusion.  This  would  be  a  possible 
future  extension  of  the  model.  The  presence  of  supercooling 
might  explain  why  the  model  results  lag  behind  the 
experimental  values  just  below  the  initial  freezing  point. 
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As  heat  of  fusion  need  not  be  removed,  supercooled  foods 
will  initially  cool  faster  than  predicted. 

Originally,  in  testing  the  model  for  still  air  thawing, 
a  constant  h  value  of  11  J/(s  .  m2  .  K)  was  used.  This 
estimate  of  h  was  suggested  by  Charm  (37)  for  slow  moving 
air.  From  Figure  23,  it  is  apparent  that  use  of  this  h  value 
resulted  in  poor  agreement  being  obtained  between  model 
predictions  and  experimental  results.  The  experimental  h 
values  in  Equations  (33)  -  (35)  were  then  employed.  As  can 
be  seen  from  Figures  24  -  26,  this  led  to  excellent 
agreement  between  experimental  and  predicted  temperature 
prof i les . 

The  experimental  h(uN+1)  data,  for  air  thawing,  was 
also  approximated  by  cubic  and  discontinuous  functions.  Use 
of  these  functions  resulted  in  T(t)  curves  which  were  as 
close  to  experimental  values  as  the  model  results  in  Figures 
24  -  26.  It  is,  therefore,  concluded  that  many  different 
types  of  approximating  functions  can  be  used  to  represent 
h(uN+1)  and  produce  accurate  T(t)  profiles. 

Model  results  (for  Tc  =  243  K,  h  =  283  J/(s  .  m2  .  K)) 
are  used  in  Figure  27  to  show  how  the  T(t)  profile  changes 
with  position,  during  freezing.  This  figure  indicates  that  a 
very  fast  freezing  rate  is  obtained  at  the  centre  relative 
to  the  outer  regions.  As  TN+1-  T^  at  the  surface  is  small 
when  the  centre  is  being  frozen  one  might  expect  the  latter 
to  freeze  slowly.  However,  the  model  predictions  agree  with 
the  experimental  results  of  Meryman  (142)  and  Luyet  (129). 
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FIGURE  27:  T(t)  VARIATION  UITH  POSITION 
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Fast  freezing  rates  are  obtained  at  the  centre  because  it  is 
easier  to  remove  the  heat  of  fusion  when  the  centre  starts 
to  freeze.  At  this  stage  most  of  the  heat  of  fusion  has 
already  been  removed  from  the  exterior  regions.  Also,  the 
thermal  conductivity  has  been  significantly  increased  as  a 
result  of  freezing. 

In  conclusion,  the  numerical  model  predicted 
temperature-t ime  relationships  which  were  within  0.1  K  of 
the  analytical  solution  for  one  phase  heat  transfer.  For 
Karlsruhe  Test  Substance,  the  model  predicted  freezing  times 
which  were  within  -2  to  +11.6%  of  experimental  values.  For 
freezing  of  sugarbeet  cylinders  in  the  243  K  bath,  the 
predicted  freezing  time  was  13.6%  greater  than  the  average 
experimental  value.  The  situation  was  reversed  for  the  253  K 
bath.  In  this  bath,  the  predicted  freezing  time  was  14.6% 
lower  than  the  average  experimental  value.  The  maximum 
difference  between  experimental  and  predicted  temperatures 
was  9.2  K  for  freezing  and  6.3  K  for  thawing.  The  greatest 
value  of  the  standard  deviation  from  the  predicted  curve  was 
4.1  K  for  freezing  and  2.4  K  for  thawing.  The  differences 
between  model  predictions  and  experimental  values ^could  have 
resulted  from  the  inaccuracies  involved  in  the  determination 
of  the  thermal  properties  or  from  failure  to  allow  for 
supercooling . 
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13.  CONCLUSIONS 


1. 


2. 


3. 


4. 


Frozen  beets  should  not  be  allowed  to  thaw  prior  to 
processing.  This  is  because  thawed  beets  are  soft  and 
decay  rapidly. 
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Two  frozen  storage  methods  appear  economically 
attractive  and  warrant  further  investigation.  The  first 
one  is  the  storage  of  frozen  beets  in  outdoor  piles.  The 
second  method  is  the  frozen  storage  of  beets  in 
mechanically  cooled  buildings. 

Use  of  enthalpy,  in  the  energy  balance  equations, 
simplifies  the  development  of  heat  transfer  models  for 
freezing  and  thawing  of  biological  materials.  No  loss  in 
accuracy  results  from  this  simplification. 
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1 4 .  RECOMMENDATI ONS 


1 •  A  numerical  model  should  be  developed  for  the  freezing 
and  thawing  of  storage  piles  of  sugarbeets  or  other 
plant  material.  This  model  could  then  be  used  to  predict 
the  effect  of  unusual  weather  conditions,  such  as 
Chinook  winds  in  Southern  Alberta,  on  the  temperature  of 
stored  beets. 

2.  Frozen  storage  experiments  should  be  performed  on  beet 
piles  in  the  Taber  area.  Cold  night  air  should  be  blown 
through  the  pile  to  freeze  the  beets  to  263  K  or  below. 
Temperatures  in  the  pile  should  then  be  monitored  to 
determine  how  long  the  beets  remain  frozen.  The  use  of 
pile  covers  in  the  prevention  of  thawing  should  also  be 
investigated . 

3.  Accurate  experimental  measurements  of  thermal 
conductivity  and  enthalpy,  in  the  freezing  range,  should 
be  obtained  for  numerous  biological  materials. 

4.  The  enthalpy  approach  should  be  used  to  develop  heat 
transfer  models  for  other  shapes,  such  as  finite 
cylinders,  finite  slabs  and  cones. 

5.  The  numerical  model,  developed  in  this  thesis,  should  be 
extended  to  allow  for  the  effects  of  both  supercooling 
and  water  diffusion. 
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15.  NOMENCLATURE 


2c  f3i  1^2  /S] 

constants  used  in  calculating  enthalpy 

a 

W 

water  activity 

A 

surface  area,  cm2 

bo  ,  .  .  ( b  i 

constants  used  in  calculating  enthalpy 

B 

constant  used  in  the  experimental 
determination  of  enthalpy,  J/(min  .  K) 

c 

specific  heat,  J/(g  .  K) 

CJ 

specific  heat  of  vegetable  juice, 

J/(g  .  K) 

c 

s 

specific  heat  of  sample,  J/(g  .  K) 

c 

u 

specific  heat  of  insoluble  solids, 

J/(g  .  k) 

Cc,C:,C2,C3 

constants  used  in  calculating  enthalpy 

do ,dlfd2fd3 

constants  used  in  calculating  thermal 
conductivity 

D 

shape  factor  used  in  Plank's  equation 

e  0  r  2 i  f  2 2  r63 

constants  used  in  calculating  thermal 
conductivity 

EMS 

error  mean  square 

^  o  f  f  i  r  •  •  r  ^  7 

constants  used  in  calculating  thermal 
conductivity 

G 

shape  factor  used  in  Plank's  equation 

H 

enthalpy  per  unit  mass,  J/g 

H 

J 

h 

enthalpy  per  unit  mass  of  juice,  J/g 

heat  transfer  coefficient  at  the  surface, 
J/(s  .  m2  .  K) 

J.(  ) 

zero  order  Bessel  function 

Jx(  ) 

first  order  Bessel  function 

k 

thermal  conductivity,  J/(s  .  m  .  K) 
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k 


d 


k  i 


k 

w 


1 


LSR 

m 

m 

s 

n 


N 

P 

Q 

r 

rep 

R 

s 

m 

Sail 
Sr  1 0 

St  8 

SSR 

t 


thermal  conductivity  of  dry  matter, 

J/( s  .  m  .  K) 

thermal  conductivity  of  frozen  food, 

J/( s  .  m  .  K) 

thermal  conductivity  of  ice, 

J/( s  .  m  .  K) 

thermal  conductivity  of  unfrozen  food, 

J/(s  .  m  .  K) 

thermal  conductivity  of  water, 

J/( s  .  m  .  K) 

thickness  of  infinite  slab,  or  diameter 
of  sphere  or  cylinder,  cm 

least  significant  range 

mass  of  copper  cylinder,  g 

mass  of  sample,  g 

number  of  time  increments  at  which 
experimental  and  predicted  temperatures  are 
compared 

number  of  elements 
Rh/k 

rate  of  heat  flow  from  block  to  sample, 
J/min 

radial  position,  cm 

number  of  replicates 

radius  of  cylinder,  cm 

standard  error  of  the  mean 

sum  of  all  the  measured  sucrose  contents 

sum  of  the  7  treatments  which  had  a 
replicate  10  measured 

sum  of  the  11  measured  replicates  of 
treatment  8 

significant  studentized  range 
time,  min 


I- 1  Spl 
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V 

trt 

T 


To 


T 

a 

Tb 

Tci 

T  c  a 

T  (t ) 

exp 

T/ 

TN+  1 

T  (t ) 
pre 

T 


U 


u 


c  i 


U  C  2 


U  c  3 


Wi 

W0 


time  for  centre  of  sample  to  reach  the 
freezing  point,  min 

number  of  treatments 

temperature,  K 

initial  temperature,  K 

temperature  of  cooling  or  heating  medium,  K 

temperature  of  block,  K 

initial  freezing  point,  K 

temperature  at  which  In  H  vs.  T  becomes 
linear  (during  freezing) 

temperature  at  which  k  vs.  T  is  no  longer 
linear  (during  freezing) 

average  experimental  temperature  at 
time  t,  K 

unique  freezing  point,  K 
surface  temperature,  K 
predicted  temperature  at  time  t,  K 
freezing  point  of  water,  K 

T  -  Ta 

temperature,  - 

T 0  *  Ta 


of  juice  which  has  become  ice 
of  soluble  solids  in  unfrozen 

juice 


dimensionless 


T0  -  Ta 


mass  fraction 
mass  fraction 


» 


s 
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w  mass  fraction  of  soluble  solids  in  residual 

s 

solut i on 

w  mass  fraction  of  insoluble  solids 

u 

w  mass  fraction  of  water 

w 

x  dimensionless  position,  r/R 

y  distance  in  one-dimensional  heat  equation,  cm 

Y  estimated  missing  value  in  Table  4,  %  sucrose 


Greek  Symbols 


a 

a 

w 

fi 

n 

AH  j(T) 
AR 
At 
AT 
AT  b 
Ax 
P 

T 


thermal  diffusity,  cm2/min 

thermal  diffusity  of  water,  cm2/min 

roots  of  fi  Ji(p  )  =  PJo (£  ) 

'n  *n  n 

Hj  (293.2)  -  Hj  (T)  ,  J/g 
element  thickness,  cm 
time  increment,  min 
Tb  -  T 

AT  without  sample 

dimensionless  element  thickness,  AR/R 
density,  g/m3 
dimensionless  time,  at/R2 


Subscripts  and  Superscripts 


d 


dry  matter 


\ 


frozen 


element  number 
time  increment 
juice 
unfrozen 
sample 
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17.  APPENDIX  A 


17.1  The  Suqarbeet 

The  sugarbeet  is  a  biennial  plant  which  produces  a 
large  tap  root  in  the  first  year  and  seed  in  the  following 
year.  Sugar  is  produced  by  photosynthesis,  in  the  leaves  of 
the  plant,  and  is  used  for  the  plant's  metabolism  and 
growth.  Excess  sugar  is  stored  in  the  tap  root.  The  top  part 
of  the  root,  from  where  the  leaves  grow,  is  called  the 
"crown".  The  sugar  content  varies  throughout  the  root.  It's 
minimum  value  is  attained  at  the  crown  where  fewer  storage 
cells  exist.  The  maximum  sugar  content  occurs  about  halfway 
down  the  root. 

The  sugarbeet  is  composed  of  white  storage  tissue  which 
is  penetrated  by  darker  vascular  rings.  The  vascular  rings 
are  composed  of  xylem  and  phloem  vessels.  The  xylem  conducts 
water  and  the  phloem  conducts  nutrients.  For  example,  sugar 
is  transported  from  the  leaves  to  the  root  by  means  of  the 
phloem.  The  storage  tissue  consists  of  thin  walled 
"parenchyma"  cells.  Small  parenchyma  cells  surround  the 
vessels  whereas  larger  parenchyma  cells  are  further  away. 

The  large  cells  are  primarily  responsible  for  water  storage 
and  have  a  much  lower  sugar  content  than  the  small  cells. 

The  "intercellular  space"  between  cells  occupies  about  0.5 
to  2%  of  the  total  volume  (67). 

It  is  important  to  remember  that  the  sugarbeet  root  is 
composed  of  different  types  of  cells  of  varying  composition. 
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Each  of  these  cells  may  freeze  at  a  different  temperature. 
Oldfield  et  al.  (157)  observed  that  sugarbeets  start  to 
freeze  at  271  R.  However,  Wyse  (230)  stated  that  the 
intercellular  region  and  some  cells  (epidermal  cells)  near 
the  surface  of  the  beets  may  freeze  at  272  K.  Different 
cells  may  also  differ  in  their  tolerance  to  freezing  when  it 
does  occur.  Phloem  cells  are  more  resistant  to  frost  injury 
than  parenchyma  cells  (214). 


18.  APPENDIX  B 


18.1  Freezing  Injury  in  Plants 

Plants  are  composed  of  small  units  called  "cells",  each 
of  which  is  bounded  by  a  "cell  wall".  Every  cell  is  divided 
into  a  number  of  different  compartments,  each  of  which  has 
its  own  specific  functions.  Separating  the  various 
compartments  is  the  "ground  cytoplasm"  which  is  a  complex 
colloidal  solution.  Changes  in  the  concentration  or  pH  of 
this  solution  can  irreversibly  damage  the  cell  (106). 

Biological  materials  are  injured  by  freezing  and 
thawing.  This  injury  results  from  the  conversion  of  water 
into  ice.  When  the  temperature  falls  below  the  initial 
freezing  point,  ice  begins  to  form  in  the  intercellular 
region.  Water  then  leaves  the  inside  of  the  cell  in  an 
attempt  to  equalize  the  chemical  potentials  inside  and 
outside  the  cell.  Thus,  a  growing  ice  layer  is  formed 
outside  the  cell.  At  the  same  time  the  solute  concentration 
inside  the  cell  is  increased  and  this  lowers  the  freezing 
point  of  the  solution  inside  the  cell  ("intracellular 
solution").  If  the  temperature  is  decreased  sufficiently 
slowly  so  water  can  leave  the  cell  then  intracellular  ice 
formation  will  be  avoided.  However,  if  the  sample 
temperature  falls  very  rapidly  then  water  cannot  leave  the 
cell  quickly  enough  and  ice  is  eventually  formed  inside  the 
cell.  This  intracellular  ice  is  considered  to  be  lethal  for 
cells.  For  frozen  foods  lethality  may  not  be  important  as 
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cells  need  not  survive  freezing  and  thawing  if  the  food  is 
consumed  within  a  short  time  of  thawing. 

The  increased  solute  concentration,  obtained  during 
freezing,  may  result  in  precipitation  of  some  solutes  (65), 
attack  on  cell  membranes  (5,6,87,125,183,191),  and 
acceleration  of  some  chemical  reactions  (119).  Damage  to 
cell  membranes  can  result  in  the  movement  of  enzymes  to 
regions  of  the  tissue  where  they  can  act  on  reactants  to 
which  they  do  not  normally  have  access.  Loss  of  water  leads 
to  cell  shrinkage,  which,  in  turn,  brings  potential 
reactants  closer  together.  This  can  result  in  protein 
denaturation  (119,120)  and  degradation  of  lipid  components 
(6,219,236)  of  cell  membranes. 

The  rate  constant  for  most  chemical  reactions  decreases 
with  decreasing  temperature.  However,  one  must  consider  that 
the  cell  solution  is  concentrated  when  water  is  transformed 
into  ice.  Thus,  the  rate  of  some  chemical  reactions  is 
actually  increased  during  freezing.  These  reactions  are 
subsequently  slowed  down  when  the  temperature  falls  below 
about  263  K.  This  is  because:  (1)  the  reaction  rate  constant 
continues  to  fall  with  temperature,  and  (2)  the  viscosity  of 
the  cell  solution  is  increasing  with  decreasing  temperature 
and  will  eventually  become  so  great  as  to  essentially  stop 
the  movement  of  potential  reactants. 

When  water  freezes  it  expands  about  9%.  This  expansion 
may  cause  the  growing  ice  layer  to  push  adjoining 
apart  or  to  even  penetrate  some  of  the  cells 
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(57,77,119,145,173,198) . 

During  thawing  intercellular  ice  melts  and  flows  back 
into  the  cell  and  this  causes  a  stress  on  the  cell  wall  and 
membranes  (117).  In  addition,  thawing  of  ice  momentarily 
results  in  large  concentration  gradients  inside  the  cell. 
Regions  of  almost  pure  water  are  adjacent  to  regions  of 
concentrated  solutes.  The  resulting  high  osmotic  gradient 
can  damage  the  cell. 

The  plasma  membrane  is  situated  just  inside  the  cell 
wall  and  it  maintains  the  integrity  of  the  cell.  It  allows 
some  molecules  to  enter  the  cell  but  stops  others.  If  the 
plasma  membrane  is  damaged  then  it  may  lose  this 
semipermeability  and  solutes  are  lost  from  the  interior  of 
the  cell.  This  results  (because  of  a  reduction  in  osmotic 
pressure)  in  a  decrease  in  the  amount  of  water  which 
re-enters  the  cell  during  thawing.  Loss  of  water  may  also  be 
partly  attributed  to  the  reduction  in  water  holding  capacity 
of  proteins  denatured  during  freezing.  A  decrease  in  cell 
water  content  causes  a  reduction  in  the  pressure  on  the  cell 
wall  and  results  in  a  softening  of  the  tissue.  Cell  walls, 
because  of  their  rigid  nature,  are  partly  responsible  for 
the  firmness  of  a  plant  tissue.  Thus,  damage  to  cell  walls 
and  membranes  during  freezing  contributes  to  the  softening 
of  the  tissue. 

As  stated  above,  cell  membranes  can  be  irreversibly 
damaged  during  freezing  and  lose  their  semipermeability. 
There  are  many  possible  causes  of  membrane  damage. 
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Lovelock  (125)  stated  that  concentration  of  electrolytes, 
rather  than  the  mechanical  effect  of  ice  crystals,  was  the 
main  cause  of  damage  in  the  freezing  of  red  blood  cells. 

When  the  cells  were  placed  in  NaCl  solutions  stronger  than 
0.8  M  the  plasma  membrane  became  permeable  to  sodium  ions. 

Levitt  (119)  considered  that  when  the  cell  contracts 
during  slow  freezing  a  "break”  occurs  in  the  lipid  layer  of 
its  membranes.  This  break  is  then  stabilized  by  the 
formation  of  disulfide  linkages  between  membrane  proteins 
and  is  not  removed  during  thawing.  For  intracellular 
freezing,  Levitt  (122)  speculated  that  submic roscopic 
crystals  "expand  into  the  plasma  membrane,  creating  small 
holes  in  the  now  solid  lipid  layer".  During  thawing,  these 
holes  cause  a  loss  of  semipermeability  in  the  plasma 
membrane . 

Meryman  et  al.  (143)  postulated  that  irreversible 
damage  to  cell  membranes  occurs  during  slow  freezing  when 
the  cell  volume  is  reduced  below  a  certain  critical  value 
(about  "40  to  50%  of  the  unfrozen  volume").  This  damage  was 
attributed  to  loss  of  lipid  or  protein  components  from  the 
cell  membrane.  The  loss  of  the  membrane  material  was 
considered  to  be  caused  by  either  a  direct  osmotic  stress  or 
physical  forces  resulting  from  the  reduction  in  cell  volume. 

Araki's  (6)  experimental  results,  for  mitochondrial 
membranes,  indicated  that  phospholipids  (particularly 
phosphat idy lethanolamine )  may  be  degraded  during  slow 
freezing  by  the  action  of  a  phospholipase.  This  enzyme  is 
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apparently  activated  by  changes  inflicted  on  the  membrane 
following  exposure  to  high  solute  concentrations. 

18.2  Effect  of  Freezing  and  Thawing  Rates 

Fast  freezing  generally  produces  food  of  a  better 
quality  than  slow  freezing  (28,77,180,220).  Fast  freezing 
results  in  small  distributed  ice  crystals.  Slow  freezing 
produces  large  ice  masses  which  tend  to  physically  disrupt 
the  cell  walls.  Another  disadvantage  of  slow  freezing  is 
that  cells  lose  a  large  amount  of  water  during  freezing.  The 
remaining  solution  then  becomes  highly  concentrated.  The 
probability  of  unfavourable  reactions  or  solute 
precipitation  occurring  is  therefore  increased. 

Chemical  reactions  can  take  place  quite  rapidly  above 
263  K.  The  faster  the  freezing  rate  the  quicker  the  food 
temperature  will  go  below  263  K  and  the  less  quality 
deterioration  will  occur.  However,  it  should  be  noted  that 
if  large  samples  are  frozen  too  quickly,  for  example  in 
liquid  nitrogen,  the  surface  will  be  frozen  solid  before  the 
interior  starts  to  freeze.  When  the  interior  expands  during 
freezing  severe  cracking  occurs. 

Slow  thawing  allows  growth  of  ice  crystals 
(recrystallization)  and  thus  minimizes  the  advantages  of 
fast  freezing.  It  also  permits  more  time  for  undesirable 
reactions  to  occur.  However,  it  should  be  noted  that  fast 
thawing  has  the  following  disadvantages  : 
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1.  Ice  is  melted  rapidly  and  neighbouring  regions  of  almost 
pure  water  and  concentrated  solution  will  exist. 

2.  Water  is  returned  very  quickly  through  the  cell  wall  and 
membrane,  thus  creating  a  shear  stress  on  these 
structures . 

The  disadvantages  of  fast  thawing  are  greater  following  slow 
freezing  as  more  water  must  be  returned  to  the  cell.  Thus, 
as  a  general  rule,  fast  freezing  should  be  followed  by  fast 
thawing  and  slow  freezing  by  slow  thawing. 

Rates  of  thawing  are  limited  by  the  fact  that  the 
surface  temperature  cannot  be  allowed  to  reach  too  high  a 
value.  If  it  is  desired  to  keep  a  plant  alive  after  thawing 
then  the  surface  must  be  kept  below  318  K.  When  dealing  with 
a  food,  higher  surface  temperatures  may  be  used  but  care 
must  still  be  taken  to  avoid  damaging  the  surface.  Thawing 
rates  are  also  limited  by  the  fact  that  heat  must  be 
transferred  through  a  growing  layer  of  thawed  food  which  has 
a  thermal  diffusity  about  nine  times  less  than  that  of  the 
frozen  food. 


18.3  Storage  Temperature 

Frozen  vegetables  are  usually  stored  at  255  K  as 
biochemical  reactions  take  place  very  slowly  at  this 
temperature.  Frozen  storage  above  263  K  is  not  recommended 
as  above  this  temperature  : 
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1.  Micro-organisms  can  grow. 

2.  Some  biochemical  reaction  rates  are  high. 

■% 

The  actual  storage  temperature  is  decided  by  economics  i.e. 
the  potential  improvement  in  vegetable  quality  or  storage 
life  must  be  balanced  against  the  increased  refrigeration 
cost  of  a  lower  storage  temperature. 

As  a  result  of  ice  crystal  growth  the  benefits  of  fast 
freezing  may  be  lost  during  storage  (180,224).  The  lower  the 
storage  temperature  the  slower  will  be  the  ice  growth  rate. 
It  is  generally  accepted  (103,180,208)  that  fluctating 
storage  temperatures  encourage  the  growth  of  ice  crystals 
and  thus  the  storage  temperature  should  be  kept  as  constant 
as  possible.  However,  some  experiments  (92)  have  indicated 
that  fluctuations  in  the  storage  temperature  are  not 
important  for  food  storage. 

Enzymatic  reactions  can  continue  at  frozen  storage 
temperatures.  Thus,  vegetables  usually  undergo  a  blanching 
treatment  prior  to  freezing.  It  consists  of  placing  the 
vegetable  in  steam  or  hot  water  (353  to  373  K)  for 
approximately  one  minute.  Blanching  deactivates  some  enzymes 
and  kills  many  micro-organisms  (159).  It  also  removes  air 
from  vegetable  tissues  and  thus  reduces  the  probability  of 
oxidation  occurring  (91). 


19.  APPENDIX  C 


19.1  Estimation  of  Sucrose  Content  of  Beet  Cylinders 

The  following  method  was  used  to  determine  the  sucrose 

content  of  beet  cylinders.  It  was  adapted  from  a  procedure 

developed  by  Karr  and  Norman  (109a). 

1.  Weigh  a  closed  container  to  4  decimal  places  (=A). 

2.  Add  3  beet  cylinders  and  reweigh  to  4  decimal  places 
( =B )  . 

3.  Transfer  the  3  cylinders  to  a  4  oz .  oster  blender  jar. 

4.  Cut  the  cylinders  into  small  pieces  and  add  100  ml  of 
distilled  water. 

5.  Add  0.4  g  of  dry  basic  lead  subacetate  (aids  in 
subsequent  filtration). 

6.  Blend  for  6  min  at  the  top  speed  of  a  CYCLOMATIC  14 
OSTERIZER. 

7.  Filter  the  solution  obtained  into  a  preweighed  (to  1 
decimal  place)  plastic  container  (weight  of  container, 
including  top,  =C). 

8.  In  this  procedure  it  is  assumed  that  no  sucrose  is  lost 
during  extraction.  Consequently,  after  pouring  the 
extract  solution  into  the  funnel  the  remaining  extract, 
in  the  blender  jar  and  on  the  cutter  blades  etc.,  is 
washed  into  the  filtration  funnel  with  distilled  water. 
The  precipitate  left  on  the  filter  paper  (WHATMAN  no.  4, 
24  cm  diameter)  after  filtration  is  also  washed  with 
distilled  water.  The  total  amount  of  distilled  water 
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used  for  washing  is  about  170  ml. 

9.  The  top  is  then  put  on  the  container  and  the  container 
and  extract  solution  are  weighed  to  1  decimal  place 

( =D )  . 

10.  Weigh  a  5  ml  volumetric  flask  (=E). 

11.  Add  about  0.9  ml  of  the  extract  solution  and  reweigh 
( =F )  . 

12.  Add  about  0.9  ml  of  a  trehalose  solution  of  known 
concentration  (should  be  about  0.007  g  trehalose/g  of 
solution)  and  reweigh  (=G). 

The  above  three  weights  are  all  recorded  to  4  decimal 
places . 

13.  The  solution  in  the  volumetric  flask  is  mixed  well  and 
then  0.8  ml  of  the  solution  is  added  to  a  labelled, 
empty,  6  ml  hypovial. 

14.  The  hypovial  is  placed  in  a  vacuum  oven  (at  333  K  and 
10.2  kPa  (3  in.  Hg))  for  1  hour  to  remove  all  the  water 
from  the  solution. 

15.  The  hypovial  is  sealed  and  0.5  ml  of 

tr imethyls i lyl imidazole  and  0.5  ml  of  dimethylf ormamide 
are  injected.  The  resulting  solution  is  left  overnight 
at  room  temperature  or  for  15  minutes  at  333  K. 

16.  Then  9.8  microlitres  of  this  solution  are  injected  into 
a  Gas  Chromatograph  (F  and  M  model  810). 

Using  the  resulting  ratio  of  the  sucrose  and  trehalose 
peak  areas,  the  %  sucrose  can  be  calculated  from  the 
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formula : 


%  SUCROSE 
IN  COSSETTES 


AREA  CONC.TREH.  (G-F)(D-C) 

X  X  K  X  -  X  100% 

RATIO  SOLUTION  (F-E)(B-A) 

( C  —  1 ) 


The  constant  K  is  the  calibration  constant  for  the  Gas 
Chromatograph  and  is  defined  by: 

trehalose  peak  area/trehalose  mass  injected 

K  =  -  (C-2) 

sucrose  peak  area/sucrose  mass  injected 

It  is  obtained  by  carrying  out  calibration  runs  using 
standard  sucrose  and  trehalose  solutions  of  known 
concentration.  These  standard  solutions  were  injected,  every 
day  samples  were  being  analysed,  to  detect  any  variations  in 

K. 


Gas  Chromatograph  Details 

A  single  column  flame  ionization  detector  was  used  to  detect 
the  sugar  peaks.  The  response  from  the  detector  was 
transmitted  to  a  HEWLETT  PACKARD  computer  (21MX  E-series), 
via  an  A/D  converter,  where  the  peak  areas  were  determined 
and  reported  on  a  telytype. 
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G . C .  Operating  Conditions 


Helium  flow  rate:  30  ml/min. 

Hydrogen  flow  rate:  45  ml/min. 

Air  flow  rate:  440  ml/min. 

Helium  inlet  pressure:  380  kPa  (40  psig). 
Hydrogen  inlet  pressure:  250  kPa  (22  psig). 
Air  inlet  pressure:  380  kPa  (40  psig). 

Oven  temperature:  533  K 
Injection  port  temperature:  573  K 
Flame  detector  temperature:  583  K 


20.  APPENDIX  D 


20.1  Est ima t i on  of  Missing  Value  in  Table  4 

Due  to  experimental  problems  the  percent  sucrose  for 
replicate  10  of  treatment  8  in  the  first  set  of  experiments 
was  not  measured.  As  a  value  was  required  for  the  Randomized 
Block  Design  it  was  decided  to  estimate  one  using  the 
formula  in  Snedecor  and  Cochran  (195): 

trt  X  St8  +  rep  X  SrlO  -  Sail 

Y  = - (D-l) 

( t rt  -  1 )  ( rep  -  1 ) 


where : 

Y  =  estimated  ^sucrose  of  treatment  8,  replicate  10. 
trt  =  number  of  treatments  =  8 
rep  =  number  of  replicates  =  12 

St8  =  sum  of  the  11  measured  replicates  of  treatment  8 
=  193.4 

SrlO  =  sum  of  the  7  treatments  which  had  a  replicate  10 
measured  =  110.8 

Sail  =  sum  of  all  the  measured  sucrose  contents  =  1624.8 
Thus , 


Y 


8  X  193.4  +  12  X  110.8  -  1624.8 
7  X  11 


(D-2) 


=  16.3% 
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20.2  Duncan 1 s  Multiple  Range  Test 

20.2.1  First  Set  of  Storage  Experiments 


Randomized  Block  Design 
a)  Rank  Treatment  (Trt.)  Means 


Trt. 2  Trt. 4  Trt. 3 

Trt . 5  Trt . 6 

Trt  .8 

Trt .  7 

Trt .  1 

16.4  16.6  16.8 

16.9  17.5 

17.5 

17.6 

17.6 

Calculate  the  Standard  Error  of 

the  Mean 

(sm) 

From  Table  4  the  Error  Mean  Square  (EMS)  =  0.1849 
Number  of  replicates  (rep)  =  12 

s  =  /  EMS/rep  =  Jo. 1849/12  =  0.1241  (D-3) 

m 

c)  Degrees  of  Freedom  (DF)  for  Error 
DF  =  (trt  -  1) (rep-1) 

=  (8-1)  (12-1)  =  77  (D-4) 

d)  5%  Level  of  Significance 

For  DF=77  a  linear  interpolation  of  Harter's  (94)  Table  1 
yields : 

No.  of  means  2345678 
SSR  2.82  2.97  3.06  3.13  3.19  3.23  3.27 

where : 

SSR  =  significant  studentized  range. 
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e)  Calculate  the  Least  Significant  Ranges  (LSR) 

LSR  =  sm  X  SSR  (D-5) 

This  yields  the  following  LSR  values: 

No.  of  Means  2345678 
LSR  .35  .37  .38  .39  .40  .40  .41 

f)  Determine  the  Treatments  which  are  Significantly 
Di f  f  erent 

Trt.l  -  LSR e  =  17.6  -  0.41  =  17.19 
Trt.5  =  16.9 

Therefore,  Trt.l  is  significantly  greater  than  treatments 
2-5. 

Trt.l  -  LSR4  =  17.6  -  0.38  =  17.22 
Trt .6  =  17.5 

Therefore,  Trt.l  is  not  significantly  different  from 
treatments  6-8. 

Trt. 8  -  LSR t  =  17.5  -  0.40  =  17.10 
Trt.5  =  16.9 

Therefore,  Trt. 8  is  significantly  different  from 
treatments  2-5. 

Trt.5  -  LSR4  =  16.9  -  0.38  =  16.52 
Trt. 2  =  16.4 

Therefore,  Trt.5  is  significantly  different  from  Trt. 2. 

Trt.5  -  LSR3  =  16.9  -  0.37  =  16.53 
Trt. 4  =  16.6 

Therefore,  Trt.5  is  not  significantly  different  from 
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treatments  4  and  3. 

Trt.3  -  LSR3  =  16.8  -  0.37  =  16.43 
Trt.2  =  16.4 

Therefore,  Trt.3  is  significantly  different  from  Trt.2. 

Trt.4  -  LSR2  =  16.6  -  0.35  =  16.25 
Therefore,  Trt.4  is  not  significantly  different  from  Trt.2. 
The  above  results  can  be  summarized  as  follows: 

Trt.2  Trt.4  Trt.3  Trt.5  Trt.6  Trt.8  Trt.7  Trt.l 


Two  treatments  which  are  underlined  with  the  same  line  are 
not  significantly  different  at  the  5%  level.  If  two 
treatments  are  not  underlined  with  the  same  line  then  they 
are  significantly  different.  The  calculations  were  repeated 
at  the  1%  level  and  the  following  results  were  obtained: 

Trt.2  Trt.4  Trt.3  Trt.5  Trt.6  Trt.8  Trt.7  Trt.l 


Factorial  Design 

The  six  treatments  in  order  of  ranking  are: 
Trt.2  Trt.4  Trt.3  Trt.5  Trt.6  Trt.7 

16.4  16.6  16.8  16.9  17.5  17.6 

From  Table  7: 

EMS  =  0.210 
rep  =  12 


s 

m 


=  j/  EMS/rep  =  /0. 210/12  =  0.1323 


(D-6) 
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DF  for  Error  =  (trt  -  l)(r-l) 


(6-1) (12-1)  =  55 


(D-7) 


For  a  5%  level  of  significance 
use  of  Harter's  Table  1  gives: 


No.  of  Means  2 


3 


4 


5 


6 


SSR 


2.84  2.98  3.08  3.15  3.20 


LSR 


0.38  0.39  0.41  0.42  0.42 


Trt. 7  -  LSR6  =  17.6  -  0.42  =  17.18 
Therefore,  Trt. 7  is  significantly  greater  than 
treatments  2-5. 

Trt. 7  -  LSR2  =  17.6  -  0.38  =  17.22 
Therefore,  Trt. 7  is  not  significantly  different  from  Trt. 6 
Trt. 6  -  LSR s  =  17.5  -  0.42  =  17.08 
Therefore,  Trt. 6  is  significantly  greater  than 
treatments  2-5. 

Trt. 5  -  LSR4  =  16.9  -  0.41  =  16.49 
Therefore,  Trt. 5  is  significantly  different  from  Trt. 2. 

Trt. 5  -  LSR3  =  16.9  -  0.39  =  16.51 
Therefore,  Trt. 5  is  not  significantly  different  from 
treatments  4  and  3 

Trt. 3  -  LSR3  =  16.8  -  0.39  =  16.41 
Therefore,  Trt. 3  is  significantly  different  from  Trt. 2 
Trt. 4  -  LSR2  =  16.6  -  0.38  =  16.22 
Therefore,  Trt. 4  is  not  significantly  different  from  Trt. 2. 
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The  above  results  are  now  summarized: 

Trt.2  Trt.4  Trt.3  Trt.5  Trt.6  Trt.7 


The  following  result  was  obtained  for  a  1%  level  of 
signi f icance : 

Trt.2  Trt.4  Trt.3  Trt.5  Trt.6  Trt.7 


20.2.2  Second  Set  of  Storage  Experiments 

The  four  treatments  in  order  of  ranking  are: 

Trt .12  Trt.10  Trt.l  Trt.ll 

17.1  17.5  17.5  17.6 

From  Table  13: 


EMS  =  0.084 

rep  =  8 

sm=  /EMS/rep  = 

/0 . 084/8  =  0.1025 

(D-8) 

DF  for  Error  = 

21 

(D-9) 

For  a  5%  level  of  significance,  use  of  Harter's 
Table  1  yields: 

No.  of  means:  2 

SSR  2.94 


3  4 

3.09  3.18 


LSR 


0.30 


0.32 


0.33 


. 
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Trt.ll  -  LSR4  =  17.6  -  0.33  =  17.27 
Therefore,  Trt.ll  is  significantly  greater  than  Trt.12. 

Trt.ll  -  LSR3  =  17.6  -  0.32  =  17.28 
Therefore,  Trt.ll  is  not  significantly  greater  than 
treatments  1  and  10. 

Trt.l  -  LSR3  =  17.5  -  0.32  =  17.18 
Therefore,  Trt.l  is  significantly  greater  than  Trt.12. 

Trt.10  -LSR2  =  17.5  -  0.30  =  17.20 
Therefore,  Trt.10  is  significantly  greater  than  Trt.12. 
The  above  results  are  summarized  as  follows: 

Trt.12  Trt.10  Trt.l  Trt.ll 


The  following  result  was  obtained  for  a  1%  level  of 
signi f icance : 


Trt.12  Trt.10  Trt.l  Trt.ll 


A 


21.  APPENDIX  E 


21 . 1  Experimental  Results  Used  in  Determining  h 


TABLE  E-l:  AIR  HEATING  OF  COPPER  CYLINDER,  T0=243  K 


TIME 

TEMPERATURE 

(min ) 

(K) 

RUN  1 

RUN  2 

RUN  3 

(T  =297.7) 

a 

(T  =297.5) 

a 

(T  =297.5) 

a 

0 

243.1 

242.9 

242.9 

1 

250.5 

249.9 

250.5 

2 

257.7 

256.5 

256.7 

4 

268.9 

266.5 

266.7 

8 

281.2 

278.3 

278.0 

12 

286.2 

283.5 

283.2 

16 

288.2 

286.2 

285.5 

20 

289.3 

288.2 

287.6 

24 

290.5 

289.8 

288.9 

28 

291.3 

290.9 

289.9 

32 

292.1 

292.1 

290.8 

36 

292.7 

292.9 

291.7 

40 

293.3 

293.7 

292.3 

44 

293.9 

294.5 

292.7 

48 

294.2 

294.7 

293.3 

52 

294.7 

295.2 

293.9 

216 


■ 


- 
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RUN  1 

RUN  2 

RUN  3 

56 

295.1 

295.6 

294.3 

60 

■» 

295.2 

315.9 

294.6 

64 

295.3 

316.0 

294.8 

' 
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TABLE  E-2:  AIR  HEATING  OF  COPPER  CYLINDER,  T0=253  K 


TIME 

TEMPERATURE 

(min) 

(K) 

RUN  1 

RUN  2 

RUN  3 

(T  =297.8) 

(T  =297.9) 

(T  =297.9) 

a 

a 

a 

0 

252.9 

252.9 

252.9 

1 

258.3 

258.2 

258.4 

2 

263.7 

263.7 

263.7 

4 

271.7 

271.4 

271.4 

8 

281.2 

280.5 

280.3 

12 

284.9 

284.3 

284.1 

16 

287.3 

286.5 

286.4 

20 

288.9 

288.3 

288.2 

24 

289.9 

289.8 

289.3 

28 

290.9 

290.9 

290.2 

32 

291.9 

291.9 

290.9 

36 

292.7 

292.7 

291.9 

40 

293.6 

293.3 

292.6 

44 

293.9 

293.8 

293.2 

48 

294.3 

294.3 

293.7 

52 

294.6 

294.5 

294.2 

56 

294.7 

294.7 

294.5 

60 

294.9 

294.8 

294.7 

64 

295.2 

295.1 

294.9 

68 

295.6 

295.5 

295.1 

72 

— 

295.6 

— 

' 

* 
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TABLE  E-3:  AIR  HEATING  OF  COPPER  CYLINDER,  T0=263  K 


TIME 

TEMPERATURE 

(min ) 

(K) 

RUN  1 

RUN  2 

RUN  3 

(T  =297.9) 

3 

(T  =297.8) 

3 

(T  =297.9) 

0 

263.2 

263.2 

263.2 

1 

267.3 

267.3 

267.2 

2 

270.9 

271.2 

271.1 

4 

276.8 

276.8 

276.7 

8 

283.4 

283.3 

282.9 

12 

286.9 

286.7 

286.3 

16 

289.3 

288.8 

288.6 

20 

290.5 

290.3 

290.2 

24 

291.7 

291.6 

291.4 

28 

292.6 

292.6 

292.3 

32 

293.1 

293.3 

293.3 

36 

293.9 

293.9 

293.8 

40 

294.7 

294.5 

294.4 

44 

295.1 

294.8 

294.7 

48 

295.5 

295.1 

295.2 

' 
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TABLE  E-4:  COOLING  OF  COPPER  CYLINDER ,  T  =243  K 

'  a 


TIME 

TEMPERATURE 

(min ) 

(K) 

RUN  1 

RUN  2 

RUN  3 

0 

297.6 

297.2 

296.1 

0.167 

284.7 

281.9 

276.6 

LO 

• 

o 

266.7 

263.1 

259.8 

1.0 

253.5 

251.3 

250.1 

1.5 

247.9 

246.7 

246.1 

2.0 

245.4 

244.5 

244.2 

2.5 

244.1 

243.4 

243.3 

3.0 

243.5 

243.3 

- 

3.5 

243.1 

- 

- 

1  ^  ■ 
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TABLE  E-5:  COOLING  OF  COPPER  CYLINDER,  T  =253  K 

a 


TIME 

TEMPERATURE 

(min ) 

(K) 

RUN  1 

RUN  2 

RUN  3 

0 

297.3 

296.8 

296.6 

0.167 

286.2 

281.3 

278.8 

0.5 

269.3 

267.5 

264.9 

1.0 

259.1 

258.4 

257.2 

1.5 

255.5 

255.2 

254.7 

2.0 

254.1 

253.8 

253.7 

2.5 

253.2 

253.2 

- 

3.0 

253.1 

- 

- 
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TABLE  E-6:  COOLING  OF  COPPER  CYLINDER,  Ta=263  K 


TIME 

TEMPERATURE 

(min ) 

(K) 

RUN  1 

RUN  2 

RUN  3 

0 

297.9 

298.6 

298.5 

0.167 

282.3 

284.3 

282.2 

o 

• 

<J1 

270.9 

272.7 

271.3 

1.0 

265.4 

266.2 

265.6 

1.5 

264.0 

264.1 

264.1 

2.0 

263.5 

263.7 

263.5 

2.5 

- 

263.3 

- 

TABLE  E-7:  AIR  HEATING  -  h  AS  A  FUNCTION  OF  u 


N  +  i 

(To  =  243  K) 


TIME 

RUN  1 

RUN  2 

RUN 

3 

(min ) 

u  N  +  i 

h 

UN  +  1 

h 

u  N  +  1 

h 

0 

1.0 

- 

1.0 

- 

1.0 

- 

1 

0.864 

20.2 

0.872 

18.8 

0.861 

20.2 

2 

0.733 

23.0 

0.751 

20.9 

0.747 

19.5 

4 

0.527 

23.0 

0.568 

19.5 

0.564 

19.5 

8 

0.302 

19.5 

0.352 

16.7 

0.357 

15.3 

12 

0.211 

12.6 

0.256 

11.2 

0.262 

10.5 

16 

0.174 

6.3 

0.207 

7.7 

0.220 

6.3 

20 

0.154 

4.2 

0.170 

7.0 

0.181 

7.0 

24 

0.132 

5.6 

0.141 

6.3 

0.158 

4.9 

28 

0.117 

4.2 

0.121 

5.6 

0.139 

4.2 

32 

0.103 

4.9 

0.099 

7.0 

0.123 

4.2 

36 

0.092 

4.2 

0.084 

5.6 

0.106 

4.9 

40 

0.081 

4.2 

0.070 

6.3 

0.095 

3.5 

44 

0.070 

4.9 

0.055 

8.4 

0.088 

2.8 

48 

0*. 064 

2.8 

0.051 

2.1 

0.077 

4.9 

52 

0.055 

5.6 

0.042 

7.0 

0.066 

5.6 

56 

0.048 

4.9 

0.035 

6.3 

0.059 

4.2 

60 

0.046 

1.4 

0.029 

6.3 

0.053 

3.5 

64 

0.044 

1.4 

0.027 

2.1 

0.049 

2.8 

h  is 

expressed 

in  units  of  J/(s 

.  m2 

.  K) 

. 

' 
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TABLE  E-8:  AIR  HEATING  -  h  AS  A  FUNCTION  OF  uN+1 


(T0 

=  253  K) 

TIME 

RUN 

1 

RUN 

2 

RUN 

3 

(min ) 

h 

u  N  +  1 

h 

UN  4  1 

h 

0 

1.0 

— 

1.0 

- 

1.0 

- 

1 

0.880 

17.4 

0.882 

16.7 

0.878 

18.1 

2 

0.759 

20.2 

0.760 

20.9 

0.760 

20.2 

4 

0.581 

18.8 

0.589 

17.4 

0.589 

17.4 

8 

0.369 

15.3 

0.387 

14.7 

0.391 

14.0 

12 

0.287 

9.1 

0.302 

8.4 

0.307 

8.4 

16 

0.234 

7.0 

0.253 

6.3 

0.256 

6.3 

20 

0.198 

5.6 

0.213 

6.3 

0.216 

5.6 

24 

0.176  . 

4.2 

0.180 

5.6 

0.191 

4.2 

28 

0.154 

4.9 

0.156 

4.9 

0.171 

3.5 

32 

0.131 

5.6 

0.133 

5.6 

0.156 

3.5 

36 

0.114 

4.9 

0.116 

4.9 

0.133 

5.6 

40 

0.094 

7.0 

0.102 

4.2 

0.118 

4.2 

44 

0.087 

2.8 

0.091 

4.2 

0.104 

4.2 

48 

0.078 

3.5 

0.080 

4.2 

0.093 

4.2 

52 

0.071 

2.8 

0.076 

2.1 

0.082 

4.2 

56 

0.069 

1.4 

0.071 

2.1 

0.076 

2.8 

60 

0.065 

2.1 

0.069 

1.4 

0.071 

2.1 

64 

0.058 

3.5 

0.062 

3.5 

0.067 

2.1 

68 

0.049 

5.6 

0.053 

5.6 

0.062 

2.1 

72 

— 

— 

0.051 

1.4 

- 

— 

• 

- 

TABLE  E-9:  AIR  HEATING  -  h  AS  A  FUNCTION  OF  uN + x 


(Te  =  263  K) 


TIME 

RUN 

1 

RUN 

2 

RUN 

3 

(min ) 

UN  +  i 

h 

UN  +  i 

h 

UN  +  i 

h 

0 

1.0 

- 

1.0 

- 

1.0 

- 

1 

0.882 

17.4 

0.882 

17.4 

0.882 

17. 

2 

0.778 

17.4 

0.769 

18.8 

0.770 

18. 

4 

0.608 

17.4 

0.607 

16.0 

0.609 

16. 

8 

0.418 

13.3 

0.419 

12.6 

0.431 

11. 

12 

0.317 

9.8 

0.321 

7.0 

0.333 

9. 

16 

0.248 

8.4 

0.260 

7.0 

0.267 

7. 

20 

0.213 

5.6 

0.217 

6.3 

0.221 

6. 

24 

0.179 

6.3 

0.179 

6.3 

0.187 

5. 

28 

0.150 

5.6 

0.150 

6.3 

0.161 

4. 

32 

0.138 

3.5 

0.130 

4.9 

0.132 

7. 

36 

0.115 

6.3 

0.113 

4.9 

0.118 

4. 

40 

0.092 

7.7 

0.095 

5.6 

0.101 

5. 

44 

0.081 

4.9 

0.087 

3.5 

0.092 

2. 

48 

0.069 

5.6 

0.078 

3.5 

0.078 

5. 

h  is  expressed  in  units  of  J/(s  .  m2  .  K) 


4 

8 

0 

9 

1 

7 

3 

6 

9 

0 

2 

6 

8 

6 


• 
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TABLE 

;  E-10: 

LIQUID 

COOLING  - 

h  AS  A 

FUNCTION 

o 

c 

TIME 

RUN 

1 

( T  ^  =  243 

RUN 

K) 

2 

RUN 

3 

(min ) 

uN+i 

h 

UN  ♦  i 

h 

UN  +  i 

h 

““  ■“  “■ 

* “ 

““  "" 

“■ 

0 

1.0 

- 

o 

• 

1 — 1 

- 

1.0 

- 

0.167 

0.764 

223 

0.718 

275 

0.633 

380 

0.5 

0.435 

235 

0.372 

273 

0.318 

288 

1 

0.194 

225 

0.155 

243 

0.135 

237 

1.5 

0.091 

208 

0.070 

155 

0.060 

225 

2 

0.046 

193 

0.029 

105 

0.024 

250 

2.5 

0.022 

203 

0.009 

322 

0.008 

327 

3 

0.011 

193 

0.007 

62 

- 

- 

3.5 

0.004 

305 

— 

- 

- 

- 

h  is  expressed  in  units  of  J/(s  .  m2  .  K) 


' 
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TABLE  E - 1 1 :  LIQUID  COOLING  - 

h  AS  A  FUNCTION  OF  uN 

(Ta  =  253 

K) 

TIME 

RUN 

1 

RUN 

2 

RUN 

3 

(min ) 

•h  i 
+  i 

Z  1 

P  i 

h 

UN  +  1 

h 

h 

0 

o 

• 

1 — 1 

- 

1.0 

- 

1.0 

- 

0.165 

0.750 

240 

0.647 

362 

0.593 

435 

IT) 

• 

O 

0.369 

295 

0.333 

277 

0.275 

320 

1 

0.140 

270 

0.125 

270 

0.098 

285 

1.5 

0.059 

240 

0.052 

242 

0.041 

242 

2 

0.027 

215 

0.021 

260 

0.018 

225 

2.5 

0.007 

385 

0.007 

305 

- 

- 

3 

0.005 

112 

- 

— 

— 

— 

.  m 


2 


h  is  expressed  in  units  of  J/(s 


K) 


228 


TABLE 

i  E-12: 

LIQUID 

COOLING  - 

h  AS  A 

FUNCTION 

OF  u. 

(Ta  =  263 

K) 

TIME 

RUN 

1 

RUN 

2 

RUN 

3 

(min ) 

UN  +  i 

h 

UN  +  i 

h 

UN  *  i 

h 

0 

1.0 

- 

1.0 

- 

1.0 

- 

0.167 

0.550 

497 

0.596 

430 

0.538 

515 

0.5 

0.222 

378 

0.268 

332 

0.229 

355 

1 

0.063 

348 

0.085 

320 

0.068 

337 

1.5 

0.023 

358 

0.025 

335 

0.025 

272 

2 

0.009 

272 

0.014 

163 

0.009 

305 

2.5 

— 

— 

0.003 

447 

— 

- 

m 


h  is  expressed  in  units  of  J/(s 


K) 


22.  APPENDIX  F 


22.1  Karlsruhe  Test  Substance 

Karlsruhe  Test  Substance  is  a  methy lcellulose  gel  which 
was  first  produced  by  Riedel  (175)  in  Karlsruhe,  Germany. 
This  gel  also  contains  77%  water,  a  preservative  called 
parachlorometac resol ,  and  enough  common  salt  to  lower  the 
initial  freezing  point  to  272.2  K.  It  is  extremely  useful  in 
testing  a  freezing  model  because: 

1.  It's  enthalpy  and  thermal  conductivity  have  been 
measured  as  a  function  of  temperature  in  the  freezing 
range . 

2.  It's  thermal  properties  are  less  variable  than  those  of 
a  food. 

3.  It  behaves  very  similarly  to  foods  (especially  lean 
beef)  during  freezing. 

Riedel  tabulated  (175)  enthalpy  versus  temperature 
data.  For  ease  of  use  in  the  numerical  model  these  values 
were  approximated  by  the  following  sets  of  equations: 

H ( T )  =  285.44  +  4 . 272 (T-273 . 2 )  T  >  272.2  K 

(F-l) 

H(T)  =  exp[6. 09566  +  0 . 51751 (T-273 . 2 ) 

+  0 . 06387 (T-273 . 2 ) 2  +  0 . 003089 (T-273 . 2 ) 3 ] 

264.2  K  <  T  <  272.2  K 

( F  —  2  ) 

H ( T )  =  exp[ 4 . 9287  +  0 . 0632 ( T-273 . 2 ) 3  T  <  264.2  K 

( F  —  3  ) 
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The  reason  for  choosing  these  particular  functions  is 
explained  in  the  chapter  on  "Determination  of  Thermal 
Properties  of  Sugarbeets".  The  above  approximations  lead  to 
a  maximum  deviation  from  Riedel's  tabulated  values  of  2.6%. 
Also,  the  enthalpy  calculated  from  the  equations  agreed  with 
the  actual  data  at  293.2  K  and  263.2  K.  Thus,  the 
approximations  should  have  a  negligible  effect  on  predicted 
freezing  times. 

Bonacina  and  Comini  (22)  published  a  graph  displaying 
the  dependence  of  thermal  conductivity  on  temperature.  The 
graph  was  approximated  in  this  thesis  by  the  following 
equations : 


k(T) 

k(T) 

k(T) 


0.495  +  0.0030(T-273.2)  T  >  272.2  R 

(F-4) 

0.1312  -  0.36456T  270.2  K  <  T  <  272.2  K 

(F-5) 

0.95271  -  0 . 10682 (T-273 . 2 )  -  0 . 00572 ( T-273 . 2 ) 2 
-  0.000098(T-273.2) 3 

T  <  270.2  K 
(F-6) 


These  approximations  lead  to  thermal  conductivity  values 
within  the  accuracy  (5%)  of  reading  data  from  the  published 
graph.  The  constants  in  Equation  (F-6)  were  obtained  by 
means  of  a  least  squares  fit  through  data  read  from 
Bonacina’ s  graph.  The  constants  in  Equation  (F-4)  were 


. 


231 


determined  by  making  k(T)  equal  to  the  published  value  at 
293.2  K  and  272.2  K.  Similarly,  the  constants  for  Equation 
(F-5)  were  obtained  by  making  k(T)  equal  to  the  published 
value  at  272.2  K  and  the  value  calculated  from  Equation 
(F-6)  at  270.2  K.  This  procedure  ensured  that  the  k(T) 
function  was  continuous  at  both  272.2  K  and  270.2  K. 


23.  APPENDIX  G 


23.1  Results  for  Freezing  and  Thawing  of  Sugarbeet s 


TABLE  G-l:  FREEZING  OF  SUGARBEET  CYLINDER,  T  =243  K 

a 


TIME  CENTRE  TEMPERATURE 


(min)  (K) 


RUN  1 

RUN  2 

RUN  3 

MODEL 

0 

297.6 

298.2 

298.2 

298.0 

0.5 

287.7 

285.9 

287.5 

288.3 

1 

274.9 

272.7 

272.9 

275.8 

1.5 

270.5 

271.1 

270.5 

271.2 

2 

267.5 

264.9 

268.4 

271.0 

2.5 

257.1 

249.9 

255.2 

263.3 

3 

247.9 

245.4 

247.4 

253.7 

3.5 

244.3 

243.8 

244.5 

247.2 

4 

243.5 

243.2 

243.6 

244.2 

4.5 

243.2 

243.1 

243.2 

243.3 

5 

243.1 

243.0 

- 

243.1 

A  constant  h 

was  used  in 

value  of 

the  model 

283  J/(s 

run . 

.  m2  .  K 

232 


' 


1 
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TABLE  G-2:  FREEZING  OF  SUGARBEET  CYLINDER ,  T  =253  K 

a 


TIME 

CENTRE  TEMPERATURE 

(min ) 

RUN  1 

RUN  2 

(K) 

RUN  3 

RUN  4 

MODEL 

0 

297.2 

295.9 

298.2 

298.7 

297.5 

0.5 

289.2 

289.1 

289.7 

291.3 

288.8 

1 

276.9 

276.9 

277.2 

278.1 

276.9 

1.5 

270.5 

270.4 

271.1 

271.1 

271.9 

2 

270.2 

269.7 

270.3 

270.2 

271.0 

2.5 

269.2 

268.2 

269.1 

269.1 

271.0 

3 

267.9 

267.2 

266.9 

268.3 

269.5 

3.5 

266.5 

265.6 

265.6 

267.5 

262.7 

4 

263.7 

262.9 

262.7 

264.9 

258.8 

4.5 

259.4 

260.4 

259.1 

261.1 

256.2 

5 

256.2 

257.7 

256.2 

257.7 

254.6 

5.5 

254.7 

255.9 

254.7 

255.5 

253.7 

6 

253.7 

254.9 

253.7 

254.7 

253.3 

6.5 

- 

254.2 

253.5 

253.8 

253.1 

7 

- 

253.7 

253.2 

253.3 

253.0 

7.5 

- 

253.5 

- 

- 

252.9 

A  constant  h 

value  of 

317  J/(s  . 

m2  .  K) 

was 

used  in  ■ 

the  model 

run . 

• 

< 
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TABLE  G-3:  FREEZING  OF  SUGARBEET  CYLINDER ,  T  =263  K 

'  a 


TIME 

CENTRE  TEMPERATURE 

(min ) 

(K) 

RUN  1 

RUN  2 

RUN  3 

MODEL 

0 

299.3 

298.9 

298.2 

298.8 

in 

• 

o 

291.9 

291.1 

290.9 

290.6 

1 

280.3 

279.7 

280.5 

279.0 

1.5 

273.5 

272.9 

274.1 

273.4 

2 

270.2 

270.2 

270.5 

271.4 

2.5 

270.2 

270.5 

270.3 

271.0 

3 

269.4 

269.3 

269.7 

271.0 

3.5 

269.2 

269.1 

269.3 

271.0 

4 

269.0 

268.9 

269.1 

270.9 

4.5 

269.0 

268.9 

268.9 

270.3 

5 

269.0 

268.9 

268.9 

269.1 

5.5 

269.0 

268.8 

268.9 

267.7 

6 

268.9 

268.7 

268.5 

266.6 

6.5 

268.5 

268.4 

268.2 

265.5 

7 

267.9 

267.9 

267.7 

264.7 

7.5 

267.5 

267.5 

267.1 

264.2 

8 

266.9 

267.1 

266.5 

263.8 

8.5 

266.2 

266.4 

265.8 

263.6 

9 

265.3 

265.7 

265.2 

263.5 

9.5 

265.1 

265.3 

264.7 

263.4 

* 


' 

' 

’  :• 

MODEL 


RUN  1 

10 

264.5 

10.5 

264.3 

11 

264.1 

11.5 

263.9 

12 

- 

12.5 

— 

RUN  2 

RUN  3 

264.7 

264.2 

264.4 

264.1 

264.1 

263.8 

264.1 

— 

263.7 

— 

263.6 

— 

A  constant  h  value  of  417  j/(s  . 
was  used  in  the  model  run. 


263.3 


m2  .  K) 


■ 
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TABLE 

G-4:  THAWING  OF  SUGARBEET  CYLINDER, 

T o  =  24 3  K 

TIME 

CENTRE  TEMPERATURE 

(min ) 

(K) 

RUN  1 

RUN  2 

RUN  3 

MODEL 

(T 

=297.1) 

a 

(T  =297.1) 

a 

(T  =296.9) 

a 

(T  =297 

a 

0 

243.0 

243.0 

243.2 

243.0 

2 

256.8 

255.9 

256.1 

256.1 

4 

263.7 

262.5 

262.7 

261.9 

8 

268.1 

267.7 

267.3 

267.1 

12 

269.5 

269.3 

269.1 

268.9 

16 

270.2 

269.9 

269.8 

269.8 

20 

270.5 

270.3 

278.7 

270.3 

24 

278.7 

277.8 

284.3 

274.0 

28 

284.1 

283.5 

286.9 

280.9 

32 

286.2 

285.9 

287.9 

284.9 

36 

287.3 

287.3 

288.6 

287.5 

40 

287.9 

288.3 

289.1 

289.2 

44 

288.5 

288.8 

289.3 

290.5 

48 

288.8 

289.5 

289.8 

291.4 

52 

289.2 

289.9 

290.1 

292.1 

56 

289.7 

290.5 

290.6 

292.7 

60 

290.1 

- 

291.1 

293.1 

64 

290.6 

— 

291.5 

293.5 

68 

290.9 

- 

291.9 

293.8 

■ 
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RUN  1 

RUN  2 

RUN  3 

MODEL 

72 

291.2 

- 

292.3 

294.1 

76 

291.5 

- 

- 

— 

80 

291.8 

_ 

— 

294.5 

Heat  Transfer  Coefficient  Used  in  Model 


h  =  0.8953  +  49.3939u  +1  -  30 . 5409 ( u+ x ) 2 

N  N 

J/(s  .  m2  .  K) 


• 
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TABLE  G 

i-5:  THAWING  OF  SUGARBEET  CYLINDER, 

T o  =253  K 

TIME 

CENTRE  TEMPERATURE 

(min ) 

(K) 

RUN  1 

RUN  2 

RUN  3 

MODEL 

(T 

a 

=299.0) 

(T  =297.9) 

a 

(T  =297.9) 

3. 

( T  =298 

d 

0 

252.9 

252.9 

252.9 

252.9 

2 

260.4 

259.7 

259.9 

259.7 

4 

264.5 

264.1 

264.1 

263.9 

8 

267.6 

267.2 

267.3 

267.7 

12 

269.3 

268.8 

268.7 

269.2 

16 

269.8 

269.5 

269.7 

269.9 

20 

270.3 

274.9 

271.2 

270.3 

24 

279.4 

281.3 

279.7 

275.9 

28 

283.7 

284.3 

283.6 

282.1 

32 

285.5 

285.6 

285.9 

285.8 

36 

287.4 

286.7 

287.7 

288.2 

40 

289.3 

287.8 

289.5 

289.9 

44 

290.8 

288.7 

289.8 

291.1 

48 

292.1 

289.6 

290.7 

292.1 

52 

293.2 

290.2 

291.6 

292.8 

56 

294.1 

290.8 

292.2 

293.4 

60 

294.7 

291.4 

292.7 

293.9 

64 

295.2 

291.9 

293.3 

294.3 

68 

295.5 

292.3 

293.7 

294.6 

. 

• 
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RUN  1 

RUN  2 

RUN  3 

MODEL 

72 

295.9 

292.9 

294.1 

294.9 

76 

296.4 

- 

294.5 

295.1 

80 

— 

— 

294.6 

295.3 

Heat  Transfer  Coefficient  Used  in  Model 


h  =  0.0363  +  39.0551un+1  -  19 . 3128 ( uN + x ) 2 

J/( s  .  m2  .  K ) 
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TABLE 

G-6:  THAWING  OF  SUGARBEET  CYLINDER, 

T o  =  263  R 

TIME 

CENTRE  TEMPERATURE 

(min ) 

(K) 

RUN  1 

RUN  2 

RUN  3 

MODEL 

(T 

=297.9) 

a 

(T  =297.9) 

a 

(T  =297.7) 

a 

(T  =297 

a 

0 

263.2 

263.2 

263.2 

263.2 

2 

265.3 

265.5 

265.5 

265.9 

4 

267.1 

267.2 

267.2 

267.4 

8 

268.9 

268.9 

268.9 

269.0 

12 

269.7 

269.7 

269.7 

269.9 

16 

270.2 

270.4 

272.2 

270.3 

20 

280.3 

280.5 

281.5 

275.9 

24 

284.3 

284.5 

284.7 

282.1 

28 

285.9 

286.5 

286.2 

285.8 

32 

287.7 

287.5 

287.5 

288.3 

36 

289.1 

288.7 

289.5 

290.0 

40 

290.2 

289.8 

290.1 

291.3 

44 

291.1 

290.9 

290.7 

292.3 

48 

291.8 

291.7 

291.5 

293.1 

52 

292.3 

292.3 

291.9 

293.7 

56 

292.8 

293.1 

292.3 

294.2 

60 

293.3 

293.8 

292.7 

294.6 

64 

293.6 

293.9 

293.3 

294.9 

68 

293.9 

294.2 

293.6 

295.2 

72 

294.3 

294.4 

293.9 

295.4 

' 

< 

241 


RUN  1 

RUN  2 

RUN  3 

MODEL 

76 

294.5 

294.5 

294.1 

296.1 

80 

294.7 

— 

294.5 

296.3 

Heat  Transfer  Coefficient  Used  in  Model 


h  =  1.9937  +  26.9908u  +a  -  9.1479(u  +  1)2 

N  N 

J/(s  .  m2  .  K) 


24.  APPENDIX  H 


24 . I  Computer  Programs 


24.1.1  Model 


for  Freezing  and  Tha 


ing  of  Sugarbeet s 
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TA=INITIAL  AMBIENT  TEMP 

THERMAL  CONDUCT IVITY=C5+C6*T+C7*T**2+C8*T**3 
THERMAL  CONDUCT  I V I T Y  =  P  +  P 1*U  +  P2*U**2  +  P3*U**3 
TIME=(JA+1 )+DTT=TIME  IN  MINUTES 
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